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Abstract 

The ultimate goal for this NASA/USRA-sponsored “Apollo Lightcraft Project” 
is to develop a revolutionary manned launch vehicle technology which can poten- 
tially reduce payload transport costs by a factor of 1000 below the Space Shuttle 
Orbiter. The Rensselaer design team proposes to utilize advanced, highly energetic, 
beamed-energy sources (laser, microwave) and innovative combined-cycle (airbreath- 
ing/rocket) engines to accomplish this goal. 

The research effort focuses on the concept of a 100 MW-class, laser-boosted 
Lightcraft Technology Demonstrator (LTD) drone. The preliminary conceptual de- 
sign of this 1.4 meter diameter microspacecraft involved an analytical performance 
analysis of the transatmospheric engine in its two modes of operation (including an 
assessment of propellant and tankage requirements), and a detailed design of internal 
structure and external aeroshell configuration. The central theme of this advanced 
propulsion research was to pick a known excellent working fluid (i.e., air or LN 2 ), 
and then to design a combined-cycle engine concept around it. Also J a structural 
vibration analysis was performed on the annular shroud pulsejet engine. Finally, the 
sensor satellite mission was examined to identify the requisite subsystem hardware: 
e.g., electrical power supply, optics and sensors, communications and attitude control 
systems. 

This is the third year of RPI’s participation in the NASA/USRA University Ad- 
vanced Design Program for Aeronautics. 
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ABSTRACT 



Chapter 1 


EXECUTIVE SUMMARY 


For this year s effort, the detailed description and performance analysis of an un- 
manned 1.4 m Lightcraft Technology Demonstrator (LTD) drone is presented. The 
novel launch system employs a 100 MW-class ground-based laser to transmit power 
directly to an advanced combined-cycle engine that propels the 120 kg LTD to or- 
bit — with a mass ratio of two. The single-stage- to-orbit (SSTO) LTD machine 
then becomes an autonomous sensor satellite that can deliver precise, high quality 
information typical of today’s large orbital platforms. 

The dominant motivation behind this study is to provide an example of how 
laser propulsion and its low launch costs can induce a comparable order-of-magnitude 
reduction in sensor satellite packaging costs. The issue is simply one of production 
technology for future, survivable SSTO aerospace vehicles that intimately share both 
laser propulsion engine and satellite functional hardware. A mass production cost 
goal of $10 3 /kg for the LTD vehicle is probably realizable. 


1.1 Introduction 

In order for laser propulsion to enable a significant reduction in the cost of certain 
critical space systems, both launch and payload costs must be reduced by an order of 
magnitude or two. Canavan [1] was first to bring this fact to light, noting that a reduc- 
tion in either category alone would have much less economic impact. This conclusion 
emerged from recent in-depth examinations of the economics for laser propulsion de- 
ployment of sensors, interceptors, and decoys [1,2]. Furthermore, Canavan affirms 
that the minimum effective system must be able to launch 60-100 kg payloads. After 
evaluating cost projections, he concludes that a system designed for payloads smaller 
than this could increase costs significantly, reducing laser propulsion’s margin with 
respect to conventional chemical rocket alternatives [1], 

Apparently, the true costs of building and launching today’s large satellite plat- 
forms are not widely known, as discussed in a recent Aerospace America article [3]. For 
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example, sensor hardware can cost upwards of $200, 000/kg ($10 5 /lb); and to boost 
that sensor into geostationary orbit (GEO) typically requires $10, 000/kg ($4, 500/lb) 
of payload. Flying the Shuttle Orbiter to a 250 km low Earth orbit (LEO) may cost 
$6, 600/kg ($3000/lb), but the actual price is really almost twice that amount because 
the shuttle itself is not amortized against the payload [4]. The projected launch cost 
goal for the Advanced Launch System (ALS) is roughly $660/kg ($300/lb) to LEO. 

Hence, for laser propulsion to play a significant role in boosting future critical space 
systems, launch costs must fall to $66/kg ($30 /lb), or at least below $100/kg. 

Canavan has raised the fascinating issue of whether and how laser propulsion 
and its low launch costs could induce a substantial reduction in satellite package 
costs [l]. In the opinion of the authors, such reductions could be facilitated by an 
exceptionally close integration of the laser-propulsive engine and satellite functional « 

hardware. Pushed to the extreme, alm ost every vehicle component could be designed 
to serve multiple functions, in both transatmospheric and orbital flight modes. 

Clearly, the final configuration of any laser-boosted machine will be strongly 
driven, if not entirely dominated by the mission it must perform, be it intercep- 
tor, decoy or sensor. A near infinite number of successful configurations could be 
alleged to exist, but it is most instructive to select a specific mission , and then to 
explore a single configuration from the initial design concept, through the preliminary 
engineering design process. 


1.2 Initial LTD Design Concept 


The advanced aerospace vehicle considered here is exemplary of a class of sensor ma- 
chines that can be derived largely from an intimate integration of propulsion and 
sensor systems. The proposed design exploits the inherent advantages of advanced 
beamed-energy sources (i.e., high power lasers) and innovative combined-cycle (air- 
breathing/rocket) engines to accomplish this goal. The authors believe that this 
unique approach could possibly enable a reduction in both launch and sensor package 
costs by two orders of magnitude below present levels. However, as pointed out by 
Canavan [1], the numbers of these sensor satellites may not be great enough to justify 
the expense of the entire laser launch facility for this application alone. Nevertheless, 
the laser launch facility is likely to be built for a completely different defense-related 
purpose, and amortized over a great number of users and dissimilar mission applica- 
tions anyway. 
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Figure 1.1: Ground-based FEL laser launch facility. 

1.3 Ground Based Laser Launch Facility 

Portrayed in Figure 1.1 is a 100 MW-class ground-based laser (GBL) that could be 
built with free-electron laser technology in the next 5 years, by assembling numerous 
(e.g., 5) smaller units into an array. Redundant units could be built into the system so 
that inoperative modules could be dropped out, with no loss in system utility during 
a boost. As shown in Figure 1.1, a beam combiner could then be invoked to link all 
the output beams together. 

All the units could be fired simultaneously to give the lowest pulse repetition 
frequency (PRF) of perhaps several hundred Hertz at the highest pulse energy (E p ); 
at the other extreme, each unit could be triggered sequentially to yield the lowest E p ’s 
and PRF’s up to 10 kHz. Near-term pulsed laser propulsion engines can be designed 
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around whatever laser pulse durations are available (e.g., from 40 to 400 nsec.). 

Configured in this manner, the GBL facility can be programmed to deliver a 
complicated pulse train sequence of PRF’s, E p s and t p ’s with the utmost of ease. This 
pulse train can be calculated to exactly match what a laser propulsion engine will 
need along a given launch trajectory — i.e., normally a direct function of flight Mach 
number and altitude. The goal is to produce an efficient thruster without introducing 
too much flight hardware, which has added so much to the cost of chemical rockets, 
as Kantrowitz [4] has observed. 

The essential point of this advanced launch scheme is to place as much of the laser 
power transmission system complexity as possible on the ground (where there is no 
weight penalty) so that it can be serviced easily. With this approach, laser powered 
thrusters can be reduced to their simplest, and most reliable configuration. 

The “straw man” GBL facility suggested here is set at the 250 MW level (peak) 
which is adequate to launch a 120 kg (dry mass), 1.4 meter diameter Lightcraft 
Technology Demonstrator (LTD) to low Earth orbit. The range of laser pulse energies 
required by its combined- cycle (airbreathing/rocket) engine is 40 KJ < E p < 70 KJ; 
PRF varies from 200 Hz to 10 kHz, and t p varies from 0.3 to 0.4 fisec — depending 
on the exact trajectory '(i.e., Mach number vs. altitude) flown to orbit. With these 
parameters, the peak flux across the 1.0 m diameter LTD primary optic will fall in 
the range of 13.0 to 30.0 MW/cm 2 . 

Kantrowitz [5] notes that the important costs for a GBL installation, are for 
capital and operating expenses (which might add another 20% of the capital cost per 
year). Refer to Ref. 5 for an in-depth accounting of the economics for a GBL launch 
facility. The Free Electron Laser (FEL) is the favored future GBL system due to its 
promise of high electric-to-laser conversion efficiency and reliability. 

Kantrowitz [5] has noted that to make laser propulsion a serious contender for 
space transportation to low Earth orbit (LEO), it is necessary to develop propellants 
that can achieve high thruster efficiency at low incident laser flux levels. It is also 
apparent that atmospheric transmission problems must be considered, especially in 
the immediate vicinity of the vehicle where the beam must propagate unhampered 
through the thruster’s rapidly expanding, and potentially absorbing exhaust. 

Adaptive transmitter optics can be invoked to successfully bring the power beam 
up through the atmosphere. The 10 m diameter beam-director mirror would allow a 
10 ftm beam to be focused on a one meter diameter vehicle base, out to a range of 
about 800 km. This performance is, of course, close to the diffraction limi t. With 
even shorter wavelengths (e.g., 1 /im), the transmitter diameter could be reduced to 
4 m (and below). 
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Figure *1.2: Lightcraft Technology Demonstrator (1.4 meter diameter). 

1.4 The LTD Combined-Cycle Engine Concept 

The Lightcraft Technology Demonstrator (LTD) is a single-stage- to-orbit transatmo- 
spheric vehicle ut ilisi ng both airbreathing and rocket propulsion modes. The LTD 
will transition into the rocket mode at about Mach 5 and 30 km altitude. 

The LTD advanced microspacecraft is in several respects literally a “flying engine.” 
The forebody aeroshell acts as an external compression surface: e.g., the airbreathing 
engine inlet. The afterbody has a dual function: it is the primary receptive optic 
(parabolic mirror) for the laser beam which provides power to the engine, and it is 
also an external expansion surface (plug nozzle) during the rocket engine mode. It 
is the opinion of the authors that focusing mirrors mounted on the laser-propelled 
vehicle are the most expedient way to permit both low flux levels in the atmospheric 
transmission link and elevated flux levels within the thruster ( which are required for 
high propulsive efficiency). The reference “point design” for the LTD is configured 
around a 1.0 meter diameter parabolic receiving mirror, as shown in Figure 1.2. 

The primary thrust structure is the annular shroud. The shroud serves as both 
inlet and impulsive thrust surface during the airbreathing engine mode; at this time, 
the LTD centerbody is treated as the “vehicle airframe” component. In the rocket 
mode, the annular inlet is closed; the afterbody and shroud combine to form the rocket 
thrust chamber. The three primary structures: forebody, shroud, and afterbody, cire 
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interconnected by an internal support frame to which all internal subsystems are 
attached. 


1.5 Focus of the LTD Design Project 

This year’s efforts have been to develop a conceptual design for a small, 1.4 meter 
diameter Lightcraft Technology Demonstrator around components derived from cur- 
rent liquid propellant chemical rocket engines, advanced composite structures, and 
high power laser mirrors. Specific areas that were addressed by the design team in- 
cluded aerodynamics, propulsion, structures, propellant management, heat transfer, 
optics for sensor and propulsion missions, lightsat/microsat subsystems, mechanical 
subsystems, detailed LTD mockup construction, and finally, the definition of future 
experimental test facilities (for next year’s effort). 


1.5.1 Aerodynamics and Propulsion 

In the area of aerodynamics and propulsion, an analytical computer model of the 
airbreathing pulsejet cycle was first developed, providing a data base of engine per- 
formance characteristics. A vehicle aerodynamic drag model was then assembled. 
Next, an optimal trajectory analysis was then performed using a computer code called 
SORT (Simulation and Optimization of Rocket Trajectories), which had previously 
been modified for light craft engines and aerodynamics. 

The focus of the propulsion analysis was to model the laser-heated blast wave 
and impulse- generation process using an axi-symmetric representation, that excludes 
radiation and convection. For the most part, air within the laser- generated blast wave 
was modeled as an ideal gas, throughout the expansion process. The engine inlet air 
state (i.e., pressure, temperature, density, etc.) prior to laser- heating, was a direct 
function of flight Mach number and altitude; hence, graphical results are presented 
versus these flight variables. The equations developed by Raizer [6,7] to describe the 
initial state of the laser-induced blast waves, were used in modeling the “line-source,” 
impulse-generation process. 

The simple LTD inlet flow model predicted the state of the inlet air (pressure, 
temperature, velocity, etc.) which refreshes the lower annular shroud impulse sur- 
face. This aerodynamic model was needed not only for the refresh air state, but also 
for defining external drag characteristics of the entire vehicle. In the latter, it is impor- 
tant to distinguish between the drag produced by engine-related vs. airframe-related 
components, so as not to over penalize Lightcraft performance. 
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1.5.2 Trajectory Analysis 

The trajectory analysis of any launch vehicle becomes a critical step in the overall 
system integration process. Many important engine/ vehicle related characteristics 
must come together for the final product: performance. Launch vehicle performance 
is typically measured by payload capability, which for the LTD is the entire dry mass 
and ullage gas, totaling 124 kilograms. The LTD must attain low Earth orbit with the 
available propellant, while minimizing total laser energy consumed along the launch 
trajectory. 

The trajectory was evaluated using a computer tool called SORT which was writ- 
ten by McDonnell Douglas Astronautics Co. for NASA, to design space shuttle trajec- 
tories [8]. The code is sufficiently general such that it can analyze any trajectory and 
model all important environmental parameters that affect vehicle dynamics. Engine 
performance, vehicle aerodynamics, guidance algorithms, and mass histories interact 
with atmosphere and gravity models. 

The capability of the SORT program is quite extensive in its use of these sophis- 
ticated vehicle and environmental models. The program can iterate on trajectory 
parameters to optimize performance, achieve a desired criteria, or constrain the solu- 
tion to avoid some specified limit. 

Even with all of the generality built into SORT certain modifications were required 
for the LTD. The most significant software modification involved the unique energy 
source, a laser. A new vehicle steering option was also encoded so that the LTD could 
always point at the ground-based laser station (or laser relay satellite). After these 
modifications were included, SORT was able to model Lightcraft performance to a 
high degree of accuracy. 

1.5.3 Structural Design 

Once the overall LTD vehicle structural configuration was defined, more detailed 
analysis of primary and secondary structure subcomponents was carried out. Consid- 
eration was given to using state-of-the-art materials technology and mass production 
techniques for low unit cost (e.g., S 10 5 /vehicle) . A close integration of laser-propulsive 
engine and satellite functional hardware was necessary. Much emphasis was placed on 
the annular shroud structure — the primary impulsive thrust surface (see Fig. 1.3). A 
finite element computer code called CAEDS™ (Computer Aided Engineering Design 
System) was used to conduct both static and dynamic structural analyses; stress and 
modal behavior of the LTD shroud structure were extracted. This work will continue 
in the future as the LTD structured model becomes more elaborate. 

Driving issues in the LTD structural design process were the component mass al- 
lowances, and propulsion system requirements; aerothermal loads were also considered 
in a generic sense, to identify the current operating limits. Specific transatmospheric 
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Figure 1.3: The LTD Shroud and support Stmt. 


flight paths must be analyzed in future studies to assess their suitability for alternative 
engine and structures technology. 

The LTD structure must meet the combined challenge of providing a lightweight, 
efficient airframe that can survive severe engine and aerodynamic heating: i.e., a 
“ thermo- structural” viewpoint. Clearly, any successful thermo-structural concept 
must meet existing materials and manufacturing limitations. 

The LTD primary structure is a large fraction of the vehicle’s inert weight; yet, 
this comprises a small structural weight fraction. Applied loads must be predicted to 
an accuracy that will allow knowledgable reductions in margins of safety, to eliminate 
every ounce of non-optimum structure. A large experimental data base of aerothermal 
loads must now be assembled, for future detailed finite element structural and heat 
transfer analyses. 


1.5,4 Structural Dynamic Analysis 

Finite element modeling techniques were used to evaluate the annular shroud and sup- 
port strut design proposed for the LTD (see Fig. 1.4). Both static and dynamic struc- 
tural analyses were implemented to extract the stress and modal behavior. Pulsed 
airbreathing or rocket engines must avoid the natural frequencies of the vehicle struc- 
ture, or catastrophic failure of the system will result. The present analysis concluded 
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Figure 1.4: Finite element models of Shroud and Strut. 


that the LTD shroud and support structure is able to withstand stresses of induced 
oscillations at the current engine design P RF s of 1 to 10 kHz. Even though resonant 
frequencies may be produced in flight, extremely small amplitudes (of displacement) 
will not cause damaging stress levels. The LTD structure seems to resonate at fre- 
quencies in the 10 to 100 Hz range. 

The static stress analysis examined the integrity of the shroud/strut section, and 
revealed areas that needed further refinement. The shroud itself was found to be 
excessively strong and can easily sustain the applied loads; thus some material can 
be removed from the shroud interior to permit further weight reduction. Also, the 
proposed strut mounting design provides sufficient strength for attaching the struts 
securely to the internal LTD support frame. 

1.5.5 Propellant Management Systems Design 

The LTD combined- cycle, laser-propulsion system employs onboard cryogenic fluids 
for reaction propellant, sacrificial coolant and feed system pressurant. Liquid nitrogen 
has been chosen as the propellant for several reasons. Although liquid hydrogen would 
be the most desirable propellant because of its low molecular weight (and higher rocket 
specific impulse), only 14 kg would fit into the LTD’s 28 in. diameter storage tank; 
in contrast, 140 kg of LN 2 would occupy the same volume (see Fig. 1.5). Heated to 
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Figure 1.5: Cutaway view of LTD showing storage tanks for LN 2 and Helium. 

high temperature and pressure, nitrogen can produce specific impulses in the range 
of 725 to 1025 seconds — depending on whether the exhaust gases are dissociated 
or not. (LN 2 is also favored because it is an inert and exceptionally clean coolant, a 
consideration which is especially important for high power laser optics.) 

Several candidate injector designs have been identified for the rocket engine mode 
(see Fig. 1.6). The injector orifices are designed to act as “acoustic- valves,” to auto- 
matically deliver propellant into the rocket engine thrust chamber. Precise control of 
injector orifice pressure is necessary for optimum propellant utilization during flight. 

High pressure, super- critical helium was chosen for the pressurant because of its 
superiority as an inert agent with a very low boiling point. The LTD pressurant 
system is modeled after one flown on the Apollo Lunar Module descent rockets. 
Storing the helium at near liquid temperatures enables maximum pressurant density. 

The computer- controlled propellant management system must: 1) deliver the 

proper mass flow rate of liquid nitrogen coolant to critical actively- cooled engine 
components during the airbreathing propulsive mode; and 2) provide uniform propel- 
lant mass flow during the rocket mode. Sensors distributed over the vehicle’s outer 
skin and engine will monitor surface temperatures. The mass flow rate of LN 2 coolant 
will be adjusted to compensate for variations in laser power absorbed by the primary 
mirror during the airbreathing mode. Laser power input is constant during the rocket 
mode. 
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Figure 1.6: Proposed LTD Injector design 


The heat transfer analysis has focused upon active cooling of the LTD’s engine hot 
structures: i.e., mirror, shroud and shroud support struts. Regenerative cooling will 
be employed by the primary mirror/plug nozzle and hot thruster impulse surfaces. 
The leading edges of the shroud support struts and the shroud will utilize transpi- 
ration cooling, and platelet technology. The shroud thruster surface and propellant 
injector array will be film cooled. A thermal protection system (TPS) such as used 
on the Space Shuttle Orbiter will cover the vehicle forebody. 

1.5.6 Optics for Propulsion 

The design and fabrication of the LTD primary receptive mirror has been investigated. 
Although some of the physics behind propagating the high energy laser beam through 
the heated exhaust plume is not fully understood, the concept is shown to be within 
reach of current technology. The primary optic could be manufactured using existing 
diamond- turning lathe facilities. Development of multi-layer high reflectivity coatings 
and microchannel cooling systems are the main technical problems encountered. The 
mirror must exhibit high figure quality and minimal surface roughness throughout 
the launch trajectory. Calculations of beam divergance indicate that a 4 m di am eter 
adaptive ground-based transmitter is entirely sufficient for a 1.0 micron wavelength. 
The effects of pointing error were evaluated using a computer ray trace routine; errors 
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of 0.5 degrees or less shift the focus asymmetricadly on either side of the annular 
engine, but the effect is self- correcting and the LTD is therefore stable in this respect. 
At incident angles up to 1 degree, the power loss and spot size increase but do not 
cause a critical loss of thrust for near, intermediate, or far field beam profiles. An 
adaptive pointing error correction system is proposed, using retroreflectors to sample 
the power beam edges (at 24 azimuthal locations) and feed data back to the ground 
power source via laser rangefinding. 

1.5.7 Optics for Sensor Satellite Mission 

The use of the LTD as a sensor satellite has been explored. The primary mirror could 
serve a dual purpose on each mission; however, since the mirror is designed primarily 
for propulsion applications, its performance in other respects is limited. A ray trace 
analysis of the primary optic indicates that it has superior light gathering ability, but 
the angular field of view is hmited to only about 0.5 degrees. Coma is the dominant 
aberration over the entire field of view; a small amount of astigmatism is present for 
larger angles of incidence, while other aberrations are essentially negligible. Several 
candidate sensors, including CCD and CID technology, are commercially available 
and meet the requirements of the LTD retina. 

A special mechanical subsystem was developed to deploy this advanced segmented 
photo-optic sensor (the “ring retina”), into the focal region of the primary mirror; in 
this manner, the propulsive optics are transformed into a powerful one meter diameter 
telescope. Functional hardware necessary for the LTD to perform satellite functions 
is located in the forward (nose) payload section of the vehicle. 


1.5.8 Lightsat Systems Design 

Once in orbit, the LTD will function as a light-weight sensor satellite, complete with 
all the requisite on-board electronics, communications, and attitude control systems 
(see Fig. 1.7). The LTD senses its attitude with a sun sensor and a Earth horizon sen- 
sor. To change attitude, it uses a combination of cold gas jets and/or magnetic torque 
bars. It will be equipped with an erectable light-weight antenna for communicating 
with launch co mman d, accessing the TRDSS (and STDN) network, and relaying data 
gathered from LTD sensors. An on board computer will manage satellite electrical, 
mechanical and attitude control systems, coolant flow rates, and solar array actua- 
tors. The solar array is mounted onto the inner forebody surface, which opens into 
four “petals” to collect solar energy for recharging on-board nickel-hydrogen batter- 
ies. The batteries are used primarily when the satellite traverses the Earth’s shadow. 
This combination of satellite components is chosen for minimal cost and weight, so 
that the LTD can fulfill its primary objective of a self-launched lightsat. 
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Figure 1.7: LTD satellite functional hardware is located in the forward payload section 
of the vehicle. 


1.6 Future Directions 

In the upcoming year, efforts will be focused on the LTD’s propulsion system and vehi- 
cle aerodynamics rather than on its mission as a sensor satellite. A 3-D computational 
fluid dynamic (CFD) model of the LTD’s external inlet will be assembled, and run 
for a variety of different shapes and flight speeds. From these computer simulations, 
the boundary layer thickness, external surface pressure distribution, and transition 
point (from laminar to turbulent flow) can be determined. A Naval Research Labora- 
tory (NRL) 2-D blast wave code may also be employed to model the impulsive thrust 
generated by laser-induced blast waves. In addition, a computer analysis of the radia- 
tion/convection heat transfer to the LTD will be conducted to find the engine/ vehicle 
thermal profile. All of this will enable the choice of appropriate spacecraft materi- 
als (specifically tailored to withstand the hypersonic transatmospheric environment), 
and also to design active thermal cooling systems for the primary optic and engine 
hot-sections. Finally, a more detailed vibrational analysis will be performed that in- 
cludes not only the shroud/strut assembly, but also the spacecraft interned support 
structure. 

Several experiments will be performed to validate theoretical vehicle and engine 
performance simulations (see Fig. 1.8), Hot and cold flow wind tunnel tests are planned, 
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with two different scale models, 5.5 in. diameter and 1.25 in. diameter. The air 
speeds will be varied from subsonic to hypersonic. Schberen photographs and pressure 
data will be taken. Laser impulse experiments will also be performed with an exact 
1 / 15th segment of the annular engine (see Fig. 1.9), in both static and dynamic wind 
conditions. These high power laser experiments will yield engine impulse and heat 
transfer data to prove technical feasibility of the propulsion concept. 

Finally, additional systems-integration questions that were not addressed in this 
final report may be explored: e.g., LTD range safety issues, manufacturing cost anal- 
yses, and high power laser attenuation by the engine exhaust plume. 
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Figure 1.9: LTD Laser Impulse Experiment Test Section (side view). 
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Chapter 2 


VEHICLE CONFIGURATION AND 
STRUCTURAL DESIGN 


2.1 Introduction 

The purpose of this chapter is twofold: first, to present an overview of the LTD struc- 
tural design influences from the perspective of the aerothermal flight environment; 
and second, to describe the overall LTD vehicle configuration and structural design. 
It is important to note that the vehicle preliminary conceptual design documented 
in this report is a bold first step in resolving conflicting design requirements and 
integrating complex vehicle syst ems 

There are two stages in the conceptual design process. First is the creation of a 
preliminary concept, in which the functions and desired features that comprise a suc- 
cessful design are identified. Alternative concepts that also meet these requirements 
should be explored. In the second stage the most promising concept is selected for 
further analysis. The structural design of the LTD has emerged from the conceptual 
design phase and is now undergoing more detailed analysis. 

Driving issues in the LTD structural design process are the component mass al- 
lowances, and propulsion system requirements; aerothermal loads have been consid- 
ered in a generic sense, to identify the current operating limits. Specific transatmo- 
spheric flight paths must be analyzed in future studies to assess their suitability for 
LTD engine and structures technology. 


2.2 The Structural Design Process 

The LTD structure must meet the combined challenge of providing a lightweight, 
efficient airframe that can survive severe engine and aerodyn ami c heating: i.e., a 
thermo-structural” viewpoint [1], Clearly, any successful thermo-structural concept 
must meet existing materials and manufacturing limitations. Initial selection of ma- 
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terials which may be considered for application to the LTD requires that the range 
of operating temperatures and the accumulated time at maximum operating temper- 
ature be known. 

The LTD primary structure is a large fraction of the of vehicle’s inert weight; yet 
this comprises a small structural weight fraction. Applied loads must be predicted to 
an accuracy that will allow knowledgable reductions in margins of safety, to eliminate 
every ounce of non-optimum structure. A large experimental data base of aerothermal 
loads must now be assembled, for future detailed finite element structural and heat 
transfer analyses. 

A transatmospheric vehicle must withstand aerothermal, aerothermoelastic, acous- 
tic, inertial, and engine thrust loads. Aerothermal loads are high because a single- 
stage-to-orbit vehicle with airbreathing propulsion follows a high dynamic pressure 
trajectory to achieve the necessary propulsion efficiency [2]. Added to this thermal 
load is the severe radiation environment imposed on the power-beam receptive optics, 
essential to laser-propulsive engines. 

A compilation of the structures and materials technology needs for the LTD are 
represented below: 

1. Propulsion /Thrust Structures: for which rigorous advancements are required. 

2. Hot Structures: to handle high operating temperatures without degradation of 
strength or structural integrity. 

3. Cryogenic Tankage: requiring simple and lightweight insulated structures able 
to handle cryogenic temperatures and supercritical fluids. 

4. Thermal Protection Systems (TPS) : Utilizing a diverse range of materials in- 
cluding ceramics, carbon-carbon and refractory metals. 

5. Dimensionally Stable Materials: which are essential for rigid optical platforms 
(carbon fiber composite structures are ideal candidates). 

6. Active Control Surfaces: that will enable the propulsion system to transition 
smoothly between engine modes during flight. 

7. Thermal Control Coatings: to reduce the severe heating loads experienced by 
hypersonic vehicles and laser-receptive optics. 


2.3 LTD Vehicle Configuration 

The Lightcraft Technology Demonstrator (LTD) is a single-stage-to-orbit transatmo- 
spheric vehicle utilizing both airbreathing and rocket propulsion modes (see Fig. 2.1). 
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Figure 2.1: 
al) view 
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Figure 2.2: The LTD engine/ vehicle external geometry 


The LTD, an unmanned launch vehicle, 1.4 meters in diameter and 1.77 meters in 
length, will transition into a remote sensor satellite upon reaching orbit. The moti- 
vation behind this vehicle concept is to provide an example of how laser propulsion 
and its low launch costs can induce a comparable order-of-magnitude reduction in 
satellite pa cka ging costs. The LTD will serve as a test bed for near term demonstra- 
tion of laser propulsion technology. In its role as a ‘lightsat’, the LTD provides a 
low cost alternative for remote sensor satellite missions. This is accomplished by an 
exceptionally close integration of the laser propulsion engine and satellite functional 
hardware. 

The LTD advanced microspacecraft is in several respects literally a “flying en- 
gine” (see Fig. 2.2). The forebody aeroshell acts as an external compression surface 
for the airbreathing engine inlet. The afterbody has a dual function: it is the pn- 
marv receptive optic (parabolic mirror) for the laser beam which provides power to 
the engine, and it is also an external expansion surface (plug nozzle) during the rocket 
engine mode. The primary thrust structure is the annular shroud. The shroud serves 
as both inlet and impulsive thrust surface during the airbreathing engine mode. In 
the rocket mode, the annular inlet is closed; the afterbody and shroud combine to 
form the rocket thrust chamber. The three primary structures: forebody, shroud 
and afterbody, are interconnected by an internal support frame to which all internal 
subsystems are attached. 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPM 


i 


2.4. LTD FOREBODY 


23 



Figure 2.3: Forebody flaps open and close to deploy inlet shutter and retina 

2.4 LTD Forebody 

As mentioned above, the LTD forebody acts as an external compression surface. In 
addition, the forebody encloses the payload compartment, which is dedicated to satel- 
lite and launch control system electronics. The forebody is assembled from three sep- 
arate sections. The advanced caxbon-caxbon (ACC) nosecap is designed to withstand 
the highest thermal loads (at the stagnation point) experienced during hypersonic 
flight. The forebody mid-section consists of four hinged petals. Once in orbit the 
petals open and reveal photovoltaic arrays which are mounted to the inner petal sur- 
faces. Finally, the base section is comprised of twenty-four “flaps” (ACC material) 
that are hinged at the forward edges(see Fig 2.3). During a flight to orbit, they are 
opened twice: first, when the annular air inlet shutters are closed, and later after 
orbit is achieved (to permit deployment of the ring retina sensors). The entire fore- 
body is hermetically sealed to prevent infiltration from the environment during the 
transatmospheric flight. 

The forebody nosecap will not require active cooling. Advanced carbon-carbon 
composite structures have been flight tested on the Space Shuttle Orbiter and have 
proven ideal for the high temperature transatmospheric environment. Carbon-carbon 
composites utilize carbon fibers in a carbon matrix. The nosecap is a shell structure 
provided with thermal barriers to minimize heat transfer into the payload compart- 
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Figure 2.4: Side (cross sectional) view of forebody 

ment. The ACC paxts must be protected from oxidation by a high temperature 
oxidation- stable coating such as silicon carbide. On the shuttle orbiter leading edges 
ACC materials have performed well at elevated temperatures of 2500° F [3]. 

The forebody mid-section forms a composite shell structure (see Fig. 2.4) and has 
a truss-reinforced primary sub-structure made of polymide graphite (PG) laminates. 
The structural temperature limit of PG laminates is superior to graphite/ epoxy, and 
they enable a forty percent weight savings over aluminum [4]. The forebody shell 
structure will be covered by a thermal protection system (TPS), consisting of tiles 
made from boron rich fibrous refractory composite insulation (FRCI) [5]. FRCI is a 
re-radiative thermal protection system; the tiles dissipate most of the heat load before 
it can be reach the vehicle primary structure. This is accomplished by a reaction cured 
glass coating (RCG) only twelve mil s thick; it forms the forebody aerodynamic skin 
and re-radiates 90% of the incident heat energy back to the environment. These tiles 
require a flexible attachment system based on strain isolation pads, and carefully 
designed gaps at tile interfaces. The reduced heat pulse will not reach the primary 
structure until after the LTD achieves orbit. The photovoltaic array is attached 
to inner petal surfaces, and is thermally protected by radiation heat shields which 
prevent overheating from the residual heat pulse (see Fig. 2.5). 

In order to provide a rigid forebody shell during flight, the mid-section petals 
are firmly secured by a ‘Super-Zip’ separation system (see Fig. 2.6), used for many 
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SUPERZIP UNEAR ORDNANCE SEPARATING DEVICE 



Figure 2.5: Super-Zip Linear Ordnance separating device 


years in the space program [6]. Petal interfaces are temporarily fastened together by 
notched plates. Within the notched plate joint is a tube containing linear ordnance 
(an exploding cord) embedded in silicon rubber matrix. When the pyrotechnic cord is 
detonated, the resulting explosion fractures the notched plates, and a clean separation 
is produced in milliseconds. The mechanical interface must be designed to attenuate 
vibration from the explosion. 

The aerodynamic contour of the LTD forebody can be described as a doubly- 
curved ogive with blunted nose. The initial design contour for the forebody was a 
thirty degree half- angle cone atop a forty-five degree half-angle conical frustrum. The 
ogive shape (see Fig. 2.7) was chosen after a preliminary bow shock wave analysis 
(using NACA 1135) indicated that the shock would be swallowed by the annular 
engine inlet before the transition from airbreathing to rocket mode (inlet design Mach 
number, M = 5). 

The forebody aerodynamics are intimately related to the vehicle propulsion sys- 
tem. At supersonic flight speeds, the vehicle bow shock wave pre-compresses the free 
stream air for the airbreathing engine cycle; hence, the inlet is formed between the 
forebody and the shroud leading edge. The ideal geometry for this inlet would have 
the bow wave attach to the inlet lip at the design mach number; all the air compressed 
by the bow wave would thereby be ingested into the engine at Mach 5 [7]. 

The pressure and density flow fields around the annular inlet station provide 
insight into the uniformity and magnitude of the external compression process. Total 
pressure recovery, process efficiency, and kinetic efficiency provide measures of inlet 
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PHOTOVOLTAICS ARRAY 



Figure 2.6: •Side (cross-sectional) view of proposed forebody aeroshell. 

performance. Optimal design of an inlet requires sophisticated 3-D computational 
fluid dynamic (CFD) codes, as well as wind tunnel experiments to validate analytical 
results [8]. 

The forebody shape was selected to minimize vehicle pressure drag characteristics. 
The resulting contour can be described as 7i-power body: 


(y/b) = (x/a) n 


(2.1) 


where n = 2/3, and the resulting shape is an approximation of the minimum pressure 
drag body for a given fineness ratio (i.e., length/diameter) [9]. This blunted, minimum 
pressure contour should also delay attachment of the bow wave to the inlet lip until 
the design Mach number is reached. 

Obviously, several criteria could be used for optimizing the aerodynamic contour 
of the LTD forebody; pressure drag, total pressure recovery, heat transfer, etc. It is 
clear that a single optimum shape necessitates compromise between several conflicting 
measures of aerodynamic efficiency [10]. 
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Figure 2.7: Bow shock wave and forebody contour at inlet design Mach n um ber 
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Figure 2.8: T-tail shroud support strut assembly (24) 


2.5 LTD Shroud/Impulse structure 

As mentioned earlier the shroud also serves as primary structure for delivering im- 
pulse to the LTD. In the airbreathing mode the ERH thruster itself is extremely 
simple, consisting only of an impulse surface (inner, or lower shroud surface) which 
can withstand high temperatures and pressures. To generate thrust, a laser- induced 
plasma is first ignited over the thruster surface. This high pressure (600 atm) plasma 
bubble then expands as it is convected down stream across the lower shroud surface. 
The shroud structure will be the most abused part of the vehicle during powered 
flight, and must be exceedingly strong. 

The shroud structure is comprised of two major subassemblies: 1) the leading edge 
section is assembled from twenty-four nickel alloy “T-tail” support struts that also 
form the forward half of the shroud (see Fig. 2.8); and, 2) a single continuous ACC 
composite ring section forms the aft half (or trailing edge) of the shroud. The nickel T- 
tail assemblies require active cooling at the leading edges (transpiration) and thruster 
surfaces (film), in order to withstand the harsh aerothermal and laser propulsive 
environments. The rear half of the shroud is an ACC hot structure (passively cooled), 
to which the T-tails are firmly attached. 

Actively- cooled structures must be used when radiation equilibrium temperatures, 
which are dependent on vehicle configuration and flight trajectory, exceed usual ma- 
terial limits, (e.g., about 2000° F for nickel-based alloys). The proposed shroud 
cooling system includes transpiration techniques and platelet technology [11]. Each 
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T-tail support strut would be manufactured by stacking numerous thin metal plates, 
into which intricate flow passages are chemically etched (see Fig. 2.9). The etch 
patterns combine to form complex flow paths and super-fine filters (for the sacrifi- 
cial coolant) that determine the appropriate coolant flow rates to all hot areas of 
the structure. Once stacked the plates are diffusion bonded (at high temperatures 
and pressures) together, providing a high degree of structural integrity and strength 
[16]. This construction technique has been successfully demonstrated in experimental 
rocket engine programs [12], and in scramjet fuel-spray struts developed for the X-30 
National Aerospace Plane (at Aerojet Techsystems). 

As indicated above, the material favored for the rear half of the LTD shroud is 
the ACC structure which has also been used in the solid rocket booster cases for the 
Space Shuttle Orbiter. Two critical design issues for candidate materials are: thermal 
shock experienced at extremely high operating temperatures without the benefit of 
active cooling; and cyclic fatigue resulting from repetitively pulsed engine operation. 

The maximum heat-flux capability of a carbon-carbon structure is currently lim- 
ited by the temperature capability (about 2800° F) of its oxidation protection coating. 
If necessary, it will be possible to actively cool the shroud’s carbon fiber structure by 
incorporating heat pipe technology into the carbon-carbon structure. A new design 
concept developed at NASA Langley Research Center [19], which NASA intends to 
patent, combines high temperature heat pipe and carbon-carbon technologies to ex- 
tend both the thermo-structural capabilities and the maximum heat-flux capability of 
leading edge structures of hypersonic vehicles. The new design uses refractory- metal 
heat pipes embedded within a carbon-carbon structure (see Fig. 2.10). The walls of 
the heat pipes are thin and contain the working fluid (lithium or sodium) of the heat 
pipe. The primary application considered for such a structure is removal of high local 
heat from a leading edge surface, for transport to cooler surfaces, where it can be 
removed by radiation and/or active cooling. 

The shroud structure is the most massive component in the LTD vehicle. This 
hollow, stressed-skin structure is comprised of an outer, load carrying shell which 
is reinforced by internal vertical ribbing. A common example of this technology is 
actively- cooled jet engine turbine blades [13]. A critical design issue is the thermal 
expansion evidenced by such components under typical load conditions. The inter- 
faces between each hot component must accomodate both radial and circumferential 
expansions. Since the LTD lower shroud surface must coincide with the laser beam 
focus, radial shroud motion (i.e., expansion) would result in the loss of focusing pre- 
cision, and degradation of thruster performance. The shroud support struts must 
therefore be attached to the control support frame in such a manner that compen- 
sates for all thermal expansions, especially in the radial and circumferential directions. 
Differential thermal expansions will also occur at the interfaces between the nickel 
alloy “T-tail” and the carbon fiber composite rear shroud ring, possibly producing 
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Figure 2.9: Typical platelet flow metering and distribution schematic 
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Figure 2.10: Metal pipes integrated into a carbon-carbon structure provide active 
cooling (from Ref. 19). 

high sheer stresses at the contact surfaces. These interfaces may require innovative 
joining techniques and fastener technology. 

2.6 LTD Afterbody 

The Afterbody of the LTD performs two functions: 1) It is the primary receptive op- 
tic (parabolic mirror) for capture and reflective focusing of the incoming laser power 
beam; and 2) it acts as an external expansion surface (plug nozzle) for the laser-heated 
rocket engine mode. The receptive mirror function is the driving issue in the geom- 
etry and structural design of the afterbody; the mirror must be dimensionally stable 
over the entire LTD operating range in order to maintain its precise focusing ability. 
High energy laser mirror technology will be pushed to the limit, with reflectivities 
approaching the theoretical maximum values (e.g., > 99.9%) attainable by multi- 
layer film coating technology. Active cooling is essential in order to keep temperature 
gradients in the mirror substrate structure within limits; thermally-induced stresses 
produced by even moderate temperature gradients would seriously degrade the mirror 
focusing ability and therefore thruster performance. Internal mirror stresses resulting 
form high coolant operating pressures must be considered. The mirror surface must 
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also survive the pulsed thermal shock during engine operation, and will require a 
durable outer coating (e.g., vapor deposited diamond) that can maintain a superior 
surface quality. 

The afterbody structural design calls for a carbon composite inner shell sandwich 
to provide the necessary strength and dimensional stability. Since carbon- carbon is a 
composite material, its mechanical and thermophysical properties can be tailored to 
specific requirements by selection of the type of the weave, and the type and condition 
of the fibers. The optimum choice of properties of the composite will result in a 
minimum- mass design that satisfies peak-temperature and thermal stress constraints 

m 

This composite shell will be bonded to the outer metallic mirror substrate; polished 
aluminum, copper or nickel are likely candidates. The parabolic mirror could be 
mass-produced at a reasonable price using an electro- forming process (see Fig 2.11). 
A highly-polished metal mandrel would first be cut using a high precision diamond- 
turning lathe; the female mold would then be constructed by electro-depositing metal 
upon this mandrel. Once reinforced, this female mold could generate thousands of 
identical, light weight parabolic mirrors. A second electro-forming and machining 
process could possibly be used to add the coolant channels on the mirror inner surface. 
As indicated above the afterbody mir ror substrate will require active cooling for 
two reasons: 1) in its role as a plug nozzle there will be convective and radiative 
heat transfer from the laser-heated exhaust gases; 2) also, as the primary receptive 
optic, even 0.1 percent absorption of the incident laser energy will result in enormous 
absorbed heat loads that will ruin the mir ror figure quality. The obvious solution 
is to provide active cooling, using onboard cryogenic fluids to flow through channels 
formed into the metallic mirror substrate structure. 

Careful attention must be applied to the design of these afterbody cooling chan- 
nels. The coolant inlet manifold is positioned at the aft end of the mirror. After 
flowing through mirror cooling channels, the LN 2 enters the oulet manifold at the top 
of the mirror, which delivers the propellant to the shroud support strut “T-tails” . The 
radius ratio between mirror outlet to inlet is 8:1. This could result in an unacceptable 
pressure drop if the coolant channels are not properly engineered. Coolant channel 
passages could vary in depth from bottom to top in order to maintain constant cross- 
sectional flow area; also individual ch ann els may split in two at each doubling of the 
afterbody radius. Thinner walls and more, smaller, channels will enhance cooling and 
chamber life [11]. 

The outer surface of the afterbody mir ror structure must survive a harsh oper- 
ating environment without degradation of reflectivity or loss of figure contour [17]. 
Mirror surface defects produced during powered flight could create localized regions 
of seriously lowered reflectivity. The resultant increase in absorbed laser flux could 
create localized hot spots and a mirror ‘bum-through’ [18]. Mission failure would be 
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Figure 2.11: Proposed technique for “outside-in” fabrication of regeneratively cooled 
primary optic (not to scale). 
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inevitable. 

A diamond film outer coating on the receptive optic surface could provide the 
required exterior hardness, and enable the mirror to withstand abrasion from parti- 
cles impacting at extreme velocities. Recent research has produced new diamond film 
manufacturing techniques, now available at a reasonable cost [14]. It uses chemical va- 
por deposition (CVD) techniques, a process normally applied in making semiconduc- 
tor films. This technique has the potential to make continuous films with controlled 
material properties. Diamond film coatings are strong, stiff, hard, and slippery; they 
can transmit light from the far infrared through ultraviolet, and are five times more 
thermally conductive than copper. Diamond films can also stand up to high power 
cycling without cracking. 


2.7 Primary Internal Support Structure 

The LTD primary internal support structure (i.e., the backbone) interconnects the 
three primary sub- structures: the forebody, shroud and afterbody (see Fig. 2.12). The 
“backbone” is an annular internal support ring with a U-channel cross-section, viewed 
sideways. Conceived as a molded construction graphite/epoxy structure, this main 
support ring must perform a multitude of functions. Internal trusses that support 
the electronics payload and propellant storage tanks will attach to the ring. Also, the 
shroud support struts will be bolted into the ring, and transmit thrust loads from the 
shroud pulsejet engine to the rest of the vehicle. The main support ring may therefore 
include active or passive damping mechanisms to attenuate engine vibrations. The 
upper surface of the ring provides receptacles for attaching the forebody shell by 
means of hinges and a “Super- Zip” separation system. Recesses must also be provided 
for storage of the “ring-retina” sensors during powered flight. The bottom portion 
of the primary structure will enclose the annular inlet shutters. The structural ring 
must exhibit near-zero deflection in flight, to prevent binding of the shutters when 
they are deployed and retracted. The afterbody will bolt onto the lower main support 
ring surface. 

Detailed mechanical design of the LTD internal support structure will require finite 
element modeling to accurately evaluate structural performance. Similar carbon fiber 
structural components have been successfully demonstrated in the optical bench for 
the High Resolution Spectrograph Experiment used on the Hubble Space Telescope, 
and wing spars for advanced fighter aircraft [15]. 
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2.8 LTD Mass Requirements 

The LTD has a dry mass of 120 kg, and is filled (just prior to launch) with 140 kg of 
liquid nitrogen and 2 kg of supercritical helium gas. One tenth of the liquid nitrogen 
(14 kg) is allocated for use as sacrificial coolant during the airbreathing engine cycle 
to cool the laser mirror and engine hot section. Of the remaining nitrogen, several 
kilograms will be trapped in the propellant feed system at engine burnout. The 
result is a mass fraction of two for the vehicle which appears to be a reasonable goal. 
Table 2.1 lists the vehicle mass breakdown by component. 
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Table 2.1: LTD Vehicle Mass Breakdown. 


Section 

Mass (kg) 

FOREBODY 


C-C Nosecap 

1.0 

Primary Structure 

4.0 

TPS(FRCI) 

3.5 

Internal Support Frame 

10.0 

AFTERBODY 


Parabolic Mirror (Faceplate) 

12.0 

Mirror Support Structure 

8.0 

Manifolds 

2.0 

THRUST STRUCTURE 


Shroud 

15.0 

Strut Support Structure (24) 

15.0 

Air Inlet Shutters 

5.0 

PROPELLANT MANAGEMENT SYSTEM 


LN2 Tank 

8.0 

He Tank 

2.5 

Insulation 

1.0 

Feed Lines/Regulators 

5.5 

SATELLITE SYSTEMS 


Mechanical Actuators 

3.0 

Electrical 

5.0 

Electronics 

5.0 

Attitude Control 

2.5 

GROWTH CONTINGENCY (10%) 

12.0 


VEHICLE TOTAL DRY MASS 120.0 














38 CHAPTER 2. VEHICLE CONFIGURATION AND STRUCTURAL DESIGN 


REFERENCES 

[1] Kelly, H.N., Rummer, D.R. and Jackson, L.R., “Research in Structures and Ma- 

terials for Future Space Transportation — An Overview,” Journal of Spacecraft 
and Rockets, Volume 20, January- February 1983, pp. 85-96. 

[2] Jackson, L. Robert, et. al., “Hypersonic Structures and Materials — A Progress 

Report,” Aerospace America, Volume 25, October 1987, pp. 24-30. 

[3] Venneri, S.L., et. al., “Future Directions in Materials and Structures for Space 

Applications,” ALAA Systems Technology Conference, 1986, AIAA Paper 86- 
1185. 

[4] Eaton, D.C.G. and Slachmuylder, E.J.,“The Use of Advanced Materials in Space 

Structure Applications,” IAF International Astronautical Conference, 1987, 
LAF Paper 87-305. 

[5] Pitts, W.C. and Murbach, M.S., “Heatshield Design for Transatmospheric Vehi- 

cles,” AIAA/ASME Joint Thermophysics and Heat Transfer Conference, June 
1986, AIAA Paper 86-1258. 

[6] “SuperZip Separation Systems,” Lockheed Missiles and Space Company, Sunny- 
vale, Ca, 1986. 

[7] Johnston, Patrick J. et. al., “Fitting Aerodynamics and Propulsion into the Puz- 

zle,” Aerospace America, Volume 25, September 1987, pp. 32-37. 

[8] Waltrup, P.J., Billig, F.S., and Stockbridge, R.D., “A Procedure for Optimizing 

the Design of Scramjet Engines,” Journal of Spacecraft and Rockets, Volume 
16, Number 3, May-June 1979, pp. 163-171. 

[9] Truitt, Robert Wesley, Hypersonic Aerodynamics, Ronald Press, N.Y., 1959. 

[10] Miele, Angelo, “On the Theory of Optimum Aerodynamic Shapes,” Thirteenth 
AFOSR Science Seminar, Vistas in Science, University of New Mexico, 1968. 

[11] Schoenman, Leonard, “High Pressure Propulsion Advanced Concepts for Cool- 
ing,” IAF 36th International Congress, 1985, IAF Paper 85-191. 

[ 12 ] LaBotz, R.J. and Vallor,H., “Performance of a Transpiration- Regenerative Cooled 
Rocket Thrust Chamber,” NASA CR-159742, September 1979. 

[ 13 ] “The Aircraft Gas Turbine Engine and its Operation,” United Technologies 
Corporation, 1988. 



2.8. LTD MASS REQUIREMENTS 


39 


[14] Brown, Alan S., “Diamonds Shine Brightly in Aerospace Future,” Volume 25, 
November 1982, pp. 12-15, 37. 

[15] “Carbon Fiber Composite Materials,” B[ercules Aerospace Company, 1988. 

[16] Ghosk, A.K., and Hamilton, C.H., “Superplastic Forming and Diffusion Bonding 
of Titanium Alloys,” Rockwell International Science Center, Thousand Pairs 
USA, January, 1986. 

[17] Porteus, J.O. and Seitel, Steven C., “Absolute Onset of Optical Surface Dam- 
age using Distributed Defect Ensembles,” U.S. Naval Weapons Center, Physics 
Division, China Lake, CA , March , 1984. 

[18] Draggoo, Vaughn G., et al., “ Optical Coating Measurement for High Power 
Laser Systems,” Lawrence Livermore National Laboratory, SPIE O-E/LASE 
’86 Technical Symposium, Los Angeles, CA, January 1986. 

[19] Camarda, Charles J. and Ransone, Philip O., “Carbon- Carbon Panels Cooled 
by Heat Pipes,” NASA Langley Research Center, NASA Tech Briefs, July, 1989. 



40 CHAPTER 2. VEHICLE CONFIGURATION AND STRUCTURAL DESIGN 



Chapter 3 

PROPELLANT MANAGEMENT 


3.1 Introduction 

The LTD propulsion system requires onboard reaction propellant after transitioning 
to the rocket mode at Mach 5 and 100 kft altitude, when the atmosphere is too 
thin for airbreathing engines to produce thrust. The annular inlet is closed, and 
the mass flow rate of propellant is increased. After being expelled by an injector 
array into the annular laser beam focus, the propellant is super heated (by the beam) 
into an incandescent gas. With the inlet closed off, the LTD afterbody and shroud 
act together to form a rocket engine absorption chamber and nozzle. The heated, 
high pressure exhaust gases expand against the LTD afterbody which functions as a 
high expansion ratio nozzle for efficient thrust production. The rocket engine mode 
continues until all usable onboard propellant has been depleted; the final portion (i.e., 
ullage gas) evaporates and is used by the cold gas jets for satellite attitude control 
maneuvers. 


3.2 Propellant Management Criteria 

A suitable propellant feed system is needed to transfer propellant from on-boaxd 
storage tanks to the laser- heated rocket thrust chamber, at specific fluid flow rates 
and pressures. The feed system design depends upon the mission parameters, vehicle 
size, and mass ratio. Due to the small propellant tank volume, a pressurized tank 
system (using stored gas for pressurant) has been chosen for its simplicity, reliability 
and minimum weight penalty [1], A schematic of the propellant management system 
is shown in Fig. 3.1. 
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The following criteria are essential to the design of the LTD pressure-fed propellant 
system [1]. 

1. Maximum operating pressures and “burst” pressures (i.e., factor- of- safety) of 
the propellant and pressurant tanks. 

2. Propellant mass and volume as a function of operating temperature and pres- 
sure. 

3. Tank mass and internal volumes for liquid propellant and gas pressurant. 

4. Ullage as a fraction of the toted tank volume. 

5. Remaining propellant and pressurant mass at engine cutoff available for satellite 
attitude control. 

6. Design of auxiliary valves, pressure regulators, etc. for the complete plumbing 
system. 


3.3 System Description and Layout 

The LTD propellant management system consists of a 28 in. diameter pressure ves- 
sel for liquid nitrogen (LN 2 ) propellant, and a 12 in. diameter vessel (with internal 
heating system) for the gaseous helium (He) pressurant. A start /shutoff valve, over- 
pressure vent valve, and pressure regulator are required for each tank, along with 
temperature, pressure and flow rate sensors. During normal operation, the LTD has 
a non- vented fluid storage system; however, safety vent valves are essential to prevent 
catastrophic failure from over-pressurization (see Fig. 3.2). 

The computer-controlled propellant management system must: 1) deliver the 

proper mass flow rate of liquid nitrogen coolant to critical actively- cooled engine 
components during the airbreathing propulsive mode; and 2) provide uniform propel- 
lant mass flow during the rocket mode. Sensors distributed over the vehicled outer 
skin and engine will monitor surface temperatures. The mass flow rate of nitrogen 
coolant will be adjusted to compensate for variations in laser power absorbed by the 
primary mirror during the airbreathing mode. Laser power input is constant during 
the rocket mode. 

Propellant will flow from the bottom to top of the mirror (while providing re- 
generative cooling), and then enter the nickel “T-tail” strut structure. The auxiliary 
flow passages (i.e., for the airbreathing mode) to the strut shroud cooling system 
will be separate from the main rocket propellant feed lines. Coolant flow rates must 
be monitored carefully to avoid excessive and wasteful use of LNj. Propellant flow 
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Figure 3.2: A Supercritical Cryogenic Gas Storage System (Ref. 4). 

passages terminate at the shroud engine injector array. The injector orifices are de- 
signed to act as “acoustic valves” , automatically controlling the injection sequence of 
propellant into the rocket engine thrust chamber. Precise control of injector orifice 
pressure is necessary for optimum propellant utilization during flight. 

3.4 Cryogenic Fluid Options 

The LTD combined-cycle, laser propulsion system employs onboard cryogenic fluids 
for reaction propellant, sacrificial coolant and feed system pressurant. Liquid nitrogen 
has been chosen as the propellant for several reasons. Although liquid hydrogen (LH 2 ) 
would be the most desirable propellant because of its low molecular weight (and higher 
rocket specific impulse), only 14 kg of LH 2 would fit in the LTD’s 28 in. diameter 
tank; in contrast, 140 kg of LN 2 would occupy the same volume. Heated to high 
temperature and pressure, nitrogen can produce specific impulses in the range of 725 
to 1025 seconds — depending whether the exhaust gases are dissociated or not. (LN 2 
is also favored because it is an inert and exceptionally clean coolant, a consideration 
which is especially important for high power laser optics [2]). 

High pressure, super-critical helium was chosen for the pressurant because of its 
superiority as an inert agent with a very low boiling point. The helium tank will 
contain a heat exchanger system (see Fig. 3.3). The LTD pressurant system is modeled 
after the Apollo Lunar module descent rockets. This pressurant system operated at 
near liquid helium temperatures (5.2 K), thereby enabling the maximum pressurant 
density [3j. 
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PHYSICAL PROPERTY 

N2 

He 

MOLAR MASS 

28.0134 g/mol 

4.0026 g/mol 

CRITICAL TEMPERATURE 

126.20 K 

5.2 K 

CRITICAL PRESSURE 

3.400 MPa (493.1 psi) 

0.227 MPa (32.99 psi) 

CRITICAL MOLAR DENSITY 

0.01 121 mol/cm 3 

0.0174 mol/cm 3 


Table 3.1: Physical properties of Nitrogen and Helium 


3.5 Cryogenic Storage and Handling Considera- 
tions 

Cryogenic fluids for the LTD will necessitate launch site storage and handling of 
extremely cold gases in their supercritical state. Table 3.1 lists several physical prop- 
erties of nitrogen and helium. Such cryogenic fluids are known to be compressible, a 
desirable characteristic for achieving maximum fluid mass within fixed storage tank 
volume constraints. Strictly speaking, cryogenic fluids are normally stored ( at the 
launch site) as compressed liquids, rather than in a supercritical state. However, 
cryogenic fluids behave as a high density gas in either regime. “Supercritical” stor- 
age simply means that the fluid pressure is elevated beyond the “critical” pressure, 
and that it remains a single phase for the duration of the mission [4]. A pressure- 
enthalpy diagram (from ref. 4) for supercritical storage of cryogenic fluids is displayed 
in Figure 3.4. 

The charging process for a typical supercritical system normally begins with a 
purge and vacuum pumping of the tank and all its lines. In Figure 3.4 state “1” 
represents the “fill” condition, after which heat is added to increase the pressure 
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Figure 3.3: Helium pressurization system using heaters in storage vessel. 



Figure 3.4: Pressure-enthalpy diagram of supercritical storage (Ref. 4). 
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while maintaining a constant density (at a constant volume). In moving from state 
“1” to “2,” the liquid expands and the ullage gas condenses until the entire vol um e 
is filled with liquid. From state “2” to “3” thermal energy is added until the fluid 
rises above critical pressure. For a nonvented system, the heat that is required to 
pressurize the vessel is given by: 

Q = M f Hi - MfPfV f + MiPiVt (3.1) 

where 


Q = 

heat energy 

{KJ} 

M = 

stored Mass 

{Kg} 

H = 

specific enthalpy 

{KJ/Kg } 

P = 

pressure 

{KP a] 

V = 

specific volume 

{">7 Kg] 

f = 

final fluid state 


i = 

initial fluid state 



During withdrawal of fluid, constant supercritical operation is maintained from 
state ‘ 3” to u 4” by additional heat input to the stored liquid. The supercritical 
state insures that the stored mass remains a single phase, and fluid expulsion is 
assured because of the high operating pressure. The thermal energy input required 
to maintain pressure in a supercritical system is given by: 

Q = -D(dH/dD) p (3.2) 

where D is density (Kg/m 3 ) and p is constant pressure. The solution of Eqn. 3.2 
yields the heat energy per Kg of fluid withdrawn, for maintaining constant pressure 
as a function of density. The total required heat input is delivered to the internal 
tank by the internal heater and from environmental, electrical, and any other heat 
sources associated with on-board LTD systems. Individual LTD launch trajectories 
will require a specific fluid flow rate schedule; the vessel heat exchanger and heat-input 
rates must be carefully designed (and controlled) to avoid excessive pressurization and 
the need for venting. 

The helium gas expansion process will be polytropic. The following correlation 
can be used to obtain the final gets temperature (in the storage tank), as a function 
of pressure: 

T,/T i = (3.3) 

with n = (C -C r )/(C -C„) 

where C is the specific heat ( which will vary during system operation) and C p and 
C v are the specific heat at constant pressure and volume, respectively. The exponent 
n can be estimated analytically and verified experimentally [5]. 
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3.6 Propellant Tank Pressurization System Require- 
ments 

Determination of the pressnrant ga s weight for a cryogenic propellant is complex and 
defies exact analytical treatment because of the interdependent transient phenomema 
of heat and mass transfer that occur simultaneously in a propellant tank. Mathe- 
matical models describing the internal thermodynamics of tank pressurization have 
been developed by various investigators [13]. Although the most accurate method of 
predicting pressurant requirements is to use a computer program that has been veri- 
fied by experiment, initial design studies require fast, reasonably accurate analytical 
methods. 

The total mass of pressurant gas required is a function of the ullage mean tem- 
perature at cutoff, derived with the gas equation of state: 

PVM 

Wtotal = (3-4) 

aRT m 

where 


Wtotal 

= total pressurant weight 

P 

= ullage pressure 

V 

= propellant tank volume 

a 

= gas compressibility factor 

R 

= universal gas constant 

M m 

= pressurant molecular weight 

T m 

= mean ullage temperature 


Therefore the total pressurant mass required may be calculated if the ullage mean 
temperature at cutoff can be deter min ed. Dimensional analysis provides a relation 
for T m as a function of twelve system design variables, seven physical properties, the 
mechanical equivalent of heat, and the gravitational constant [14]. For the design 
of the LTD pressurization system, vehicle parameters such as tank volume, engine 
flow rate, tank material, etc., determined by mission profile, are fixed input values. 
However, there are various controllable parameters in the LTD pressurization system 
that can be used to optimize the system without affecting basic vehicle characteristics 

[14]- 

3.7 Pressure Vessel Design 

The LN 2 and He pressure vessels will utilize filament- wound, prestressed composite 
fiber (and metal liner) construction (see Fig. 3.5). This construction technique con- 
sists of a pretensioned fiber that is wound over a stretch-formed meted liner that acts 
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Figure 3.5: Single wall composite fiber cryogenic storage tank. 

as a vapor barrier. Proper choice of prestress and operating stress states enables 
ultra- lightweight, high pressure storage vessels. The metal finer is prestressed with 
compressive load to increase its strain limit to the point of rupture; this permits both 
the metal finer and fiber wrap to achieve their ultimate strength capabilities simulta- 
neously at rupture [7]. Candidate fibers include Kelvar-49 for the nitrogen tank, and 
graphite for the helium tank. 

The LN 2 tank design is simple: helium pressurant flows in from the top, expelling 
LN 2 out the bottom. Several problems may arise: 

• solubility of pressurant into propellant. 

• evaporation of propellant at interface into ullage gas. 

• propellant sloshing during flight. 

One well proven solution to all these problems is the “positive-expulsion” tank, 
which contains an integrally welded stainless steel ring-reinforced diaphragm (see 
Fig. 3.6). This tank design inherently provides compatibility for a wide range of fluids, 
excellent center of gravity control (during the expulsion mode), and high resistance to 
dynamic loads [8]. Application to the LTD may require alternative designs for fluid 
flow lines. 
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PRESSURANT 



Figure 3.6: Positive expulsion tank; successive positions during diaphragm reversal 
(Ref. 8). 


3.8 Thermal Insulation Requirements 

External thermal insulation is an essential requirement for the cryogenic storage tanks. 
Cryogenic insulation minimizes propellant boiloff and prevents liquefaction of air and 
the condensation of water vapor/ice onto the outer tank wall (cryopumping) [9]. The 
LTD mission includes four phases for which effective cryogenic insulation is necessary: 

• ground hold period 

• compressed air cannon launch 

• laser boost phase 

• orbital phase 

Two candidate materials have been evaluated for external insulation: rigid closed 
cell polyvinyl chloride foam (PVC), and polymethacrylimide foam (Rohacell 110WF). 
Table 3.2 lists typical material properties. The PVC foam has a lower density, but 
higher thermal conductivity. Insulative performance of the two options is similar 
when equal weight layers are compared. Both foams are widely used on LOX and 
LN 2 storage tanks and can be expected to perform well in the LTD [10]. 
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POLYURETHANE FOAM 
INSULATION 



Figure 3.7: Schematic of single wall cryogenic gas storage tank for heat transfer 
analysis • 

In addition to foam insulation, thermal radiation barriers can also reduce heat 
leaks from the forebody and afterbody exteriors. Aluminized mylar layers applied to 
the exterior of the foam insulation will function as both vapor barrier and thermal 
radiation barrier. Figure 3.7 illustrates an insulated single wall cryogenic storage tank 
schematic for future heat transfer analyses, that may focus on, amongst other issues, 
the optimization of the propellant tank heat leak factor, Q/A. 

The tank support design also has a great influence on the tank’s heat leak factor 
because conduction heat transfer through the truss and tank straps is inevitable. 
Each storage tank is attached to a truss support structure by two belts that wrap 
completely around the tank and act as a cradle. This proven method for tank support 
evenly distributes flight dynamic loads onto the tank surface. Fluid lines are routed 
over the outer surface of the tank supports to reduce conduction heat tr an sfer 

3.9 Conclusions 

The physical processes which take place during the expulsion of a liquid propellant 
from a storage tank by a gas are numerous and difficult to analyze. Applicable data 
for the LTD system is limited. Thus the basis for an analytical approach is uncertain. 
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Figure 3.8: Analytical model of pressurized propellant tank (Ref. 15). 


As a result, the initial design calculations for the propellant management system must 
be considered approximate. Fluid properties will undergo wide variations from initial 
to final states, and cannot be assumed constant. A more detailed analytical model of 
the physical process is required. 

A logical next step in the propellant management system design is the application 
of numerical methods for the pressurization and expulsion processes (see Ref. 15). 
Figure 3.8 illustrates an analytical model of a propellant tank internal volume, and 
algorithms for pressurization and expulsion processes are shown in Figure 3.9. 

These algorithms may be implemented in two computer programs: first, to predict 
helium pressurant gas requirements for the pressurization (ramp) phase of bringing 
the propellant tank up to its operating pressure; then, second, to predict system 
characteristics during expulsion of the LN 2 propellant at a uniform pressure. This 
method of analysis will involve a numerical solution of the temperature and velocity 
functions for the tank ullage at a discrete set of points in time and space [15]. 
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Figure 3.9: Logic diagrams of pressurization and expulsion process (Ref. 15) 
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PVC FOAM INSULATION 


MASS DENSITY 

50 kg/m 3 

THERMAL EXPANSION COEFFICIENT 

36 XI O' 6 mm/mm-K 

TEMPERATURE LIMIT 

120°C 


THERMAL CONDUCTIVITY : 

W/m-K 

20 K 

0.2 XI O' 2 

72 K 

1 .2 X 1 0 2 

300 K 

3.5 X 10 ' 2 


Table 3.2: Physical Properties of External Foam Insulation 
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Chapter 4 


INJECTOR DESIGN: THE LASER HEATED 
ROCKET MODE 


4.1 Introduction 


The purpose of this chapter is to describe the basic components of four candidate 
LTD propellant injector arrays, with an emphasis on the fundamental injector design 
methodology. Essential features of the four proposed injector arrays are outlined; 
each configuration has particular advantages and disadvantages. An experimental 
testing program is needed to generate performance data and evolve the final injector 
configuration. 

The function of the injector is to feed a controlled flow of propellant to the LTD 
thrust chamber during the repetitively-pulsed rocket propulsion mode. Liquid nitro- 
gen propellant is supplied to a linear array of injector elements by a series of flow 
passages routed through the nickel shroud/strut “T-tail” structure. Streams of su- 
percritical nitrogen rapidly vaporize as they expand out of numerous small orifices in 
the shroud lower surface, and arrive as a uniform flowfleld at the laser beam focus. At 
the beginning of the laser pulse, laser-induced breakdown is triggered, whereupon the 
resulting plasma strongly absorbs the beam, creating a laser-supported detonation 
wave. The wave propagates through the N 2 gas and heats it to a very high temper- 
ature and pressure. After the laser pulse is terminated, the high pressure gas then 
expands against the shroud and afterbody mirror surface. Thrust is generated by 
converting the input laser thermal energy into exhaust kinetic energy. The parabolic 
afterbody mirror acts as a large area ratio “plug” nozzle to maximize the efficiency 
of this process, in part by giving the nitrogen gas sufficient time to cool between 
laser pulses. The nozzle is large enough that the exhaust gases can re-combine and 
again become transparent (i.e., to the next incoming laser pulse), and still exert force 
against the plug nozzle surface. 
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4.2 Underlying Technologies 

The objective of a propellant injector is to enable high rocket engine performance 
and stable operation, without adversely affecting either the injector or thrust surface 
durability. In a conventional chemical rocket engine, injector design specifications 
must encompass the entire flow system geometry, including the injector element pat- 
tern or “array,” individual orifice geometry, and flow distribution system upstream of 
the orifice (i.e., manifold) [1]. The LTD injector design attempts to exploit the large 
data base of knowledge accumulated by rocket designers; however, several underlying 
technologies specific to pulsed laser propulsion must be considered. An operational 
injector design for the LTD will require research in the following areas: 

1. Propellant Technology: For a given propellant, the physical/ chemical processes 
that occur during injection, kinetics of stream breakup and mixing, vaporization 
and heat transfer must be clearly understood. The goal for the LTD is to opti- 
mize laser energy absorption by the propellant; chemical combustion processes 
are not relevant. 

2. Pulsed Flow Control: A repetitively-pulsed engine will require high frequency 
acoustic valving technology. This implies a strong coupling between propellant 
injection flow dynamics and the impulsive thrust generation environment. In 
conventional liquid rockets, similar “coupling” effects can lead to low frequency 
combustion instability (chugging), and are strenuously avoided [2]. A successful 
rocket design must eliminate all such dangerous instabilities. 

3. Structures Concepts: The shroud structure contains the injector array and 

manifolds. Since the shroud is essentially a thin annular airfoil, its narrow 
cross section imposes certain volumetric restrictions on the inject or /manifold 
design. Advanced manufacturing techniques such as platelet technology and 
diffusion bonding [3,4] offer promising solutions to the LTD’s structural design 
constraints. 

4. Laser Optics: Injected propellant mass flux must be delivered with high preci- 
sion to the primary mirror focal region. The injector array can simultaneously 
serve as a secondary optic to provide further focusing of the incident laser radi- 
ation. Alignment of the secondary optics can, of course, influence laser energy 
absorption into the propellant. 

The LTD injector designer must combine conventional liquid chemical rocket tech- 
nology with unconventional laser propulsion requirements, to produce practical solu- 
tions that exploit the unique repetitively-pulsed engine operating environment. 
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LN2 



Figure 4.1: Side view of like-impinging doublet injector element. 

4.3 The Injector Array 

For optimum gas-phase absorption of incident laser energy, it is necessary to first 
atomize the propellant; the great increase in liquid surface area accelerates vapor- 
ization. Atomization involves the breakup of liquid jets and/or sheets into droplets. 
By directing two or more streams of propellant to impinge upon each other, droplet 
and gas bubble formation occurs (see Fig. 4.1). For the LTD, the atomized propel- 
lant vaporizes very quickly upon exiting the injector array; the LN 2 will be heated 
above critical temperature and flash into a gaseous phase as the pressure drops below 
critical. 


4.4 Injector Element Design 

To design the injector element pattern, one must select the injection element (i.e., 
doublet-impinging, triplet-impinging, etc.), and element distribution (i.e., the com- 
plete annular injector array) across the lower shroud surface. Since the LTD is a 
‘monopropellant’ rocket system, the injector elements are termed ‘like-impinging’. 
The choice of injector element type is influenced by the desired mass flux contours. 
Figure 4.2 shows typical mass flux contours for a like-impinging doublet element [1], 

These injector elements should be distributed across the lower LTD shroud surface 
so as to insure the delivery of a uniform propellant mass flux at the laser optic focal 
region, a narrow (1 cm wide) circular strip above the lower shroud surface (Fig. 4.3). 
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Figure 4.2: Mass flux contour for doublet injector element (Ref. 1). 

Therefore it makes sense to place a single-file row of elements (i.e., strip injector) 
at some position near the primary mirror focus. Element orientation can be either 
‘in-line’ edge touching, or ‘parallel-overlapping’ opposed doublet (see Fig. 4.4). The 
position of one element with respect to the others should provide for maximum inter- 
element mixing and uniform mass distribution [1]. 


4.5 Injector Orifice Design 

Options in the design of injector element orifices include the geometric shape, orifice 
dimensions and the impingement angle between opposing jets of injected propellant 
(see Fig. 4.1). The impingement angle determines the position where stream mixup 
and droplet formation occurs above the injector surface. The most common impinge- 
ment angle is 60°. At impingement angles greater than 60°, propellant backsplash 
onto the injector face increases surface heat flux, resulting in injector face burnout 

[5] - 

The diameter of the orifice directly affects mixing and atomization levels produced 
by an injector element. Smaller orifices have higher performance than larger orifices, 
because the smaller droplet size results in increased vaporization rates. Impingement 
distance is another factor influenced by orifice diameter. A performance parameter 
often used is the ratio of free stream length, L it , divided by orifice diameter, D^.. 
Values of L, t /D OT from 5 to 7 are good. Larger values tend to produce misimpinge- 
ment of streams, because instabilities within a free jet are able to cause directional 
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Figure 4.5: Unconventional orifice shapes improve mixing performance. 


wandering [1]. 

Flashing injection is another process which may find application in the LTD rocket 
engine. In this process the propellant is expanded through an injector orifice until the 
supersaturated state is reached, whereupon a portion of the fluid flashes into vapor 
[11]. This is a well-known method for improving atomization properties, i.e., reducing 
droplet sizes and spray penetration, and increasing the rate of spray spread. Several 
methods have been used to promote flashing injection, depending on the nature of 
the injected liquid. Further evaluation of this concept for use in the LTD engine is 
necessary. 

The circular orifice shape has historically been the standard for most injector 
designs because passageways were drilled. With the advent of platelet technology 
(see chapter two), rectangular orifice shapes are now possible (Fig. 4.5). Experiments 
with rectangular passages have shown a measurable increase in mixing performance 
over the circular geometry [6]. 


4.6 Injector Manifold Design 

In the LTD, LN 2 propellant flows upward through the primary mirror coolant passages 
to the annular manifold at the top. From here the LN 2 enters the T-tail support struts 
by means of a two-way valve. As mentioned earlier, the twenty four T-tail support 
struts combine to form the forward half of the shroud structure. Each T-tail strut 
has its own two-way valve, which initially is set for the airbreathing mode, then later 
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Figure 4.6: Two-way valves control propellant flow to each T-tail strut assembly. 


switches to the rocket mode position. In the airbreathing mode, it will feed LN 2 to 
the T-tail leading edges and shroud impulse surface for transpiration cooling and film 
cooling, respectively. At the start of the rocket mode, it then switches the flow to the 
rocket injector array. The twenty four two-way valves must be 100% reliable during 
operation (see Fig. 4.6). 

Within each T-tail strut, LN 2 propellant must be uniformly distributed to the 
individual injector orifices by way of a mini-manifold (see Fig. 4.7). This injector 
manifold is sized to deliver a specific propellant flow rate through a minimum cross 
sectional flow area while incurring a small pressure loss. The total pressure losses 
suffered by the flow from the propellant tank to the orifice exit must be minimiz ed; 
otherwise substantial weight penalties will result for the pressurized prop ellan t tank. 
Within each T-tail strut manifold system, a central ‘downcomer’ will feed propellant 
to a distribution ring (Fig. 4.8), which contains ports designed to reduce static pres- 
sure variations in the flow region immediately behind the injector face [1], Finally, 
prior to entering each injector, the flow is put through an ultra-fine mesh filter (also 
constructed with platelet technology). 


4.7 Film Cooling for the Shroud Surface 

The lower surface of the shroud fulfills three functions: impulsive thrust surface, in- 
jector face, and (optional) secondary optic surface. The surface will be subjected 
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Figure 4.7: Typical platelet flow metering and distribution schematic for a T-tail 
strut assembly. 
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Figure 4.8: Injector face and manifold 


design details (Ref. 1). 
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to direct ablation and erosion by the pulsed laser beam and resultant high temper- 
ature/pressure blast waves (i.e., exhaust gases). The surface will also sustain large 
pulsed mechanical loads and severe thermally- induced mechanical stresses. Temper- 
atures within the laser induced blast waves can reach 10,000 K. Reduction of adverse 
thermal effects is readily accomplished by regenerative and film cooling methods [7]. 

Auxiliary film cooling for the lower shroud surface is accomplished by injecting a 
portion of the propellant to flow along the external surface and act as a buffer region 
against the high temperature exhaust gases (Fig. 4.9). These film injection sites will 
be placed near the shroud leading edge; they are completely independent of the main 
injector array . The film coolant layer thus aborbs and carries away a portion of the 
radiant energy transmitted to the shroud walls. Injected as a supercritical liquid, the 
film rapidly flashes to a gas state. The gas film temperature will rise largely as a 
result of convective heat transfer with the shroud wall. A thermal balance will be 
attained between the radiant heat input rate and coolant temperature rise. Once this 
equilibrium condition is reached, no further heat transfer occurs, and the shroud is 
maintained within structural- temperature limits [5]. The shroud and injector array 
material temperatures will be substantially lower than those in the laser-induced blast 
wave. 

Temperatures on the lower shroud surface will increase in the downstream direc- 
tion, as a function of the flow distance from the film injection site. The forward half 
of the shroud (“T-tail” strut assembly) will be metallic and have a lower temperature 
limit; the aft half of the shroud will be carbon-carbon composite, with a very high 
limit. 

4.8 Transpiration Cooling of Leading Edges 

Again using platelet technology, transpiration cooling ports are formed into the 
shroud/strut leading edges, as shown in Figure 4.7. This sacrificial coolant coun- 
teracts aerodynamic heating at the air inlet stagnation regions during hypersonic 
flight speeds. Inddently, transpiration cooling is also planned for the scramjet fuel 
spray bars on the X-30 NASP. 


4.9 Pulsed Flow Control 

To achieve efficient LN 2 mass utilization high frequency valving of the propellant is 
required [8]. These valves must operate at frequencies of 1000 Hz or greater, with 
closing and opening times of 10 -4 sec. Specially engineered propellant orifice arrays 
are known to interact with the pulsed rocket engine chamber environment in a self- 
regulating manner. In this “acoustic-valving” process, pressure fluctuations at the 
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Figure 4.9: Proposed models for active cooling. 
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injector orifice exits serve to start and stop the propellant flow, thereby metering the 
desired propellant mass [9,10]. 

The propulsion sequence for a repetitively-pulsed laser powered thruster is com- 
prised of four steps: 

1. Cold propellant is injected into the laser beam focus. 

2. Laser pulse electrical breakdown in the propellant ignites a Laser Supported 
Detonation wave (LSD wave) which produces high temperature plasma and 
throws off a strong blast wave. 

3. Shocked exhaust gas expands (generating thrust) and shuts off propellant flow, 
because of the high over-pressures (up to 600 atm). 

« 

4. Blast wave pressure decays to a low level at the orifice exit while the exhaust 
flow expands across the large plug nozzle; the propellant flow then resumes in 
preparation for arrival of the next laser pulse. 

The sequence then repeats. Propellant consumption is controlled by the peak 
injector mass flow rate, blast wave characteristics and laser pulse repetition rate. 

An important parameter in the design of self-regulating valves is the propellant 
jet velocity, as it emerges from each orifice. This jet velocity is determined by orifice 
geometry and injector operating pressure. Sufficient propellant mass must leave the 
injector face to arrive at the laser focal region in the short time between blast wave 
decay and the arrival of the next laser pulse. This time interval during which propel- 
lant will flow is governed by the blast wave decay time, ‘setback’ distance between the 
injector face and the the primary mirror focus, and the laser pulse repetion frequency. 
The decay of blast wave pressure to sea level pressure can be as fast as 10 ~ A seconds. 

Several questions arise as to the degree of propellant flow stability and uniformity 
that will result from acoustic-valving. Unstable oscillations of propellant flow (and 
resultant thrust), could cause destructive vibrations in the LTD shroud structure. 
Another design issue is exactly what critical pressure differential will stop the flow 
of injected propellant, without causing flow reversal inside the injector. Clearly, 
analytical models must be developed to predict pulsed flow control dynamics, and 
experiments are needed to validate how each design parameter affects overall engine 
performance. 


4.10 Injector Design Alternatives for the LTD 

As mentioned earlier, the LTD rocket engine needs an injector array to supply liquid 
propellant into the thrust chamber; four different injector schemes were identified. 
The rocket engine’s initial impulse surface (i.e, the lower surface of the shroud) was 
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recognized as the proper place for a circular- strip injector array. The analytical 
propulsion system model dictated the position of the injector array at the 1/3 chord 
point, measured aft from the shroud leading edge. The first design variable to be 
revealed in this injector scheme was the height, if, of the laser beam focus above the 
shroud’s impulse surface (see Fig. 4.10). 

Design No. 1 consists of a ‘fiat-plate’ impulse surface tilted at 45° to the vehicle 
centerline, with the P.O. annular focus positioned at the 1/3 chord point and lying 
exactly on the lower shroud surface (H = 0). This initial configuration contained no 
provision for injectors or active cooling of the hot structures. The central motivation 
for this early design was to illustrate the interdependence between primary optic and 
shroud geometry. 

A Super Ignitor Array (SIA) concept successfully tested at LLNL (ref. 10) pro- 
vided insight for injector design No. 2 (see Fig. 4.11). This design integrates an 
injector array into a honeycomb-patterned, one centimeter wide annular secondary 
optic (S.O.), whose function is to provide final focusing of the incident laser energy 
to intensities high enough to initiate air breakdown and plasma formation. The face 
of this super ignitor array will experience extremely high pressures and temperatures 
during operation. There exists the possibilty that unwanted plasma formation may 
occur along the many sharp * edges of the dimpled surface. The structural cooling 
requirements may tax the limits of the LTD propellant management system. The 
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geometry of the P.0, and laser beam structure may also distribute the incident laser 
energy in an unsymmetric manner across the surface of the SIA, thereby complicating 
cooling requirements. Generally, the injector array and secondary optics are pointed 
directly into the focused laser beam. 

Design No. 3 completely separates functions performed by the secondary optics 
and the injector array (see Fig. 4.12), by physically removing the injector orifices 
to a distance l d\ away from the laser focus (i.e., site of the cylindrical blast wave 
center). Note that the 1 cm wide, flat plate S.O. is inclined at an angle of 20° to the 
shroud impulse surface and that the sharply focused laser beam reflects off the S.O. 
to focus roughly 1/2 cm above this surface. The orifice exit pressure will be elevated 
above propellant critical pressure to insure that a single phase supercritical fluid flows 
throughout the LN 2 feed system. For a given injector pressure, the value for ‘d’ may 
be obtained from Sedov’s scaling law which relates the blast wave pressure versus 
radius: 

p = (cri p "t (4-D 

A plot of equation 4.1 is given in Figure 4.13. The results are based on the 
following conditions, which represent typical operating values during the rocket mode: 


altitude 

= 30 km 

speed 

= Mach 5 

T r*f 

= 0.5 cm 

P~f 

= 500 atm 


By moving the injector array away from the region of high blast wave pressures, 
injector orifice pressures effective for pulsed flow control can be lowered to levels 
conventionally used in pressurized propellant tank systems. By varying the distance 
l d\ the optimum acoustic- valving effect can be obtained. The annular S.O. is aligned 
to an angle of 25° with the vehicle centerline; a ray trace from the .707 m b eam 
diameter (on the P.O.) is normal to this S.O. surface (see Fig. 4.14). Assuming a 
‘top-hat’ beam intensity distribution exactly aligned with the vehicle central axis, 
laser energy will be uniformly focused upon this S.O. surface. Design No. 3 still faces 
the potential pitfall of high radiation heat transfer to the S.O. surface. 

Design No. 4 removes the final focusing function of the secondary optic until after 
the beam has passed through the focus (see Fig. 4.15). The final P.O. focus lies in 
the plane of the lower shroud surface, as in design No. 1. Hence, the S.O. serves only 
to reflect that portion of the laser radiation, which passes completely through the 
line- breakdown region unabsorbed, back into the P.O. focus. Radiation heat transfer 
to the S.O. should be significantly reduced from the other designs. 
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Figure 4.11: Shroud /Injector assembly, Design No. 2. 
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Figure 4.14: Shroud/Injector assembly, Design No. 4. 


4.11 Conclusions 

An examination of injector design fundamentals as they relate to laser propulsion has 
enabled identification of critical technology issues facing the LTD shroud /injector 
assembly. The next step is to develop a detailed analytical injector model that 
integrates propellant system design parameters with the vehicle configuration and 
laser-propulsive physics. 

An important issue that must be examined in this modeling effort is the integrated 
LN 2 propellant mass per laser pulse as it relates to engine PRF, thrust duration, 
and toted propellant consumed in the rocket mode; the latter must match the laser- 
heated rocket specific impulse projections. A detailed analysis of the LN 2 mass flow 
rate versus blast wave pressure over time is needed. Another critical issue is the 
injected stream impingement site: should it be coincident with the laser beam focus, 
or perhaps at some shorter distance, to insure that a wide-angle, vaporized spray 
arrives at the focus. 

Satisfactory performance for the shroud/injector assembly requires a high degree 
of reliability, propellant flow stability, and compatibility with the pulsed laser-heating 
dynamics. Although analytical performance projections will aid in the development 
process, a rigorous experimented verification program is mandatory. Finally, with 
the advent of advanced fabrication techniques (i.e., platelet structures and diffusion 
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bonding), and high temperature materials, unconventional injector designs for laser 
propulsion engines might be manufactured with relative ease. 
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Chapter 5 

FINITE ELEMENT MODELING AND 
ANALYSIS 


5.1 Introduction 

Finite Element Analysis techniques were used to evaluate the annular shroud and 
support strut design proposed for the Lightcraft Technology Demonstrator. Both 
static and dynamic analysis were implemented to extract the stress and modal be- 
havior. Pulsed airbreathing or rocket engines must avoid the natural frequencies of 
the vehicle structure, or catastrophic failure of the system will result. 

During previous research efforts, analysis of a similar ann ular shroud (i.e., for 
the Apollo Lightcraft) has been accomplished using a highly simplified analytical ap- 
proach to determining shroud natural frequencies. Although beneficial in determining 
‘ballpark’ concepts of shroud dynamics, the previous analysis simply does not provide 
sufficient useful information. 

However, by utilizing available computer based 3-D Finite Element Analysis tools, 
a much more accurate natural frequency analysis was obtained. In addition to its 
much greater accuracy, the Finite Element Modeling output facility allows for the 
animated viewing of mode shapes, as well as providing numerical values for solutions. 
This allows the designer to evaluate which components are dominating the model 
characteristics, and to redesign accordingly. 


5.2 Finite Element Mo deling/ Analysis Package 

The Finite Element Modeling (FEM) and Finite Element Analysis (FEA) package 
is a design and analysis tool that has recently become practical for use by the en- 
gineer/designer. It is a computer software package that allows static and dyn ami c 
(i.e. stress and vibrational) analysis of structures, as affected by operating conditions, 
including physical and thermal loads, boundary constraints, geometry, and material 
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Figure 5.1: Finite element model of the Shroud and its support structure. 


properties. Specific structures can be analyzed for stresses, deflections, and natural 
frequencies of vibration. This package is a powerful tool for analysis of existing de- 
signs, as well as testing of conceptual designs: it can reveal problem areas without 
necessitating creation of physical models or prototypes. 

The FEM package implemented for the LTD shroud mechanical analysis was the 
IBM Computer Assisted Engineering Design System (CAEDS™), presently ru nnin g 
on Rensselaer Polytechnic Instituted IBM 4341 computing facility. 

It was determined that a full three-dimensional Finite Element Model of the 
shroud and its 24 support struts would be the most accurate and efficient finite 
element model for this analysis (see Fig. 5.1). The shroud was modeled by solid lin- 
ear brick elements; thin shell elements of varying thickness were used to model the 
tapered struts. 

Due to computer limitations, it was necessary to use an interactive analysis pro- 
cess, starting with the global modes for the shroud and strut system. Then, to extract 
more detailed information about the modal behavior of the struts themselves, a single 
section (l/24th of the shroud structure) was examined. This allowed a much more 
dense finite element mesh with higher-order parabolics to be applied to the strut 
structure, thus yielding a more accurate and complete analysis of strut behavior. 

Two different analyses were performed on the “single” sections. One included 
the entire section of the shroud and its associated strut (see Fig. 5.2). This second 
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Figure 5.2: T-tail Shroud/Strut assembly including carbon-carbon lower half (“entire” 
section). 


analysis removed the trailing edge of the shroud, which would be made of carbon- 
carbon composites (see Fig. 5.3). The reasons for considering the second case is 
that the strut and leading edge of the shroud structure will be manufactured as 24 
individual sections; however, the carbon-carbon trailing edge will be a continuous 
ring. Thus, the leading edge of the shroud and its associated strut sections can be 
considered a single structure. It is important to note that the smallest single structure 
will exhibit the highest modal frequencies; hence, it could become the critical case at 
comparatively high engine frequencies. 

The finite element models of the struts and shroud leading edge were given the 
properties of Nickel-200; the trailing edge was assigned the properties of carbon- 
carbon composite material. These material properties include density, Elastic mod- 
ulus, and Poisson’s ratio. The struts are considered fixed at their internal vehicle 
mounts (i.e., the inside edge for the global analysis; see Fig. 5.4), and at the mount- 
ing holes for the strut section analysis. 


5.3 Advantages/Disadvantages of the 3-D Model 

The three dimensional FEM has many advantages, and several disadvantages. The 
most significant disadvantages are: 
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Figure 5.3: T-tail Shroud/Strut assembly less carbon-carbon ring lower half (“cutoff” 
section). 



Figure 5.4: Strut mounting detail. 
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1. Time consuming - the 3-D modeling process takes an extreme amount of time 
to implement. For example, the finite element analysis of the Space Shuttle’s 
main fuel pump took twelve man-years to complete. 

2. Computer Costs - much numerical computation is needed to reduce the three 
dimensional stress matricies. CPU time required for the shroud analysis was 
37 minutes. This translates to a cost around $70.00; it is not uncommon to see 
several thousand dollar runs for complex 3-D models. 

while the most significant advantages are: 

1. Very accurate — the 3-D analysis allows for up to six degrees of freedom for each 
structure node; 2-D models allow for only three degrees. Therefore the actual 
behavior of the structure at the node points is much more closely represented. 

2. Excellent Analysis Output Capabilities — the 3-D CAEDS™ System model 
has high powered output capabilities, including graphic and numerical displays 
of deflections, stress concentration, and mode shapes. This permits the engineer 
^or designer to more fully understand the true behavior of the structure. 

3. Modifiable — once the 3-D model has been created, it can be modified to 
study the effects of design changes. It is possible to change material properties, 
structure geometry, and operating conditions such as external loading and ac- 
celerations. This gives the designer great freedom to test different theories or 
ideas without having to build a new model each time. 

From the above it should be clear that the advantages of the 3-D model far 
outweigh the disadvantages; it was the correct choice for the LTD shroud stress and 
and vibrational analysis. 

5.4 Results 

5.4.1 Volumes and Masses 

The detailed section modeling also enabled a very accurate measure of volumes and 
masses for the various LTD structural components. The initial dynamic analysis 
assumed completely solid elements and thus had component masses which were con- 
siderably higher than the design limits (see Table 5.1). However, the static analysis 
of the strut-shroud section (discussed subsequently) revealed very low stresses in the 
shroud structure, as compared with the strut. Consequently, the shroud will utilize 
a semi-hollow structure which will permit the design mass limits to be reached. 
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COMPONENT 

VOLUME (cm 3 ) 

SOLID 
MASS (kg) 

DESIGN 
LIMIT (kg) 

STRUTS 

1,964 

15.8 

16.0 

LEADING EDGE 

6,074 

48.9 

10.0 

TRAILING EDGE 

7,877 

14.3 

7.0 


Table 5.1: Volumes and Masses for Shroud Structure. 


5.4.2 Global Analysis 

The global analysis involved the implementation of a 3-D linear solid finite element 
model for the entire shroud-strut structure; from this, the first 10 natural modes of 
global vibration were extracted. Results included both numeric data for the modal 
frequencies and graphic depictions of the various mode shapes. 

The numerical results are presented in Table 5.4. Note that the frequencies for 
modes 2, 3, 4 and 5, as well as 6 and 7, (respectively) are identical. This phenomenon 
is principally due to the orthogonal, or coupled modes. Note in Table 5.4 that the 
structural component (shroud or support) is the primary oscillator for a given mode. 
Figure 5.5 presents plots of the various mode shapes (deformations are exaggerated). 
Only the lower numbered mode is included for each coupled mode, since the orthog- 
onal mode is identical in shape to its counterpart. 


5.4.3 Local Modes 

The first fifteen natural vibration modes were extracted for both the “entire” section 
and “cutoff” section finite element models. The natural frequencies of the “cutoff” 
section structure are considered the critical frequencies to be compared against engine 
pulse repetiton frequency (PRF); the“entire” section model was considered primarily 



Figure 5.5: Various mode shapes for the entire LTD shroud (deformations are exag- 
gerated). 


c 
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MODE No. 

COMPONENT 

FREQUENCY (Hz) 

1 

STRUT 

3.0 

2 - 3 

STRUT 

14.8 

4-5 

SHROUD 

15.1 

6-7 

SHROUD 

20.7 

8 

STRUT 

34.0 

9 - 10 

SHROUD 

34.5 


Table 5.2: Modal Data for the Shroud Global Analysis. 


to obtain its qualitative behavior. 

The modal frequency change from the entire section model to the cutoff version 
was found to be insignificant; i.e., by less than a factor of 1.5. The resulting frequencies 
of the first 15 modes of vibration are presented in Table 5.3. 

Again, the component column refers to the structural element which is found to 
demonstrate the most pronounced oscillatory behavior for the given mode shape. The 
“solid” model and “hollow” model results in Table 5.5 refer to the solidity options 
for the shroud leading and trading edges. To approximate the hollow shroud option, 
the densities of the shroud elements were just reduced to bring their masses within 
design limits (i.e., the design goal in Table 5.1). 


5.4.4 Overview of Results 

The first measure of accuracy for an analysis is whether the output makes intuitive 
sense. For the shroud vibrational analysis, both the numerical and graphical data ap- 
pear reasonable, and relatively predictable. When evaluating finite element analysis, 
the primary concern is the graphical output. A major source of error in finite element 
analysis is improper geometry definition, which is easdy identified by observing the 
graphic output. No such improper geometry was found to exist in the shroud FEA; 
thus the data is considered accurate. 
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MODE No 

COMPONENT 

SOLID MODEL HOLLOW MODEL 

FREQUENCY (Hz) FREQUENCY (HZ) 

1 

ST 

2.5 

5.1 

2 

ST 

4.2 

8.1 

3 

ST 

21.6 

40.8 

4 

ST 

37.9 

59.6 

5 

ST 

103.3 

121.4 

6 

SH 

104.5 

186.2 

7 

SH 

168.5 

197.2 

8 

ST 

188.3 

309.3 

9 

SH 

229.8 

312.8 

10 

ST 

296.4 

414.9 

11 

SH 

341.1 

454.8 

12 

ST 

393.7 

527.9 

13 

SH 

475.3 

580.2 

14 

ST 

522.5 

586.1 

15 

ST 

540.4 

703.3 

ST = 

STRUT SH = SHROUD 



Table 5.3: Modal Data for LTD Single Strut Section 
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5.5 Frequency Information and LTD Operation 

The major concern which necessitated the shroud modal analysis was the interaction 
of the LTD’s pulsejet engine “source frequency” with the response of the shroud 
and inner support structure. The shroud and support struts comprise the main 
engine/ vehicle structure; thus structural resonant frequencies must be avoided. If the 
input frequency of engine pulses correspond with destructive resonant modes of the 
structure, catastrophic failure would be probable. 

The cause of failure in structures undergoing resonant vibrations is due to the 
inherent stress associated with unusually large displacements of the oscillating struc- 
ture. Induced stress has a linear relationship with displacement. As the mode number 
increases, so do the frequencies needed to excite them. However, the maximum dis- 
placement of a given particle of mass is inversely proportioned to the input frequencies, 
as follows: 

For a mass-spring system: 

Excitation Force F = F max cos(u>i) 

Displacement x = (F max /G)sm(ujt ) 

Max. Displ. X = {F max /G) 

UTn pro 

G = yj(m u/) J + (M 2 

m = mass 

b = Damping Coefficient 
Assuming unit coefficients (i.e., m 2 + ^ Fmax !)• 

X = l/w X — Normalized Max. Displacement 
= 1/2/ w = Ang. Frequency of Excitation 
/ = Frequency of Excitation (Hz) 

Figure 5.6 clearly indicates the inverse relationship between displacement of a 
particle in an oscillating body and input frequency. Although the derivation of this 
relation is for a mass-spring system, the basic behavior may be generalized to all 
structures. This generally holds true because any structure can be approximated by 
an infinite number of mass-spring systems. 

Although large displacement responses may be obtained from lower harmonics at 
low frequency excitation, the response quickly diminishes with increasing frequency. 
Even at 50 Hz, the normalized response has fallen to only 2% that of a 1 Hz excitation 
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EXCITATION FREOUENCY (Hz) 


Figure 5.6: Normalized Displacement versus Frequency. 

as shown in Figure 5.6. Clearly, at an input frequency of 10,000 Hz, the norm aliz ed 
response will be nominal for the high mode numbers obtained. Therefore, it may 
be concluded that the LTD shroud and support strut structure is able to withstand 
stresses of induced oscillations at the current engine design PRFs of 1 to 10 kHz. Even 
though resonant frequencies may be produced in flight, extremely small amplitudes 
(of displacement) will not cause damaging stress levels. 

For next year’s research effort, it has been proposed to lightweight the entire 
craft, including the shroud and strut assembly in order for the LTD to fly on 10 x 
less laser power. With such a lightened lightcraft engine, PRF’s would be lowered 
significantly, perhaps to the range of 100 Hz to 1 kHz. As observed in the lightened 
shroud analysis (Table 5.5), modal frequencies increase with decreasing shroud mass. 
This could put destructive modes in the range of engine PRF, an effect that will need 
further attention and analysis. 

The destructive modes will be of consequence, however, mostly at liftoff, or initial 
engine startup, when the PRFs are low. The “ramp” function used on engine PRF to 
suddenly apply power will produce a wide spectrum of source excitation frequencies; 
a large amount of power will fall in the destructive frequency range (e.g., 3-50 Hz). 
In the future the structural response to different “ramp” functions should be carefully 
examined using the basic structural mode shapes and resonant frequencies. 

In the “sectional” analysis, significant strut-dominated modes were found to exist. 
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The struts were found to oscillate in and out of the strut plane, with virtually no mo- 
tion of the shroud section. These modes began to appear with mode 8 (see Table 5.3), 
at a frequency of 309.3 Hz for the lightened version. Higher frequency harmonics of the 
original “flap” mode were interspersed with shroud dominated modes (see Fig. 5.7). 

These so-called “flap” modes could prove to be dangerous. It is quite likely that 
pulsejet frequencies capable of stimulating flap modes will be present during vehicle 
operation. Another aggrevating effect would be due to air flow fluctuations around 
the strut. Hence, rigorous wind tunnel testing must be included as a part of the strut 
design validation process. 


5.6 Mode Shapes 

The graphic representation of mode shapes permits the designer to ascertain exactly 
what is happening at each of the mode frequencies. In the first analysis of the shroud 
vibrational characteristics, six unique mode shapes (and four orthogonal mode shapes) 
were obtained. 

The greatest problem with earlier analytical solutions for the shroud vibrational 
dynamics was in discerning the relative behavior of the support struts and the annular 
shroud. With this “first order” analysis, two situations could be assumed: (1) Support 
struts have zero stiffness in relation to the shroud, or (2) support struts have infinite 
stiffness in relation to the shroud. These two assumptions made it possible to find 
global modes and local modes, independent of each other. 

Predictably, the actual behavior was somewhere in between the two extremes, as 
evidenced by the intermingling of mode types. In addition, natural modes of the sup- 
port struts themselves were readily evident. Two basic types of modes were evident. 
They were defined by the structural component which acts as the primary oscillator, 
either (1) the support struts, or (2) the annular shroud. Presented in Table 5.4 is a 
listing of the 6 unique modes, their structural component, and description of their 
behavior. After examining the results of the finite element analysis, one may note 
the absence of “in-plane” elliptical modes. Such modes are not present due to the 
inherent symmetry of the shroud and support struts, about the radial plane. Thus 
the first global modes of the annular shroud are “out-of-plane” bending modes. This 
again is another diagnostic feature that would not have been revealed by an analytical 
approach. 


5.7 Conclusions of Dynamic Analysis 

The main purpose of the dynamic analysis was to ascertain the frequencies of the 
natural vibration modes for the annular shroud/ strut support structure. This is 
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Figure 5.7: Various mode shapes for the 
exaggerated). 


T-tail Strut Assembly (deformations are 
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MODE No. 

COMPONENT 

PLANE 

DESCRIPTION 

1 

ST 

IP 

TWIST MODE 

2, 3* 

ST 

OOP 

ROCK MODE 

4, 5* 

SH 

OOP 

1st BENDING MODE 

6, 7* 

SH 

OOP 

2nd BENDING MODE 

8 

ST 

OOP 

HEAVE MODE 

9. 10' 

SH 

OOP 

4th BENDING MODE 

* = ORTHOGONAL MODE 

ST = STRUT 

SH = SHROUD 

IP = IN-PLANE 

OOP = OUT-OF-PLANE 


Table 5.4: Mode Shapes and Description 


critical because of the pulsed nature of the LTD engine, for which the shroud/support 
strut is the major structural component. 

It was determined that the first six unique global modes ranged from 3.0 to 34.5 Hz. 
Maximum structure displacement (and therefore maximum induced stress) in an os- 
cillating system is inversely proportional to the excitation frequency. Therefore, the 
induced stress at kilohertz “source” frequencies projected for the LTD pulsejet engine 
will not excite destructive vibration modes in the engine/ vehicle structure. 

Although the current engine PRF will not present a design problem for the shroud- 
support strut structure, future lightened versions of the craft (with their higher modal 
frequencies) could potentially be damaged by destructive modes of vibration excited 
by engine PRF’s in the 100 to 1000 Hz range. The lower-frequency spectral content 
of engine excitation frequencies must be ex ami ned to insure that strut dominated 
“flap” modes do not occur. Such modes exhibit extremely high lateral displacements 
which are highly stressing to the strut member. Fatigue effects must also be taken 
into consideration. 

The shroud and support struts appear to have relatively equivalent stiffnesses; 
thus, their frequencies of natural vibration are interspersed. Since various mode 
numbers for the struts and shroud are of the same magnitude, the response to a given 
input frequency are qualitatively equivalent. For example, destructive modes will be 
excited in both the struts and shroud at excitation frequencies between 3 and 50 Hz; 
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however, neither will undergo destructive modes between 10,000 and 20,000 Hz. This 
is found to be a desirable characteristic of the structure. 

A basic and sound understanding of the strut dynamic behavior has now been 
obtained. Several key design points have been clarified, and the groundwork has been 
laid for future structural analysis. Static (stress) analysis may now be accomplished 
using the existing model. 


5.8 Static Analysis 

As mentioned earlier the shroud and support structure of the LTD acts as the pulsejet 
engine of the craft. The high pressure, laser- heated gases push on the underside of 
the shroud to produce thrust. This repetitively-pulsed thrust is transmitted to the 
craft via the twenty-four radial support struts. 

Both the shroud and the struts must be able to withstand the forces induced during 
maximum engine thrust and vehicle acceleration. Under this load, the integrity of 
the structure must be maintained, and the struts must not permit the shroud lower 
surface to be deflected away from the focal region of the^ncident laser beam. 


5.9 Static Analysis Model 

The purpose of the static finite element analysis was to evaluate the strut design 
(see Fig. 5.8) to uncover stress concentration areas. Since the load case is axisym- 
metric, only one individual section was used for the analysis. This allowed use of a 
denser mesh with parabolic elements. For this initial stress analysis, the shroud was 
considered solid, and the strut was fixed at the mounting holes. 

Loads were defined as distributed pressure forces upon the flattened lower face of 
the shroud (i.e., the ERH thruster surface), the maximum thrust on each strut was 
found to be 4167 Newtons, which translates into a pressure load of 3.16 MPa across 
the shroud surface. 


5.10 Results of Static Analysis 

Of primary concern was the prediction of the stress fields throughout the strut struc- 
ture. The Von Mises stress was considered such that, for no plastic deformation: 

c o’vm = Von Mises Stress 
o'vm < y S y = Yield Strength of Nickel 200 


for no plastic deformations (failure). 
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Figure 5.8: Static Stress Analysis of Strut Assembly. 


Figure 5.9 displays the resulting stress contours for the given load case. The 
maximum stress was found at the outer edge of the strut to be 1457 MPa. The 
deflections under load were found to be significant. The maximum deflection at the 
trailing edge of the shroud was 4.625 Trim ; at the beam focus, the average deflection 
was 2.01 mm. 

As in most engineering analysis, the actual behavior of the strut system could 
only be approximated, due to inherent limitations of the approach. For example, in 
reality the trailing edge of the strut is attached to a carbon-carbon ring that encircles 
the LTD vehicle centerbody. However, for this analysis it was necessary to consider 
it discontinuous. Had a continuous structural ring been used, the bending moment 
imposed on the strut would have been significantly lowered. This, in turn, would 
have greatly reduced the stress concentration that appeared on the outer edge of the 
strut and the strut /shroud interface. 

Two very important results became evident from this analysis. First, the shroud 
section is very rigid and unstressed compared to the strut. This will permit removal of 
the excess, unloaded material at the center of the shroud; design mass limits can then 
be met. Secondly, the analysis shows that the current strut mounting hole pattern 
(for two bolts) is entirely capable of handling engine operating stresses. The critical 
stress concentration areas are located at the outboard edge of the strut, and at the 
point where the strut joins the shroud. 


I 
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Since the strut is loaded in alm ost pure bending, the stress concentration at the 
outer edge of the strut may be easily reduced by increasing its cross sectional area 
(i.e., moment of inertia, I). This can be accomplished by either thickening the strut 
or by increasing its height. While increasing strut height increases I faster than by 
increasing strut thickness, the former alternative will tend to lower modal frequencies. 
Thus it will be necessary to evaluate dynamic response as well as stress handling 
abilities. Finally, it has been proposed that the stress concentration at the intersection 
of the strut and the leading edge of the shroud will be reduced by adding a “stress 
reliever” . 


5.11 Conclusions of Static Analysis 

The static stress analysis examined the integrity of the strut/shroud section, and 
revealed areas that needed further refinement. The shroud itself was found to be 
excessively strong and can easily sustain the applied loads; thus, some material can 
be removed from the shroud interior to permit further weight reduction. Also, the 
two mounting holes provide sufficient strength for attaching the struts securely to the 
internal LTD support frame. Two areas are in need of further refinement: (1) the 
outer strut edge, and (2) the intersection of the strut and the shroud. A revised strut 
section with a higher moment of inertia and a stress reliever must be considered in 
future design and analysis. 
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Chapter 6 


OPTICAL SYSTEM FOR SENSING MISSION 


6.1 Introduction 

The Lightcraft Technology Demonstrator (LTD) represents not only a novel form 
of spacecraft propulsion, but also a different concept of launch vehicle integration 
with payload. Each part of the craft is designed with multiple functions in mind, 
and the same components may serve various applications at different times during 
the flight. For example, the liquid nitrogen is carried primarily as propellant for the 
laser heated rocket boost to orbit; however, a small amount will (i.e., 10%) also be 
used to cool both the primary mirror and the impulse thruster faceplates during the 
airbreathing propulsion mode. Also, where the cooling system is not regenerative, 
some liquid nitrogen could be vented over the impulse surfaces thereby providing a 
“film” cooling effect and perhaps helping to refresh the surrounding air before the 
next laser pulse arrives. It is this close integration of launch and mission components 
which makes the LTD so cost-efficient; this design approach is a logical extension of 
the so-called “lightsats”, a class of small, lightweight satellites designed to perform 
specialized functions [l]. 

One of the most expensive elements of the LTD design is the high quality primary 
mirror, which is essential to the launch phase. Once the vehicle has reached orbit, it 
seems reasonable to apply the onboard optical system to other mission applications. 
It has been proposed that the LTD could serve as an orbiting satellite telescope; 
a similar mission was proposed for the Apollo Lightcraft [2]. Once the LTD has 
established orbit, a sensor array (or retina) will be deployed using hydraulic arms as 
illustrated in Fig. 6.1. The telescopic arms will extend a sensor array to the focal plane 
of the primary mirror, which will measure 1 cm in height and extend around the full 
circumference (1.2 tt m) of the craft. The LTD may be used for imaging applications 
at a variety of wavelengths to observe ground targets, the upper atmosphere, or outer 
space. 

However, the main mirror has been designed primarily for propulsion applications. 
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Figure 6.1: Deployment mechanism for LTD retina. 

While the mirror is of large diameter and diffraction limited surface quality, the 
LTD telescope consists of a single refractive element with no corrective optics. The 
performance of this system is limited by geometric aberrations (reflective systems 
are free of chromatic aberration). The resolution limits of the primary optic under 
various conditions will be calculated, and the extent of significant aberrations will be 
evaluated. Several candidate devices for the sensor retina will be considered. 


6.2 Characteristics of the primary mirror 

The primary mirror is designed to be highly reflective in the infrared, although it 
should also exhibit good reflectivity in the visible, UV, or microwave regions of the 
spectrum. Multi-layer dieletric films may be designed for high reflectivity across a 
wide spectral range; thus, the mirror could serve as an antenna for microwave or 
laser communications. For such applications, the optical characteristics of the mirror 
become important. Several parameters of interest will be calculated, following the 
analysis given in a previous Apollo Light craft report [2]. The resolution of an imaging 
system, in general, is determined by the Rayleigh criteria [3]. The minimum angular 
separation between two points which can just be resolved is given by: 

A 6 = 1.22A /D (6.1) 
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where A is the operating wavelength and D is the input aperture diameter. The 
quantity (A0) is referred to as the resolving power of the system, while 2(A0) gives 
the diffraction limit angular spot size (to be compared with spot size in the presence 
of aberations). The Rayleigh criteria may also be expressed in terms of the minimum 
separation between two points which can be resolved at a distance l, 

d(l) = 1.22X1/ D (6.2) 

This quality is taken to be a measure of the system’s resolution limits when it is 
evaluated at l = f , the focal length of the system. 

Another important property is the numerical aperture of the system, N.A., which 
measures light gathering power [3]; it is related to the /-number (//#) by: 

//# = f/D = 1/2(N.A.) (6.3) 

The //# determines the amount of light which reaches the final image plane; its 
reciprocal is known as the relative aperture. A smaller //# allows more light to reach 
the final image; the / /# squared is a measure of exposure time. It is known as the 
speed of a system and smaller //#’ s imply higher speeds. 

Tjie LTD mirror has a focal length of 7.916 cm and is characterized by a numerical 
aperture of 3.158 and //0.158; the system is therefore very fast and exhibits high light 
gathering power, due to its short focal length. By comparison, the 200 in. diameter 
telescope at Mount Palomar Observatory has a focal length of 666 in., a numerical 
aperture of 0.150 and // 3.33 [4]. Similarly, one of the world’s largest refracting 
telescopes at Yerkes Observatory has a 40 in. diameter main mirror and focal length 
of 63 ft., yielding a numerical aperture of 0.026 and //18.9 [5]. However, the LTD will 
not out-perform either of these systems because it is poorly corrected for aberrations 
that limit image quality. 

Table 6.1 gives calculated values of resolution at various ranges for upper and 
lower limits of the visible spectrum, as well as two infrared wavelengths. If the 
system is to be used in flight, the image could be blurred by any mechanical vibration 
or flight turbulence greater than these values. In practice, resolution limits will be 
set not by the diffraction limit of the optics, but by the mechanical stability of the 
optical support structure. Pointing accuracy from orbit will be limited by the LTD 
attitude control system. The calculations of Table 6.1 do not account for geometric 
aberrations. In the following section, a general description will be given of the five 
monochromatic aberrations involved; a more detailed description may be found in 
Ref. [3,4]. Aberrations are more pronounced for rays far from the optic axis and 
for oblique angles of incidence. Coma and astigmatism are expected to dominate in 
the formation of a final image; the analysis begins with a description of the weaker 
aberrations. 
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Resolution limits for LTD mirror at various wavelengths 
(all angles in radians, all distances in meters) 




400 nm 

700 nm 

1 P 

10.6 p 

A6 

4.88 x 10-7 

8.54 x 10- 7 

1 .22 x 1 0' 6 

12.93 x 10*6 

d(f) 

38.63 x 10-9 

67.60 x 10-9 

96.57 x 10-9 

102.37 x 10-6 

d(36,000 km) 

17.56 

30.74 

43.92 

465.55 

d(360 km) 


0.17 

0.30 

0.44 

4.65 

d(180 km) 


0.08 

0.15 

0.22 

2.33 


Table 6.1: Resolution limits for LTD mirror at various wavelengths (all angles in 
radians, all distances in meters). 

Spherical aberration corresponds to a dependence of focal length on aperture size 
for rays far from the optic axis. This effect is well corrected by parabolic mirror 
sections [4,5]. Off-axis rays may also experience distortion, because transverse mag- 
nification is a function of distance from the optic axis. A ray trace indicates that 
the mirror has a very narrow angular field of view, so distortion is not expected to 
be a significant problem. Finally, the principal focal plane of the mirror is actually 
only a plane in the paraxial region; for off-axis rays, it takes on a finite amount of 
curvature. This is not significant because of the narrow field of view. It would be 
possible to correct for this effect by using a curved detector surface to match the 
wavefront profile, although our analysis suggests that this is not necessary. 

The dominant monochromatic aberration for the LTD mirror is coma, or comatic 
aberration. In an optic system, the effective focal lengths and the transverse magnifi- 
cations will be different for rays passing through off-axis regions of the system. This 
effect appears even for very small angles of incidence, affecting object points only a 
short distance from the optic axis. The formation of a comatic image is illustrated 
in Fig. 6.2; the final image forms a circular focus for paraxial rays, but off-axis rays 
form a comet-like tail in the image (hence the name coma). Slightly more than half 
the optical energy appears in the roughly triangular region between points 0 and 3. 

Coma is defined as positive if the marginal rays focus farther away from the optic 
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Figure 6.2: The geometric comatic image formed by an off-axis object point (Ref. 3). 

axis than the paraxial rays; if they focus closer to the axis, it is termed negative coma. 
A conservative estimate of the image spot size due to coma may be calculated from 
the angle of incidence, 0, and the primary mirror //# as follows [5] : 

fie = 0 / 16(//#) 2 ( 6 . 4 ) 

where f3 c is measured in radians. The comatic spot size for the LTD mirror at various 
angles of incidence is given in Table 6.2; for comparison, the diffraction limited spot 
size is also shown. Image quality will be limited by the aberration, even for a mirror 
with ideal figure and surface quality. The ray traces indicate that coma will be positive 
for one half of the LTD mirror and negative for the other. 

A detailed analysis of comatic aberration must account for the separate effects of 
tangential or meridonal coma (the distance between image point 0 and 1 in Fig. 6.2), 
and sagittal coma (the distance between points 0 and 3). [3]. The ray traces apply 
only to meridonal rays (a cross section cut through the primary mirror). From such 
traces, we may evaluate a meridonal error plot [4]; Fig. 6.3 illustrates the graphical 
interpretation of such plots. 

A representative meridonal error plot for the LTD is given in Figure 6.4. An 
equation which fits the curve and characterizes the aberrations is given; coma is the 
dominant effect. A 1 cm high detector array was assumed for the LTD retina design; 
this means that most of the rays will miss the detector plane entirely for incident angle 
much grater than 0.5 degrees, which is the upper limit on angular field of view. While 
a larger detector surface would intercept more rays at larger angles, the coma becomes 



98 


CHAPTER 6. OPTICAL SYSTEM FOR SENSING MISSION 


angle of incidence 
(degrees) 

coma spot size 
(radians) 

astigmatism spot size 
(radians) 

0.1 

4.37 x 10-3 

9.64 x 10-6 

0.5 

21.84 x 10-3 

0.24 x 10-3 

1.0 

43.67 x 10-3 

0.96 x 10*3 


Table 6.2: Spot size due to aberrations in coma at oblique incidence (wavelength = 
1 micron; diffraction limit = 4.88 x 10 “ 6 radians). 


dramatically worse and the final image would be distorted beyond recognition. Since 
additional corrective optics can not be included in the LTD mirror system, the field 
of view limit must be tolerated; it is assumed that the pointing systems for the LTD 
in orbit are sufficiently accurate to track a desired object within this limitation. 

Because of the limited capabilities of the ray trace software, no sagittal ray traces 
were performed and it was not possible to generate a sagittal error plot. However, 
theoretical limitations due to sagittal coma may be calculated; Kingslake [4] gives the 
following relation for sagittal coma as a function of off-axis distance, /i, for a small 
aperture mirror: 

coma, = /i/16(//#) 2 (6.5) 

Using this expression, a sagittal coma can be predicted for the LTD mirror of 50 mm 
for a point 20 mm from the axis. By comparison, the sagittal coma for the Mount 
Palomar telescope at the same distance off-axis is 0.115 mm [4], 

Despite its large diameter and high quality surface, the LTD does not perform as 
well as conventional telescope systems because of aberrations. It is possible to reduce 
coma by using an aperture stop at the focal plane [3]; the rays far off-axis are re- 
moved from the image. The image intensity will decrease sharply as the object points 
move off-axis; this is known as vignetting [3]. Since the LTD retina is composed of 
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(a) an ideal optical system 

(b) inward-curving field 

(c) positive coma 

(d) negative coma 

(e) spherical aberration 


Figure 6.3: Sample meridonal error plots (Ref. 4). 



□ ray height.cm 


Tan. angle in image space 


Figure 6.4: Meridonal error plot for right side of the LTD primary mirror at 0.1° 
angle of incidence. 
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electronic detectors, the sensors could be apertured electronically to reduce coma; 
sensitive detectors could decrease the effects of vignetting. Electronic image process- 
ing techniques may significantly increase the fidelity of images generated by the LTD; 
aberrations could be compensated by electronic post-processing of the image. 

Astigmatism occurs when an object point lies a considerable distance from the 
optic axis, so that a cone of incident rays strikes the mirror asymmetrically. The 
tangential and sagittal rays focus at two different points; the separation between 
these points is known as the astigmatic difference. The spot size of an image due to 
astigmatism may be calculated [5] from the relation: • 

0 . =« ! / 2 (//#) ( 6 - 6 ) 

where 3 a is measured in radians. Values are given in Table 6.2; astigmatism is only 
pronounced for large angles, and coma still dominated image formation over the entire 
field of view. 

The proposed mirror design represents a compromise between several possible 
mirror contours, and was selected for its propulsion applications. Different contours 
for the parabolic mirror will possess slightly different focal lengths [6], as shown 
in Fig. 6.5. The //# for each of these designs is given in Table 6.3, which also 
shows the calculated spot sizes due to coma and astigmatism. Both aberrations may 
be improved by using larger //# mirrors; the remaining aberrations would not be 
substantially affected. It is feasible to consider designing customized LTD satellites 
which trade off propulsion efficiency for improved imaging performance on certain 
mission applications. 

The ray tracing program defines a figure of merit called root mean square (rms) 
blur. It is defined as the perpendicular rms distance between the mirror focal point 
and the intercept of actual rays with the focal plane. The rms blur is thus an estimate 
of the final image spot diameter. The normal distance from each ray to the focus and 
the corresponding rms blur is given in the Appendix; the blur increases linearly with 
incident angle. 


6.3 Sensors for the LTD retina 

Current designs call for a 1 cm high sensor retina to be deployed from the LTD 
forebody to the focal plane of the primary mirror. Several different types of sensors 
could be deployed around a single vehicle, enabling it to perform multiple missions. 
The devices would be lightweight, rugged, compact, and require minimal power; this 
suggests a solid state detector array. Such a device must be sensitive enough to ob- 
serve faint, more distant objects or scan brighter objects in less time. This translates 
into high quantum efficiency (ratio of electrons generated to number of incident pho- 
tons) combined with a low readout noise level. The detection of faint objects can be 
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Figure 6.5: Alternative primary optic contours and focal lengths. 



Table 6.3: Aberration spot size for various primary optic contours at incident an- 
gle = 0.5° (diffraction limit at 1 micron length = 4.88 x 10 -6 radians). 
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enhanced by increasing the exposure time, assuming that the mirror system could 
be held sufficiently stable; this technique is also used to observe video images of a 
changing scene. Electronic detectors generate noise over prolonged exposures. The 
flow of thermally generated electrons (dark current) saturates the detector over time, 
even when no light is present. Low dark current is thus an important requirement; 
since dark current is smaller at reduced temperatures, cryogenically cooled detectors 
are sometimes used (especially for infrared sensors). 

Solid state detectors saturate when exposed to very bright objects; the optically 
generated signal overflows into adjacent pixels, causing an effect known as blooming. 
The LTD sensors will be exposed to space with very little shielding; cosmic rays or 
high energy particles may pass through the detector and blind large areas of the retina. 
Other factors to consider include wide spectral response and a large number of pixels 
for high resolution. The sensors should collect as much raw data as possible, since 
many of the details could be recovered by electronic post-processing of the image; it 
should also be possible to electronically aperture the detector, in order to concentrate 
on details of a scene or reduce coma. 

Only commercially available sensors such as the charge coupled device (CCD) 
and charge injection device (CID) [7,8] will be considered. Charge coupled devices 
represent a mature technology for space systems applications; a CCD system was 
included as part of the optical instrumentation on the Voyager spacecraft. The same 
company that developed these sensors, EG&G Reticon, has recently demonstrated a 
new CCD detector for astronomy applications [7]. It consists of a 400 x 1200 matrix 
of pixels, each measuring 27 microns square. Readout noise has been reduced to 
4 electrons rms. Dark current is very low, only 20 electrons per pixel per hour at 
- 100 0 C. Operating unshielded in space, the detector could reach low temperatures 
which would improve performance. The device has superior anti-blooming perfor- 
mance; only three adjacent pixels saturate when enough light is applied to generate 
1.5 million electrons (over three times the saturation level). Figure 6.6 shows the 
effect of a cosmic ray passing through pixel 792 of the device; only a few adjacent 
pixels are blinded. The number of electrons generated by the interaction is plotted 
for each column of the device. The temperature of the device was -111 0 C (ref. 7). 
The spectral response extends from 400 to 1100 nm. Small portions of the device can 
be scanned independently, making it possible to aperture the device. 

Another candidate device offered by CIDTEC Corp. [8] does not transfer accumu- 
lated charge from site to site. Collected charge is shifted between capacitors within 
each pixel, and the displacement current (which is proportional to the stored signal 
charge) is detected. Unlike the CCD, which clears itself whenever a new image is read 
out, the charge remains intact within a CID detector’s pixels so that readout is non- 
destructive. In order to clear the array for a new image, the pixel electrodes switch 
briefly to ground, “injecting” the stored charge into the device substrate. The nonde- 
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Figure 6.6: Effect of a cosmic ray on Reticon CCD. 

structive readout capability makes it possible to control a time exposure accurately, 
and to check the image as it develops without erasing it. 

Since charge does not transfer between individual pixels as in a CCD device, there 
axe no paths between pixels along which charge overloads can propagate. This makes 
the device very tolerant to blooming, since illumination overloads are confined to each 
pixel and excess charge is drawn off into the substrate; cosmic rays only blind indi- 
vidual pixels. Because there is no need for opaque electrical contacts between pixels 
to transfer charge, the CID has a continuous pixel structure for higher resolution. 
Since each pixel can be addressed individually, electronic apertures are possible. A 
signal to noise ratio of 46 dB has been achieved, and detectors of512 x776 pixels are 
available in sizes suitable for the LTD. 


6.4 Conclusions 

The use of the LTD as a sensor satellite has been explored. The primary mirror 
could serve a dual purpose on each mission; however, since the mirror is designed 
primarily for propulsion applications, its performance in other respects is limited. A 
ray trace analysis of the primary optic indicates that it has superior light gathering 
ability, but the angular field of view is limited to only about 0.5 degrees. Coma is 
the dominant aberration over the entire field of view; a small amount of astigmatism 
is present for the larger angles of incidence, while other aberrations are essentially 
negligible. Several candidate sensors, including CCD and CID technology, are com- 
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mercially available and meet the requirements of the LTD retina. 
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7.1 Introduction 


The following chapter describes design considerations for the Lightcraft Technology 
Demonstrator (LTD) propulsive optical train. Design and fabrication techniques axe 
considered for the LTD primary optic (PO), as well as characteristics of the ground 
based laser uplink. Using a ray trace of the PO, the effects of beam steering error on 
the propulsion system performance are evaluated for near, intermediate and far field 
intensity profiles. A proposal for the control system necessary to accurately track 
pointing error will be given. 

The LTD represents a novel concept in satellite self-deployment, although laser 
powered launch vehicles have been under study for several years. Small, lightweight 
(ten to several hundred kg), unmanned missions for specific applications have re- 
ceived increased attention as both equipment and launch costs continue to rise [1], 
These so-called “lightsats” or “microspacecraft” have been considered for many di- 
verse applications including communications satellites, lidar systems, and satellite 
laser rangefinding, among others [1,2]. As microelectronics technology continue to 
reduce payload size and weight, lightsat systems have been recognized as the means 
to fill a vital market niche in space science [2]. Laser propulsion itself has also received 
increased attention; NASA has even considered manned spacecraft based on huge (10 
GW) orbital laser power stations (the Apollo Lightcraft [3]). Lawrence Livermore Na- 
tional Laboratory (LLNL) is currently engaged in a program to develop the LTD for 
launch within the next six years, funded under the SDIO Laser Propulsion Program. 
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7.2 Design and manufacture of the LTD primary 
mirror 

The design of the 1.0 m diameter primary mirror is critical to the successful launch of 
the LTD. This mirror will be an off-axis paraboloid with a short focal length, designed 
for high light gathering capability. Such mirrors are co mm only made by cutting 
several off-axis sections from a large on- axis parabola mirror blank; however, the cost 
increases with the off-axis mirror diameter [4,5]. Since the LTD calls for a mirror cut 
45 in. off-axis, this technique would be prohibitively expensive, and other methods 
must be employed. Large diameter, high quality mirrors may be manufactured by 
diamond-turning, a technique using a specially designed high precision lathe with 
a single crystal diamond cutting tool. Computer controlled milling procedures are 
employed to machine certain metals (e.g., aluminum, copper, etc.) to exceptionally 
close dimensional tolerances. Pioneering efforts at LLNL have developed a 96 in. 
diameter di am ond- turning lathe, which would be able to cut the LTD primary optic 
in a single processing step. An aluminum substrate would be used, because of its 
machining properties and its high reflectivity in the infrared. The LLNL equipment 
has already produced annular resonator optics for the Alpha ground-based chemical 
laser experiment to a surface figure accuracy of 12.5 nm rms [6]. The design of the 
LTD mirror is shown in Figure 7.1. 

In order to achieve the necessary optical quality, the mirror figure (large-scale 
flatness) should be kept to at least 1/10 of the operating wavelength. Reflectivity 
exceeding 99.9% must be maintained due to the high laser intensities required during 
launch; this cannot be achieved with bare metal mirrors because their intrinsic ab- 
sorption is too high. Instead, high reflectivity multi-layer dielectric coatings must be 
applied over the mirror substrate, as illustrated in Figure 7.2. While these coatings 
are currently available, they set limits on the surface quality which can be realized. 
For example, a figure of 1/40 wavelengths requires coatings with thickness uniformity 
of 0.1 nm, which is less than the diameter of an atom. The outer surface roughness 
should be kept within 1 nm rms to minimize scattering of the incident beam. 

Another important issue is the stability and survivability of such a precision op- 
tic during a hostile launch and orbital environment. Since mirror systems can be 
supported from behind (as opposed to transmissive optics, which must be supported 
only along the edges), this could be achieved by using a carbon composite backing 
for structural support. The high intensity laser tends to produce isolated spot dam- 
age in the optical coatings; laser heating also causes thermal distortion in the mirror 
figure, and accelerates mechanical creep effects [7]. To address these problems, heat 
exchangers may be incorporated in the mirrored face plate, as shown in Fig. 7.3. A 
face plate on the order of 30 mils (0.5 mm) thick would incorporate coolant channels 
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1/4 

wavelength 


Multi-layer dielectric coating for high reflectivity; thickness of 
each layer is one-quarter wavelength of the illuminating light, 
so that light is reflected in phase from all interfaces in the 
coating. Dark layers represent higher refractive indices, light 
layers represent slightly lower refractive index 

Figure 7.2: Multi-layer dielectric coating for high Reflectivity. 

about 20 x 60 mils square. The performance of these channels at high pressures and 
flow rates is still under investigation, although the theoretical potential exists for an 
active cooling system to improve the mirror performance. 

This same cooling system would also circulate coolant behind the shroud impulse 
surface, where laser intensities of 500 MW/cm 2 must be maintained for useful thrust 
generation. The system would be largely regenerative, but could also used for tran- 
spiration cooling at the shroud and support strut leading edges. Also, a small amount 
of coolant could be vented over the impulse surfaces between laser pulses to provide 
a film cooling effect and perhaps assist in the convection of refresh air [8,9]. 

The mirror surface could be protected from dust and impacting particles with a 
severed hundred angstrom coating of diamond-like carbon; this coating would also 
discourage thermal damage, and act as an electrical insulator to shield the metallic 
mirror from high electric fields of the laser beam [10]. Diamond coatings are also 
highly transmissive in the infrared. 

Research into the physics of high energy beam interactions is continuing, and may 
have implications for the LTD experiments. For example, intense laser interactions 
with metal targets may lead to photoelectric effects. This would cause the PO to 
steadily lose electrons during flight and accumulate a net positive charge, which could 
attract small particles damaging to the mirror surface. Assuming single photon /atom 
interactions to be the dominant mechanism, the critical wavelength of light required 
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Figure 7.3: Proposed models for mirror face cooling channels. 

for photoelectric emission is given by [11]: 

A = hcjip (7-1) 

where Plank’s constant, h = 4.136 x 10~ 15 eV/s, c is the speed of light, and if is 
the work function of the metal target. Since the mirror substrate may consist of 
aluminum, whose typical work function is 4.28 eV, photoelectric emission does not 
occur below wavelengths of 0.29 microns. Thermionic emission of electrons, however, 
does not depend on the incident wavelength but on laser heating of the mirror; the 
total current per unit area generated, J, is given by [11]: 

J = 2AT 2 e~' fi/kT (7.2) 

where T is the temperature of the target in Kelvin, Boltzman’s constant, k = 1.3805 x 
10 -23 J/K, and the constant A = 120.4 amp/(cm 2 deg 2 ). For a worst case analysis, 
one could take T = 600° C, which is just below the melting point of aluminum; for 
this case, the emitted current density is 17.9 x 10 -18 amps/cm 2 , which is negligible 
for our purposes. 

It has also been reported [7] that multi-photon interactions may play a significant 
role in high energy laser interactions. For laser wavelengths of 1 micron, a single 
infrared photon energy is 1.24 eV; thus, a 4 photon interaction is required to overcome 
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the work function and produce electron emission. While this is highly unlikely in the 
infrared, recall that approximately 50 x lO 30 photons /laser pulse are incident on the 
mirror during flight, so multi-photon interactions may become significant. Recent 
experiments [12] have demonstrated 7 to 12 photon interactions with xenon gas at 
wavelengths of 10.6 microns and beam intensities of the order 10 13 W/cm 2 . This 
suggests that multi-photon absorption during atmospheric propagation of the power 
beam may also be significant. Both of these effects are dependent on factors such 
as the laser power beam intensity profile, incident angle, and polarization; thus, it is 
expected that such factors could be controlled in the LTD system. 


7.3 Power beam system analysis 

Several transmission models have been proposed for free space high power laser prop- 
agation [3,131. To be presented here is a first-order analysis for the LTD system based 
on some aspects of these models, assuming a ground-based laser and using the pri- 
mary optic of the LTD as the receiver antenna. The power beam diameter will spread 
as it propagates through the atmosphere; the significant effects have been identified 
as diffraction from the source aperture, atmospheric turbulence, pointing jitter, and 
thermal blooming [13]. In practice, the effects of thermal blooming can probably 
be reduced or eliminated by using adaptive optics at the transmitter [13] and the 
judicious choice of a laser frequency that coincides with an “atmospheric window 
The three remaining sources are expected to contribute almost equally to total beam 
spread. The maximum tolerable beam spread, <7, (in steradians) is given by: 

2 a=r/R (7.3) 

where r is the radius of the receiver aperture (0.5 m) and R is the range between 
transmitter and receiver. The maximum effective range encountered by the LTD 
depends on the path taken to low Earth orbit; several possible trajectories have been 
calculated [3, 8, 9]. Assuming a trajectory of 30° inclination with respect to the 
horizontal, the maximum range is 560 km. This leads to a maximum tolerable spread 
of 446.42 x 10~ 9 steradians, if we are to prevent excess spillage of light around the 
vehicle. The contributions of diffraction, turbulence and jitter are all approximately 
258 x 10 “ 9 steradians. 

Given the maximum beam spread due to diffraction, we may estimate the size of 
the transmitter aperture from: 

transmitter diameter = 0.45A/9/cr (7-4) 

where A is the optical wavelength and /? is the beam quality (assumed equal to 1 
for the diffraction limited case). The minimum required transmitter diameter is thus 
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Figure 7.4: Graphic map of laser upink propagation parameters (ref. 3). 


1.74 m. The relationship between transmitter and receiver size depends in detail on 
factors such as the operating wavelength, beam quality, orbit range, and tolerable 
spread due to jitter and turbulence; these relationships have been calculated [13] and 
are displayed graphically in Figure 7.4. An example showing how a given transmitter 
aperture and wavelength might map into different receiver diameters at different 
ranges is given (diffraction limited performance is assumed). In a similar manner, it 
is possible to map a 1.7 m diameter source into a 2 to 3 m diameter beam in low Earth 
orbit, depending on the assumptions made regarding jitter and turbulence. Current 
SDIO plans involve the construction of a 4 m diameter adaptive primary mirror for a 
ground based telescope system which would also track a laser boosted launch vehicle; 
an adaptive 4 m mirror recently produced by Itek Corp. [14] may be considered as a 
point design for the transmitter. 

For the present LTD vehicle, the propulsive power beam is assumed to deliver 
severed hundred megawatts; an ultralight LTD would require 10 times less power. 
These figures represent time-average levels, including both the pulse duration and 
the interval between successive pulses. Each laser pulse is expected to last only 25 
- 50 ns, with an energy of 30 - 40 kJ per pulse. A square pulse is assumed at the 
transmitter; for this case, each pulse is transmitted peak power, P, of 1 x 10 12 W. 
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The intensity received at the LTD primary mirror, I , may be calculated from [13], 

I = PTJ2tt{(tR) 2 (7.5) 

The quantity T r is the atmospheric transmittance, which is about 0.85 for an operating 
wavelength of 1 micron [13]. This yields an intensity of 2.16 x 10 W/cm at the 
primary mirror. Assuming a pulse duration of 35 ns, the mirror must be designed to 
tolerate a fluence level F of 

F = /At = 7.57 x 10 4 J/m 2 (7.6) 

Although this is higher than currently available mirrors have obtained, the use of 
cooled mi rrors should make it possible to achieve this level of fluence without inducing 
substantial laser heating damage to the main mirror. * 


7.4 Effects of pointing error 

In order to determine the performance of the LTD receptive primary mirror under 
v ari ous conditions, a computerized ray trace study of the primary mirror has been 
generated, as described in the Appendix. Characteristics such as the mir ror figure 
have been estimated previously; the primary purpose of this analysis was to determine 
the effects of pointing error and oblique incidence on the propulsive laser beam. Three 
possible beam intensity profiles were considered: near, medium, and far field. Actual 
ranges corresponding to these three regimes are not given, since they would depend 
on the specific trajectory to orbit [3, 8, 9]. Since the laser will be directed by the 
laxge aperture telescope, with the minor beam steering performed by adaptive optics, 
the near field irradiance is taken to be a torus. The central obscuration due to a 
secondary telescope optic is assumed to be of maximum radius 50 cm, which would 
obscure the central 25% of the PO area. Not performed, was a detailed fractional 
active area analysis for the LTD, as reported for the Apollo Lightcraft [3], because of 
the small diameter primary optic and narrow field of view permitted. The remainder 
of the toroidal beam is assumed uniform in intensity. 

At intermediate ranges, this central portion fills in by diffraction to yield a uniform 
or ‘top hat’ intensity distribution. This case may also be encountered in the near field 
of a telescope with no central obscuration; the total power delivered is assumed to 
be the same for both toroid and uniform distributions. In the far field, the beam 
assumes a Gaussian profile, which is assumed to be truncated when the intensity in 
the Gaussian tails has fallen by e ~ 2 (or 0.1353) from its peak value [13]. This allows 
for 95% of the power beam to be intercepted by the LTD mirror. A small amount of 
light spillage is thus tolerated, and may even prove to be useful as shall be illustrated 
later. 
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If there is a small pointing error in the power beam, the rays focus at a differ- 
ent point from the ideal focus; some rays miss the focus entirely for large angles of 
incidence. This corresponds to a loss of incident optical power. The spot size near 
the focus spreads out over a larger area, reducing intensity at the focus. These two 
effects combine to decrease the effective thrust generated by the LTD; the effect is 
worse for the larger incident angles and for the rays further off the parabolic mirror 
centerline. The effects of oblique incidence are shown in the Appendix. Each of these 
figures represents the case of uniform illumination; in the near field, rays 1-5 would 
not exist. In the far field, rays 1-5 would carry over 68% of the incident power, while 
the remaining rays convey about 27% (recall that about 5% beam spillage around 
the ship is permitted). The intensity at the focal point may be calculated at differ- 
ent angles for each of the three intensity distributions; the results are summarized in 
Table 7.1. Note the different intensities generated on either side of the craft. For the 
cases considered, a sufficient intensity to provide useful thrust was always maintained. 
The focal point moves in opposite directions on either side of the ship for small angles. 
This creates asymmetric thrust, which tends to push the craft back towards normal 
incidence; thus, for errors of 0.5 degrees or less, the LTD remains stable. However, 
it is still important that the LTD detects errors of this magnitude (or greater) and 
takes corrective action. 

A system to sense and correct for large pointing errors is also proposed for the 
LTD. Even for the torus and uniform intensity distributions, a small amount of light 
leaks around the vehicle due to scattering; in the far field, the light in the Gaussian 
tails of the pulse becomes significant. A set of comer cube retroreflectors could be 
placed around the shroud trailing edge to intercept this light and reflect it back to 
the source. A detector array distributed around the transmitter source aperture (i.e., 
on the ground) could sense these reflected signals, and the distance to different parts 
of the shroud could be calculated by laser rangefinding techniques. Adaptive optics 
at the source would then be able to compensate for pointing errors before subsequent 
laser pulses were d eli vered to the vehicle. A total of 24 retroreflectors could be placed 
on the LTD (one for each shroud support strut), each with an area of about 1 cm 2 . 
The maximum laser intensity assumed to exist at this point is about 6.08% of the 
peak intensity in the far field; this should not damage the highly reflective comer 
cubes, and still provides a sufficient signal for rangefinding. If the ranging signal 
proves to be too weak in the near or medium field, the power beam can be tailored 
to spill slightly more light around the LTD for this purpose. Typical satellite laser 
rangefinding systems exhibit range accuracies of 1 cm. If the retroreflectors are placed 
at the shroud edge, 0.7 m from the vehicle centerline, then a difference of 1 cm between 
either side of the ship corresponds to a pointing error of 0.477°, which is adequate 
resolution for our purposes. 

Assuming the same maximum range used in the beam spread calculations, the 
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power beam would take 1.867 ms to reach the vehicle at maximum range. This 
provides a round trip time for the ranging signal of 3.733 ms. If the laser power 
beam is assumed to operate at 100 Hz, this provides 6.267 ms between successive 
laser pulses to adjust parameters of the beam. Recall that current designs call for 
a 4 m diameter transmitter aperture, which is estimated to require of the order 10 
individual adaptive elements, each of which could be adjusted at rates of 1 to 3 kHz 
[3]. Accounting for the actuator response time, this leaves between 5.267 and 5.937 ms 
for the rangefinding and beam control systems to operate. This is within the realm 
of adaptive computer control systems; however, if the laser power beam operates at 
1000 Hz the control system only has time to adjust every 4 th pulse. This time is 
set by the round trip propagation time for the rangefinding signals, and cannot be 
overcome by more efficient adaptation algorithms or faster actuators. However, an 
adaptive computer controller might estimate the correction required for successive 
(i.e., future) laser pulses; this would allow the continuous adjustment of the source, 
with fine adjustments provided by rangefinding data whenever it was available. 


7.5 Conclusions 

‘ This chapter describes the design and fabrication of the LTD primary receptive mirror. 
Although some of the physics behind propagating the high energy laser beam through 
the heated exhaust plume is not fully understood, the concept is shown to be within 
reach of current technology. The primary optic could be manufactured using existing 
diamond-turning lathe facilities; development of multi-layer high reflectivity coatings 
and microchannel cooling systems are the main technical problems encountered. The 
mirror must exhibit high figure quality and minimal surface roughness throughout 
the launch trajectory. Calculations of beam divergence indicate that a 4 m diameter 
adaptive ground based transmitter is entirely sufficient for a 1.0 micron wavelength. 
The effects of pointing error were evaluated using a computer ray trace routine; errors 
of 0.5 degrees, or less, shifts the focus asymmetrically on either side of the annular 
engine, but the effect is self- correcting and the LTD vehicle is therefore stable. At 
incident angles up to 1 degree, the power loss and spot size increase but do not cause 
a critical loss of thrust for near, intermediate, or far field beam profiles. An adaptive 
pointing error correction system is proposed, using retroreflectors to sample the power 
beam edges at 24 azimuthal locations and feed data back to the ground power source 
via laser rangefinding. 
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Chapter 8 


MODELING OF LASER-GENERATED 
IMPULSE 


8.1 Introduction 

The focus of this chapter is to model the laser-heated blast wave and impulse- 
generation process using an axi-symmetric representation, that excludes radiation 
and conduction. Chapter 9 will review the particular methods used to remove some 
of the more restrictive assumptions of this case, as applied to the Lightcraft Technol- 
ogy Demonstrator engines in the airbreathing propulsion mode. 

For the most part, air within the laser-generated blast waves is modeled as an 
ideal gas, throughout the expansion process. Note that the initial state (i.e., pressure, 
temperature, density, etc.) prior to laser- heating, is a direct function of flight Mach 
number and altitude (see chapter 9); hence, the graphical results in this chapter are 
also presented versus flight Mach number and altitude. 


8.2 Initial State of an LSD Wave 

The equations developed by Raizer [1,2] to describe the state of the laser-induced blast 
wave, were used in modeling the “line-source” blast wave expansion and impulse- 
generation process. The derivation is briefly presented below. 

Raizer begins with the control volume in Figure 8.1, and assumes that the laser 
supported detonation (LSD) wave moves steadily towards the left at constant strength 
with no losses. He then applies the mass, momentum, and energy conservation equa- 
tions in their integral forms, as reproduced below: 


PqVlsd = P u 


(8.1) 


Po + PqVlsd — P + pu 


( 8 . 2 ) 
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Figure 8.1: The control volume used by Raizer in modeling an LSD wave. 


PO V LSD + ^LSP 

p o 2 PqVlsd 



(8.3) 


where 4>lsd IS the laser radiation intensity, Vlsd IS the LSD wavefront velocity, u is 
the gas velocity behind the LSD wavefront (measured relative to that front), p is the 
pressure of the fluid, and p is the density. Also, e represents the internal energy of the 
gas, which is defined as e = RT / (7 — 1) = p/[(7 “ l)p]> where R is the gas constant 
and T is the temperature (i.e., ideal gas assumption). The subscript ‘O’ denotes the 
ambient values of a given quantity. These equations may also be obtained in other 
ways (i.e., see Ref. 3). 


The shock adiabat is then found to be: 



Using the equation of state, (p = pRT ), an explicit expression for the shock adiabat, 
p(p), of the wave can be obtained: 



(8.5) 
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30 km 
20 km 


Flight Mach number 

Figure 8.3: Peak LSD wavefront velocity vs. flight Mach number and altitude. 


For a full discussion of the shock adiabat plotted in Figure 8.2, see Raizer [2]. 

The curve B-C-0 in Figure 8.2 represents Eqn. 8.5. The point J gives the pressure 
for a strong detonation wave, according to Chapman- Jouguet theory (see Thompson, 
Ref. 4). The point C represents a deflagration wave. LSD waves are expected to be 
strong detonation waves. 

With the CJ values much greater than ambient, and c, CJ = y^TPj p)cj = V LSD ~ 
v rCJ , it is found that c, CJ = 7 v rCJ . However, with 7 close to 1, c tCJ = v rCJ so that 
c *cj = Vlsd! 2, in agreement with Chapman- Jouguet theory (i.e., a resultant local 
Mach number of M C j = 1). Also, the velocity of the wavefront is found to be: 


Vlsd — 


2(7 2 -l) 


$ 


LSD 

Po 


and the pressure directly behind the LSD wavefront is: 

P oVIsd 


Vdw — 


7 + 1 


(8,6) 


(8.7) 


The LSD wavefront Mach number (through ambient air, and right after initia- 
tion), can be defined as M LSD = V LSD / ■sj') RT amincnt , using the local ambient gas 
temperature. This Mach number can be as high as 250. 
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Figure 8.4: CJ pressure vs. flight Mach number and altitude. 
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Figure 8.5: CJ temperature vs. flight Mach number and altitude. 
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While plasma temperatures can be very high, these conditions exist only for cycle 
times on the order of a millisecond. Peak plasma temperatures occur over shorter 
time intervals within small volumes that are somewhat removed from material sur- 
faces. Due to the pulsed nature of the impulse generation process and the feasibility 
of including active cooling systems for material surfaces, these thermal problems ap- 
pear surmountable. This report does not examine extemal-radiation-heated (ERH) 
thruster heat transfer characteristics, but such aspects certainly need further study. 

It is noted that Raizer’s equations do not incorporate real gas behavior, which 
is complicated because of the dissociation and ionization effects. Hence the extreme 
plasma temperatures predicted in Figure 8.5 would never be realized; peak values of 
10,000 K to 30,000 K are more in agreement with reality. 


8.3 Application *of the Method of Characteristics 


The method of characteristics may be used with Raizer’s equations to describe the 
blast wave expansion process. The following is based on such a model as similarly 
discussed by Holmes et al (Ref. 12), and is covered in greater depth by Liepmann 
and Roshko [5]. 


The key equations are: 


dp 

dt 


dp dv pv T 

+ u T — + P~K~ + = 0 

dr dr r 

(8.8) 

dv T dv T 1 dp 

dt JrVr dr p dr 

(8.9) 

II 

fU q* 

(8.10) 


which are basically the same expressions presented earlier. Information from the 
energy and state equations is is contained in Eqn. 8.10. As before, ideal gas behavior 
and isentropic flow conditions are assumed. 

If the radius, r, of the shock front is much greater than its thickness, the wave can 
be assumed to move out with a planar front; hence, the 1/r term can be neglected. 

If Eqns. 8.8 and 8.9 are multiplied by the speed of sound, c,, then Eqn. 8.10 can 
be used to remove p from these expressions. After addition and subtraction of the 
resulting equations one obtains: 


which is equivalent to 


8 , . d 

— + K ±c,)— 


d , ^ ,8 

rn +{v ’ ±c ‘h r 


v T ± c. 


V r ±2 


O 

II 

* 

U 

(8.11) 

1 c ’ \ = 0 
7-1/ 

(8.12) 
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using the ideal gas assumption of 7 = constant . 
From this, it is reasonable to write 


d ;l±2 \)=„ 

Dt \ 7 ~ 1 / 

(8.13) 

as a differential if 

dr 

Tt =v ’ ±c - 

(8.14) 

From these relations, it can be said that 


c, 

v T ± 2 = constant 

7 - 1 

(8.15) 

along the paths 

dr 

Tt =v ’ ±c - 

(8.16) 


Equation 8.15 yields the Riemann invariants, so called because these do not vary 
along the characteristic paths given by Eqn. 8.16. The Chapman-Jouguet condition 
which detonation must also satisfy, includes the requirement that (v r -t-c,)^ = V^sd 
(see Courant and Friedriechs [6]). 

Thus, the characteristic path is the path of the shock; The ” characteristic 
is the path of a rarefaction wave which travels at the velocity, v r — c,. The latter 
carries information from the blast wave initiation point towards a direction opposite 
that of the shock, if v T < c,. 

As illustrated in Figure 8.6 (from Liepmann and Roshko), both shock and rar- 
efaction waves start out with a high pressure region bounded at, say, r = 0 and t = 0; 
the gas starts to move for t > 0. 

Holmes uses the above results to calculate the pressure at the surface, where 
v T = 0. The “surface” pressure has been misleadingly denoted by many researchers 
as Plsd i which is different from the other parameters associated with CJ conditions 
and using the subscript LSD . The actual surface pressure is found to be: 


Plsd = 




Pdw 




Pcj 


7 + l \ T " 1 PqVlsd 
27 / 7 + I 


(8.17) 


which is plotted in Figure 8.7. 

Barnard and Bradley [7] have noted that real chemical detonations depart from 
ideal theory in that pressure, temperature and density distributions may appear as 
in Figure 8.8 (planar detonation wave front). Note that the gas conditions remain 
constant over a short induction period immediately behind the detonation front. For 
an LSD wave, it is possible to have the pattern pictured in Figure 8.9. 
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Figure 8.6: Shock tube gas dynamics. 
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Flight Mach number 

Figure 8.7: Peak surface pressure, before cylindrical blast wave decay. 


8.4 Sedov’s Solution of the Equations 

The self-similar solutions for the above-mentioned differential equations (under axi- 
symmetric and isentropic assumptions) were analyzed by Sedov [8] and also reviewed 
by Whitham [9], These solutions are applied in the laser-generated blast wave and 
impulse model for the LTD engine. 

The term “self-similar” describes a type of flowfield in which the value of a certain 
dependent variable (e.g., p, p, T, v r , etc.) is scaled by its value at some specific 
reference point. This reference point may coincide with the given initial conditions, 
or another point where the value must be determined. In any case, a relationship 
must exist between the dependent variables such that the flow parameters exhibit a 
similar pattern along the flow field. 

The method of characteristics is used only for that portion of the blast wave 
evolution which is assumed planar (for t % f p ); the Sedov equations are applied 
afterwards, when the blast wave has evolved into a cylindrical shape. The method of 
characteristics is unnecessarily complicated for the simple one-dimensional problem 
being studied, since one must track the rarefaction and compression waves with this 
approach. Although much could be learned about the blast wave structure using the 
method of characteristics, only one specific item (surface pressure) is pursued here; 
full knowledge of the effects of the rarefaction waves’ motions and reflections must 
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Figure 8.8: A real chemical detonation (Ref. 7). 
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Figure 8.9: Dynamics of a propagating Laser Supported Detonation. 
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await a future study. 

Sedov had also considered the case of a one-dimensional, axially propagating blast 
wave and a radially propagating, sphencally-synimetric blast wave. These additional 
cases are mentioned here for comparison; Sedov’s scaling laws are all presented in 
Table 8.1 (Reilly [10], and Sedov [8]). Sedov is generally recognized as the first to 
derive closed-form, self- si mil ar solutions for the blast wave problem. 

Because of the scaling relationship used in self- similar solutions, it is assumed 
that variables can be arranged into a non-dimensional form. Sedov s approach uses 
dimensional analvsis. He claims that since the dimensions of p and p contain mass , 
there is at least one constant, a, whose dimensions also contain mass . This constant 
must then be a characteristic parameter of the problem. Hence, the arrangement of 
the dimensions for a can be arbitrarily chosen, and Sedov assumes 


[a] = ML k T 9 


(8.18) 


This is a definition which outlines o s dependence on mass . Here mass is is 
denoted by M, T is ‘time’ and L is ‘length’. The exponents k and s are to be 
determined. 

The dependent variables are non-dimensionalized by setting 


V r 


-V 

t 


P = 


r k+$t 


~P , V = 


r fc+ 1£«+2* 


(8.19) 


with the non-dimensional, ‘prime’ variables as yet unknown. However, they are gen- 
erally dependent on two non-dimensional variables. 

If another characteristic parameter b is introduced with dimensions independent 
of a, the number of independent variables which can be formed by combining a and b 
is reduced to one. Since the dimensions of a include M, it can be excluded from the 
units chosen for i; hence, 

[I b } = L m T n (8.20) 

The single non-dimensional independent parameter in this case is thus r m t n /b 
which for m/0 can be written as 


c = 


r 

b m t& 


( 8 . 21 ) 


where f3 — — ti/ttl . 

In generalizing his analysis, Sedov writes Eqns. 8.22, 8.23 and 8.24 as 

at or r 


dv r dv r _ 1 dp 

dt r dr p dr 


( 8 . 22 ) 

(8.23) 
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d_ 

at 




= o 


(8.24) 


where v = 1 for plane flow, v = 2 for flow with cylindrical symmetry, and v = 3 for 
flow with spherical symmetry. 

Sedov claims that the solutions of Eqn’s 8.22-8.24, may contain a number of 
arbitrary constants, depending on the independent variable. When the characteristic 
parameters include two constants with independent dimensions (in addition to t and 
t), he shows that the partial differential equations for v r , p and p can be replaced 
by ordinary differential equations for v ' T , p ' and p'. The solutions of these ordinary 
differential equations can, for some cases, be found in closed form while others can 
only be solved approximately. 

The ordinary differential equations are very complicated and need not be solved 
to obtain the information required. It is noted that the total energy between two 
moving surfaces represented by r ' and r" respectively, is given by 


E oc 




pr v 1 dr 


(8.25) 


and that the dimensions of E are ML V ~ 1 T~ 2 . The change in fluid energy between 
the two surfaces mentioned equals the pressure work, which is proportional to 

(p v r (r ) — pv r (r ) ) 


A possible functional form of E is 

E = ab v - 1 - k t 0{v - l - k) - 2 -f 


(8.26) 


Assuming also that r'(0 and r"(t) are determined from (' = const and £" = const , 
then generally 

J JP P 

d ±=\j3(v-l-k)-2-,]j (8.27) 

From the energy integral 


J J r' 
(8.28) 

which can be made independent of /(£", C\ . . .) if 

s -/3(i/ - 1 -k) = -2 (8.29) 


[s+2-0(u-l-k)}f(C,C,---)^U 


2 


p'v' r T ( V' r ~0) 


f /2 
pv r 


7 - 1 . 


)d elding 


- i /+2 


p'v' T + (u' r -0) 



= const 


(8.30) 
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Table 8.1: Sedov’s Scaling Laws 



Unpowered 

Powered 
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P/Prtl = 
z/z„, = 

P/Pre/ = COnst. 
z ! z r ef = {t/tref) 

CYLINDRICAL^ 

P/Vrtf = {t/tref) 1 

r l r mf =(t/t„ f y /2 

P/Pr.1 
r/r„, = 

SPHERICAL 

vlvref = [t/tref ) ^ 

T l r n:f = (^/^re/ 

PlPrtf = 413 

r/r„, = ((/l„,) 3/s 


The condition of constant energy may be used to evaluate the constants in Using 
p 0 (as defined earlier) and Ej po as the fundamental constants in the exponents and 
constants in £ are evaluated and thus, 


c = 


r 


(8.31) 


A reference radius can be used with a specific reference time such that £ remains 
constant, — all of which leads directly to the equations of Table 8.1. Note that this 
outcome allows for variations in 7 within the partial differential equations. In fact, 
with the self-similar solution, all points throughout the blast wave expansion behave 
the same, so that the ! re/’ point can be arbitrarily chosen. Naturally, this would not 
occur in a flow regime with complex reflections. 

The powered scaling laws are supplements to Sedov’s original equations, and are 
useful in modeling laser energy deposition. The powered scaling laws assume an 
energy deposition that is linear in time (see Reilly et ah, [10]). Although the exponents 
are different, the problem is still self-similar. The 'powered’ laws are valid only while 
the laser is on; when the laser is off (i.e., after t p , the laser pulse duration), the 
unpowered scaling laws become valid. 

With Sedov’s scaling relations, the ordinary differential equations for n r , p and 
p‘ do not have to be solved. Since p = 2V Xack! (7 + 1) anc ^ V shock = drfdt , Raizer s 
equations give the initial values needed to calculate the impulse. An expression for 
the laser- generated impulse imparted to a thrust surface, based on a model which 
employs both Sedov’s and Raizer’s equations will now be presented. The analysis ap- 
plies Sedov’s equations to the LTD’s cylindrically symmetric ERH thruster geometry. 
Table 8.1 lists pressure relationship with time and radius of the blast wave for the 
unpowered cylindrical case. Given these relations, the impulse can be determined in 
a straightforward manner. 
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Figure 8.10: Surface Pressure vs. time using Sedov’s equations($ l$d — 500 MW/cm 2 , 
t lsd = 0.5 cm, Moo = 0 . 1 , 0 km, with Plsd at £ 2 d-) 


For example, using p LS D = Pr*/> r LSD = = 0.5 cm and t^ f = £ 2 d,( for the 

conditions of <f> = 500 MW/cm 2 , M = 0.1, and h = 0 km), Sedov scaling relationships 
give the values plotted in Figs. 8.10-8.13. Note that c # , the speed of sound, is cal- 
culated as V t hock{t)l 2, in accordance with the literature (e.g., see Reilly, [ 10 ] ) while 
V 9 hock{t) — dr(t)/dt : p — yp/c] and T = p/(pZR ). The values of 7 = 1.2, Z = 3 
were assumed in this calculation. As before, the plasma temperatures indicated in 
Figure 8.12 greatly exceed the reasonable peak values of 10,000-30,000 K. 
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Figure 8.11: Surface Pressure vs. radius using Sedov’s equations ( $lsd = 500 
MW/cm 2 , tlsd = 0.5 cm at t 2 D , = 0.1, 0 km, with Plsd at *2 d)- 

8.5 Calculating the Impulse with Sedov’s Scaling 
Laws 

To enable the derivation of the impulse, a control volume is first drawn around the 
thrust surface as sketched in Figure 8.13. This semi-cylindrical control volume is 
taken to the outer edges of the ERH thruster surface and is centered on this plate. 
On the blast wave side of the surface, the pressure is set equal to p,(t) which is a 
function of the peak LSD wave pressure (given by Eqn. 8.18) via the pressure scaling 
relationship. On the other side of the thruster surface, the pressure is assumed to be 
at ambient value. Figure 8.14 shows an end view of the control volume drawn around 
the thrust surface. Within this context, the following equation can be derived to relate 
the time-resolved blast wave pressure at the surface, p,(t), the affected thruster area 
and the resultant impulse. 

\p.(t)- P .]dA<a = 

Jo j A(0 

rto [w r r {t ) 

/ / / \p,{t) -p a }drdzdt = 

Jo Jo Jo 
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Figure 8.12: Wavefront temperature vs. time using Sedov’s equations ($lsd = 500 
MW/cm 2 , r LSD -- 0.5 cm at i 2 £>, M* — 0.1, 0 km, with T^sd at t 2 zj). 

ft o f r (t) 

W Jo Jo ^*(0 -Pa)drdt (8.32) 

Equation 8.32 is the total impulse integral which depends on the time for the blast 
wave pressure to fall to ambient (i.e., t final = t 0 = t a ), and on the radius of the blast 
wave, r(t). The blast wave radius relation, r(f), corresponds directly to the pressure 
scaling law in Table 8.1. The width, w, of the ERH thrust surface in the streamwise 
flow direction (see Fig. 8.15) is given by: 

w = 2ir(r ceol + Vagi) (8.33) 

where r ceoi and r aol are the thrust plate exit outer lip and inlet outer lip , respectively. 

The LSD wave is assumed to start at the line focus of the laser beam (i.e., virtually 
at r = 0),and then propagate radially outwards at a fairly constant pressure (Plsd) 
until the laser pulse is ended (at r^sD and t p ). In reality, a very high laser intensity 
may be used to trigger air breakdown and ignite the LSD wave; then a lower intensity 
beam would sustain a relatively constant LSD wave pressure, during which the LSD 
wave moves at constant velocity, V LSD (given by Eqn. 8.6). At t = £ p , (where t p is the 
laser pulse duration), the wave reaches a radius defined by r LSD . Hence, r LSD is a 
convenient reference radius which is then used for the blast wave decay process (after 



134 CHAPTER 8. MODELING OF LASER-GENERATED IMPULSE 



Figure 8.13: Blast wave radius vs. time using Sedov’s equations ($lsd — 500 

MW/cm 2 , r LSD = 0.5 cm at t 2D , = 0.1, 0 km). 


the laser is shut off at t p ), which is modeled with the unpowered cylindrical scaling 
law. If, for example, the LSD wave propagates to a height h off the surface, then 

t p = h/VisD (8.34) 


In this case, h = r LSD (although h may be small enough to be neglected in the 
semi-cylindrical blast wave calculation). 

Shortly after the laser pulse terminates, a rarefaction wave comes in from the edges 
of the blast wavefront of the wave, according to method of characteristics analyses 
[12], at an average plasma speed of c, lsd , gi ven h T- 


c tLSD 


Vlsd 

2 


(8.35) 


When this rarefaction fan first arrives at the centerline of the impulse surface in 
Figure 8.14, (at a time t 2D ), the blast wave geometry is assumed to evolve into a 
completely cylindrical geometry. The time t 2 q is given by: 


, !D = (8.36) 

c i LSD 

Therefore, t 2D may be taken as the reference time in the cylindrical unpowered scaling 
law. Once t 2D is known, the expansion (or impulse generation) time, £ 0 , of the blast 


i 





Figure 8.15: End view of the control volume drawn around the thrust surface. 
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wave mav be calculated as: 


t 0 — t 


Plsd 


2D' 


(8.37) 


while the fully expanded radius, r o, of the wave at to is: 


T’o = r LSD < 


(8.38) 


'2D 


Prior to the time t 2£) , impulse delivered by the ‘powered’ planar LSD wave is 
treated the same as by Pirri, (11). Thus, the first ‘powered portion of the impulse 
delivered to the square plate in Figure 8.14 is given by: 


Ipwrp - j 2D J J LS Wlsd -p a }drdzdt = 27r(r ceoi + r ciol )r LSD (p LS D ~Pa)t2D (8.39) 

Without subtracting the ambient pressure (p Q ) term, one has, 

Ipp = 2ir(r ceoi +T„ 0 i)r LSD p L sDt2D (8.40) 


The final “laser-off” contribution of the impulse, J^, from the ‘unpowered’ cylindrical 
phase, delivered to the l x l square plate, can now be calculated. 


to /r(t) 


I uc = W P LSD T LSDt 2D / / 

Jtj D J r 


LSD \Plsd Plsd / Wsd 


) d (—)d (—) 

J '■ r LSD ' \t-LD' 


= 2 1, 


Plsd 


vp 


-11 - 


Plsd 


(8.41) 


The total impulse is then: 

Pa _ Pa \ 

Plsd J Plsd J 

(8.42) 

A computer program is used to compute through the procedure outlined above. 
Further details are presented in the next chapter. Note that the maximum radial 
extent of the thrust region (over the impulse plate) is denoted by Z/2. If the Tq calcu- 
lated exceeds 1/2, then the impulse calculation is terminated. This is a conservative 
assumption since the blast wave could be convected downstream across the impulse 
plate by the inlet airflow, effectively yielding larger thrust area in the upstream di- 
rection (i.e., r- cioi ); on the other hand, r ceol would be shortened. Since the wave is 
assumed axi- sy mm etric and the inlet air velocity is much smaller than that of the 
blast wave, this effect was not considered here. However, the experimental results 
obtained by Woodroffe et al suggest that significant beneficial effects may exist with 
a supersonic inlet airstream (see Ref. 13). 


I total = 27r(r ceoi +r ciol )r LSD p L sDt2D 0 + 2 


I Plsd 

Pa 


-1 - 


i 
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The peak laser power delivered into the LSD wave across a square plate of dimen- 
sions l x l is: 

P = 2irTLSD( r ceol + r dol)$LSD (8.43) 

From this, the peak laser energy can be obtained as follows: 

E v = Pt p = 2-Kr L s D {r ceol + t 'doi^LSDtp (8.44) 

The time required for fresh, unprocessed air to convect into the impulse region (in 
preparation for the next blast wave) is defined as t re f n , h , given by: 

, */ 2 r ° ,o 

trefre.h = Or — (8-45) 

initial initial 

whichever is smaller, because this fresh air is assumed to enter after the blast wave 
expansion process is complete. Thruster performance is optimized with the smaller 
refresh time. A total engine cycle time can then be defined as: 

t cycle — ^0 (8.46) 

The engine pulse repetition frequency, PRF , is: 

PRF = — = — - (8.47) 

c cycle ^0 f Lrefrcth 

Next the engine duty factor, DF, which can be interpreted as that portion of the total 
cycle time, over which laser energy can be applied, is: 

DF = t p PRF (8.48) 

Then, the time-averaged laser power spent over the cycle is: 

P = Pave = P X DF (8.49) 


Finally, the time-averaged thrust delivered to the l x l square plate is: 


T — T ave — I total P EF COS 8 C — D total 


(8.50) 


where 6 C is the orientation angle for the shroud lower impulse surface of the LTD. An- 
other important parameter concerns the efficiency of impulse delivery to the thruster 
surface: the Coupling Coefficient, CC , or efficiency of thrust production from the 
available input laser power, defined as: 



(8.51) 



138 


CHAPTER 8. MODELING OF LASER-GENERATED IMPULSE 



Figure 8.16: Cross sectional view of LTD ERH thruster surface 

As shown in Figure 8.16, the LTD primary optic focuses the laser beam to an 
angle of 85° upon the ERH thruster surface, and into a LSD wave cross-section that 
is roughly se mi -circular. Hence, it is necessary to scale the actual LSD-wave plasma 
geometry and input laser pulse energy to the one-dimensional assumption treated 
above (i.e., E v in Eqn. 8.44). Thus, only ((180/360 )tt)/ 2 = 78.54% of the 1-D model 
laser pulse energy is really absorbed into the gas; this sc alin g ratio is incorporated in 
the computer code. 

The principal results of the comprehensive engine performance analyses for the 
Lightcraft Technology Demonstrator, are presented in the next chapter after the 
LTD aerodynamic inlet model is introduced. Finally, note that the important param- 
eters described above are all plotted for the sample case of r LS D — 0.5cm, and $lsd — 
bOOMW/cm 2 . 
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Flight Mach number 


Figure 8.19: The time required for blast wave expansion to ambient pressure (<F 
= 500 MW/cm 2 , t LS d = 0.5 cm) 


10 km 


/ 20 km 


Flight Mach number 


Figure 8.20: Gross impulse produced by the laser- generated blast wave for one com 
plete engine cycle (with units of N — sec), for various vehicle flight conditions ( $lsi 
= 500 MW/cm 2 , r LS D — 0.5 cm) 
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Chapter 9 


THE LTD ENGINE PERFORMANCE 
ANALYSIS 


9.1 Analysis of LTD Inlet Flow Aerodynamics 

A simple LTD inlet flow model was developed in order to determine the state of 
the inlet air (pressure, temperature, velocity, etc.) which refreshes the lower annular 
shroud impulse surface. This aerodynamic model model ‘was needed not only for 
the refresh air state, but also for defining external drag characteristics of the entire 
vehicle. In the latter, it is important to distinguish between the drag produced by 
engine-related vs. airframe-related components, so as not to double penalize overall 
vehicle performance. 

9.1.1 Investigative Procedure 

Several basic assumptions in modeling the flow over the LTD vehicle forebody were 
made to simplify the aerodynamic inlet analysis. The primary objective of this model 
was to adequately specify fresh air flow conditions over the ERH thruster impulse 
area. The flow around the vehicle was assumed to be inviscid. It is also dominated 
by compressible flow effects since the flight regime of interest can extend to Mach 8. 

A number of important considerations enter in the design of the vehicle forebody, 
shroud and inlet configuration. First, the inlet must provide a reasonable compnmise 
between subsonic, supersonic and hypersonic performance. Conventional wisdom 
for efficient hypersonic inlets dictates a slender conical spike with a sharp point; 
potential pitfalls for this approach include an excessive inlet weight penalty, large 
boundary layer buildup and extreme heat transfer rates at the point. For certain 
hypersonic inlets, the boundary layer can grow thick enough to fill the entire inlet 
gap; temperatures within this boundary layer can be very high. 

The LTD geometry in Figure 9.1 is designed to automatically bleed this hot hy- 
personic boundary layer; also, the shroud inlet gap is placed at a sufficiently large 
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radius to prevent the bow shock from being swallowed below the design Mach number 
(e.g., Mach 5). In addition, the inlet is fairly short so that the boundary layer does 
not grow large at hypersonic speeds. Although the total pressure recovery will be 
somewhat lower than a slender cone inlet could provide, aerodynamic heating will be 
much reduced for the blunt geometry, and passive thermal protection systems should 
be adequate; in contrast, a pointed cone would require active cooling of the tip at hy- 
personic speeds. Finally, the blunt inlet should give superior performance in subsomc 
flight. Clearly, the best compnmise for the above mentioned inlet requirements can 
only be derived from a thorough experimental testing program. The LTD’s engine is 
designed for the acceleration role and must exhibit a good measure of performance 
over a wide Mach number range; it attempts this goal with a fixed inlet geometry, 
but a minimally variable geometry (e.g., a translating shroud, or variable “kicker” 
plates) could easily be included. 

Perhaps the most important inlet parameter is the velocity of air convected across 
the thrust generation site. To determine the flow velocity at this location, the inlet 
gap velocity is determined by the aggregate flow conditions over each of the cone- 
shaped segments that make up the LTD inlet forebody. For example, in subsomc 
flight the conical LTD forebody will cause the external flow to accelerate into the 
inlet gap at a velocity greater than the freestream, which is typical of cones and/or 
spheres (e.g., see White [1])* Here, the ratio of inlet gap to freestream velocity was 
set equal to 5/3 — up until some inlet conditions are reached; thereafter, choked inlet 
flow was assumed to remain until a flight Mach number of 1. For subsonic flight, this 
inlet Mach number schedule and a fixed inlet gap area directly determines inlet mass 
flow rate and ram drag, required in thrust calculations. 

In supersonic flight, a strong bow shock is formed over the LTD inlet forebody. 
The Mach number after this oblique conical bow shock was calculated using the 
conical shock relations for a semi- vertex angle of 30° (see Shapiro [2], Kuethe and 
Chow [3], Oates [4] and NACA 1135 [5]). These relations were also used to compute 
the stagnation pressure ratio, (pti/pt«)> static pressure ratio (pi/p a ), a ^d bow shock 
angle, 0 X . 

Once the gas conditions behind the conical bow shock are known, a similar process 
is applied to obtain the flow state behind the shock which forms over the truncated 
second cone which has a semi- vertex angle of 51°. This flow transition point occurs 
sufficiently far back on the first cone, at which the oblique shock relations for a 21° 
wedge can be applied — to obtain M 2 , Vt*IVt\ an( ^ P 2 /PI) an d 02 - 

Figure 9.3 shows the results of this analysis for M 1} M 2 and M/;v (i.e., at the inlet 
gap). Note that M [N is assumed to be some beyond a flight Mach number of 0.6. 

The overall pressure recovery of the LTD inlet is the ratio of inlet gap stagnation 
pressure to the freestream (p t JJV /p too ), which is obtained by multiplying the stagnation 
pressure ratios across the first and second shocks. Note in Figure 9.2 that the inlet gap 
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Figure 9.1: LTD forebody flow parameters 
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0 12345678 


Flight Mach Number 

Figure 9.3: LTD forebody /inlet Mach number schedule assumed for supersonic flight 
speeds. 30° and 51° cones (semi-vertex angle) model the forebody geometry. 


is choked beyond a flight Mach number of 0.6. The inlet pressure recovery schedule 
is plotted in Figure 9.4. No real gas effects have been included in this analysis. 

With the above information, the inlet gap velocity can now be calculated as: 

(9.1) 

where p in is the inlet gap static pressure and pj n , the density. The inlet static pressure 
and temperature are plotted in Figures 9.5 and 9.6 (respectively) as a function of LTD 
flight Mach number and altitude. 

First the stagnation temperature is calculated from 



— = 1 + -Ml 

2 


The stagnation temperature is constant across all shocks throughout the inlet flow; 
hence, (T too = T tin ). This fact is used with the following relationship to calculate T in . 


= 1 + 



refresh Ptniet (^*0 



It oO 
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Figure 9.6: Inlet static temperature vs. Mach number and altitude. 


In the free stream region, 


Ptoo/Poo = (TtoJT^r'^ 


With the assumed inlet pressure recovery schedule, p t ^ n is then computed as 


Ptin — 



PtOQ 


The air density at the inlet is simply 



In the analysis of overall LTD flight performance, the LTD centerbody is con- 
sidered the “vehicle” component, and the annular shroud is treated as the “engine.” 
The total drag of the LTD vehicle includes viscous (or frictional) drag, pressure drag, 
spillage drag and inlet ram drag. The viscous drag is neglected because the effect 
is small. Pressure drag results from the differences in static pressure applied across 
the fore and aft surfaces of the vehicle centerbody. Spillage drag occurs when all 
the inlet air flow passing through the frontal “capture” area (i.e., within the circular 
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Figure 9.7: Refresh time vs. flight Mach number. 

stream-tube cross-section defined the shroud leading edge, r^) cannot be swallowed 
by the inlet gap and must be diverted out over the shroud upper surface. 

One of the most important engine-related drag components modeled here is the 
inlet ram drag which is defined as: 


Dram — m in u in 


(9.2) 


where 

^in PinV'xn^-iniet P xn^iri^ {T ciol COS S cone 2 (0-3) 

assuming that M refn , h = M inlet , u inlet = u nfreih and that the distance from the 
inlet gap to the blast wave ignition point is 1/2, the refresh time is simply given by 


1/2 T*o 

tre/re*h — Of 

u refresh ^ re freak 


(9.4) 


whichever is smaller, as discussed earlier. The refresh time, based on this expression 
is plotted in Fig. 9.7 as a function of LTD flight Mach number and altitude. 

Plotted in Figure 9.8 is the total ram drag coefficient, where c D — 2 D ra M / P oouloA totai 
note Atotoi = is the maximum frontal “capture” area of the inlet. 

An important parameter that checks accuracy of the aforementioned inlet assump- 
tions is the ratio of the inlet gap mass flow rate at the choked flow condition, M/jv 
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Flight Mach number 


Figure 9.8: Ram drag coefficient for the LTD ERH thruster. 
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Flight Mach number 

Figure 9.10: LTD inlet mass flow rate vs. Mach number and altitude. 


= 1. to the maximum “capture” mass flow rate defined by the inlet lip (i.e., r CTO /); 
the latter is given by = p^irr^u^. This ratio, which physically cannot exceed 
1, is plotted in Figure 9.11. Note that present LTD inlet assumptions indicate a 
design Mach number of about 5. The actual LTD inlet characteristics are best left 
for experiments to confirm. 


Another critical ERH thruster parameter is the resultant pressure at the inlet 
gap, produced by a decaying blast wave after it has expanded across the shroud lower 
impulse surface. Cylindrical blast wave characteristics are calculated (as described 
earlier) with Sedov and Raizer’s equations, and plotted in Figure 9.12. Note that 
blast wave pressure does not decay fast enough to avoid exceeding the inlet static 
pressure under certain flight conditions. Since fluids will generally tend to move from 
regions of high to low pressure, the inlet flow might reverse direction and separate 
from the forward lip. This might limit the engine operational flight regime to the 
higher dynamic pressure trajectories where the parameter in Figure 9.12 is always 
less than the stagnation pressure at the inlet gap. Again, experimented confirmation 
of this assumption is necessary. 
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* 



Figure 9.11: Ratio of inlet gap to maximum “capture” mass flow rate, vs. flight Mach 
number. 
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Figure 9.12: Blast wave pressure at the inlet outer lip (i.e., r = 10.8 cm), in bars, 
vs. flight Mach number and altitude, with $lsd = 500 MW/cm 2 using Raizer and 
Sedovas equations. 
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Figure 9.13: Model of subsonic flow over a cone. 

9.1.2 Alternate Subsonic Inlet Model 

This alternate procedure attempts to justify the above assumptions made in modeling 
the subsonic inlet flow to the LTD engine. The main objective is to obtain forebody 
pressure distributions for inlet and vehicle drag calculations. 

The LTD forebody flow model for the subsonic case is based on the known aerody- 
namic flow field that occurs over a sphere; the inlet geometry is shown in Figure 9.13. 
The fluid velocity along a spherical surface can be obtained from invisdd potential 
flow theory (see Schlichting): 


u/uoa = 2 sin {x/R) 

where the variables are noted in Figure 9.13. A factor of 3/2, instead of 2, is suggested 
for the cones. This gives a critical flight Mach number of 0.6 for the choked inlet 
condition (i.e. , Mj„; et = 1), in accordance with previous discussions. Note that the 
subscripts oo, a and 0 all represent the same ambient location far upstream of the 
vehicle inlet. 

From Figure 9.13, it is observed that x = \A* 2 + Z 2 and R = —r* 2 , where 

l cma x i s defined as Z c = r c . Note also that x ~ rj sind), and R % r/ sin<p. From 
these considerations, one can express the velocity, u cl , at any point along the front 
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forebody cone as 


^cl 3 . | Id" ( ^ cone 1 /^* cone l) 

Uoo 2 \J 1 + (1 — lconel/ T conel) 


For the second cone, 


Uc2 

«el 


cone 2 ^*conel)] 


3 . 1 + [Zcone2/ ( r 



2 \( 1 (1 - Icot.j/Ii'otk! -rame 1)] J 


(9.5) 


(9.6) 


By introducing the definition for the Mach number at any point along the cone, 
M = u/c, where c, = yj'iopj p, one obtains 


,, ;7oPo 3 . 
u = Moo. x sm (- 

V Po L 


or 

. M = M sin(---) 

V V Po 2 

Also, from isentropic, steady, compressible flow equations, 
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which is simply: 


M = Moo 


! 1 + ( 7o -1)M 2 /2 3 . 
\| 1 ( 7o - l)A/i/2 2 S ‘ 


After solving for M, 


M = 


l+iTj-ljM*, p. *oo 



(9.10) 


(9.11) 


By replacing u/u* with either u cl /u „ or M 2 (right in front of cone 2), or 

Mi n (right in front of the inlet gap) can be calculated. Equation 9.11 may be used 
for any part of the flow that is subsonic. 

Note that the results presented in Figure 9.14 basically serve to confirm the earlier 
subsonic inlet analysis, and may aid in the study of other forebody cone angles. 
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Figure 9.14: Alternate forebody /inlet Mach number schedule; 30° and 51° cones 
model the forebody. 

9.2 LTD Performance Projections 

The results in this section present a complete analytical performance projection for 
the Lightcraft Technology Demonstrator airbreathing engine. The impulse coupling 
coefficient, measured in units of Newtons of thrust per megawatt of input laser power, 
remains the principal measure of engine performance. 

A base line case which assumes an inlet gap height of 3 cm, is the central focus of 
the analysis. Several engine- related parameters were varied in search of the optimum 
coupling coefficient regime and to understand the performance sensitivity to those 
parameters. The principal baseline case parameters consist of the following: tlsd = 
0.5 cm and $lsd = 500 MW/cm 2 and = 1.0 when Moo > 0.6. As described 

earlier, Raizer and Sedov equations are used in the impulse model, and a simple fore- 
body/inlet model determines the air refresh conditions in the thruster. The baseline 
net coupling coefficient performance is displayed in Figure 9.15. 

Also for the baseline case, Figure 9.16 gives the laser pulse duration; Figure 9.17, 
time-averaged input laser power; Figure 9.18, time-averaged thrust; Figure 9.19, pulse 
repetition frequency; and Figure 9.20, pulse energy. In figures 9.21 and 9.22, the 
gross time- averaged thrust and gross coupling coefficient performance are displayed, 
respectively. Note that gross time-averaged thrust minus inlet ram drag (Fig. 9.9) 
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Figure 9.15: Net coupling coefficient vs. Mach number and altitude, for the baseline 
case (inlet gap = 3 cm, $£,s£> = 500 MW/cm 2 , r^sD = 0.5 cm, and = 1.0 

when Moo > 0.6). 
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Flight Mach number 

Figure 9.16: Laser pulse time vs. Mach number and altitude, for the baseline case. 



Figure 9.17: Time-averaged power vs. Mach number and altitude, for the baseline 
case. 
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Flight Mach number 

Figure 9.18: Net time-averaged thrust vs. Mach number and altitude, for the baseline 
case. 


3 km 



Flight Mach number 


Figure 9.19: Pulse Repetition Frequency vs. Mach number and altitude, for the 
baseline case. 
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Flight Mach number 


Figure 9.20: Pulse energy vs. Mach number and altitude, for the baseline case. 

equals net time-averaged thrust (Fig. 9.18); also, time-averaged net thrust divided by 
time-average input power (Fig. 9.17) equals net coupling coefficient (Fig. 9.15). 

Finally, Figure 9.23 displays an important parameter that gives some insight into 
the dynamics of the refresh process for the baseline case: t refreth /t expan , lon , where 

^ expansion — ^ 0 - 

Note from Figure 9.23 that it may be possible for the LTD to fly along a boost 
trajectory where this parameter is close to 1. At this condition, the time required for 
fresh air to refresh the ignition point ) is exactly equal to the time necessary 

for the blast wave pressure to expand down to the local ambient level. This criteria 
has been employed to size the streamwise dimension (i.e., “Z”) of the annular shroud 
impulse surface; hence the heated air-plasma annular “bubble” would relax to ambient 
pressure, just as the aft bubble edge meets the trailing edge of the shroud. At this 
point, the next laser pulse would be delivered, because fresh air has simultaneously 
reached the ignition site. It bears mentioning that the plasma bubble is assumed 
to convect aft, immediately following ignition; this contradicts limited experimental 
evidence to the contrary, which suggests the existance of an initial delay time — 
until the plasma bubble pressure drops alm ost to the local stagnation pressure level 
(i.e., the sum of both dynamic and static pressure components) — when the ‘bubble’ 
begins to move aft. All of this has direct implications on how large the annular 
shroud impulse area must be; however, the resolution of this issue must wait for an 
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Figure 9.21: Gross time-averaged thrust vs. Mach number and altitude, for the 
baseline case. 



Figure 9.22: Gross coupling coefficient vs. Mach number and altitude, for the baseline 
case. 
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experiment. 

At the present state of knowledge, these combined air refresh flow and laser- 
heating effects are difficult to model analytically; hence, four variations in engine input 
assumptions (i.e., M^y,^ , inlet gap height, t LS d and $lsd) were run to determine 
their separate impact upon coupling coefficient performance. 

First to be studied was the effect of refresh Mach number (M^y^,/, = 0.75, 0.5 and 
0.25) at flight Mach numbers beyond 0.6.; the flow velocity is assumed fixed from the 
inlet gap to the ignition site. All other parameters were held constant at the baseline 
values, including the aforementioned inlet flow velocity schedule from 0 < M <0.6 
(wherein, inlet Mach number increases to 1.0 at flight Mach 0.6). In comparing the 
coupling coefficient ( CC ) performance for M^y,*,* = 0.75 with the baseline case, CC 
values are found to be significantly improved between 0.6< < 4, with no decrease 

in top end performance Also note that CC goes to zero at about M,* = 5.6 and 30 km 
altitude, just like the baseline case; these cases agree well with the inlet design Mach 
number 5.2 suggested in Figure 9.11. As M^y,^ is further reduced 0.5 in Figure B.2, 
low end CC performance is seen to increase somewhat beyond the = 0.75 

case, but the top end performance begins to suffer (with a zero CC at Moo = 5.5 and 
30 km). This situation becomes unacceptably exaggerated for the M^y^,^ = 0.25 
case in Figure B.3. 

From this inlet Mach number sensitivity analysis, one concludes that the ERH 
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thruster will accept a wide range of M refreth values from 0.5 to 1.0; the lower refresh 
Mach numbers seem to indicate improved CC performance, but might degrade time- 
average thrust levels and PRF (this impact was not examined, however). 

Examined next was the sensitivity to inlet gap height on CC performance. Fig- 
ure B.4 illustrates the effect of reducing the gap height by 1 cm below the baseline 
value of 3 cm. Note that the CC performance increases significantly at the top end, 
where the CC at 30 km altitude goes to zero at M«, = 7.4; this may in part be due to 
a large reduction in ram inlet drag — associated with a reduced inlet mass flow rate. 
However, this inviscid analysis does not take into account the forebody boundary 
layer thickness, which might grow to a significant fraction of the inlet height, and vis- 
cous effects will reduce the air velocity there; note also that the top end performance 
greatly exceeds the Mp of 5.0 assumed for the inlet. If, on the other hand, the gap is 
increased by 1 cm to 4 cm total, the net coupling coefficents (see Fig. B.5) suffer at 
the top end, where they fall to zero at M TO = 4.75 at 30 km altitude. 

The next variable to be vaned was r LSD which had a baseline value of 0.5 cm. 
The additional values of 0.1, 0.25, 0.75 and 1 cm were investigated. As shown in 
Figures B.6 through B.9, shortening r LS D decreases the high flight Mach number 
(‘high end’) performance, but increases the low Mach number (‘low end’) performance. 
This is probably because the blast wave impulse dies off faster with shorter laser 
pulse times and smaller initial plasma volumes. Also, the laser energy available for 
the thrust cycle is less with shorter pulse times, which in turn is less with lower tlsd- 
For the lower flight Mach number range, however, the faster cycle times (which result 
from shorter laser pulses) allows more laser pulses to be delivered per second. Thus 
the effect of a higher PRF dominates the low end performance where the air density 
is higher (i.e., lower altitudes) with the lower Mach numbers. A r LSD of 0.5 cm seems 
to yield a good compnmise between high end, and low end engine performance for 
the baseline $lsd °f 500 MW/cm 2 . 

The effect of different laser pulse intensities can also be very important. As men- 
tioned earlier, 5 x 10® W/ cm 2 was chosen as the baseline case; 3 x 10® and 10 10 W/ cm 2 
were also run. While breakdown occurs at roughly 10® W / cm 2 for 10.6 /i m wavelength, 
the breakdown threshold is a direct function of wavelength: $ Breakdown oc A -3 / 2 . By 
using lasers in the shorter infrared (e.g., 1 fi m) and visible (e.g., 0.7 fim to 0.3 /im) 
ranges, $lsd can be greatly increased beyond the baseline value of 3 x 10 8 W/cm 2 . 
However, higher $lsd’ s mean even greater plasma temperatures, which push the 
credibility limits of this straightforward ideal gas model. Hence the results of this 
sensitivity analysis on $lsd which are presented in Figures B.10 through B.27, should 
be viewed in a qualitative sense only. 

The principal effect of increasing laser intensity ($£,sd) appears to be a large 
associated reduction in laser pulse times; also, as risD is varied at a new higher $lsd-> 
the optimum CC performance is shifted to a shorter t^sd- Interestingly enough, there 
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is a value of r LSD for which the original baseline case of r LSD =0.5 cm and 4 > LSD = 
500 MW/cm 2 seems to be almost duplicated: for example at $lsd = 3 x 10 9 W/cm 2 
it is tlsd — 0.25 cm (see Fig. B.ll); at $lsd = 10 10 W/cm 2 it is t^sd — 0.2 cm 
(see Fig. B.20). Whereas the baseline laser pulse duration is roughly 0.3 ps along a 
typical launch trajectory (see Fig. 9.16), at 3 x 10 9 W/cm 2 it shortens to 70 ns; with 
10® W/cm 2 , t p falls to 40 ns. 

In general, the time-averaged input laser power (see Figs. 4.17, B.18 and B.24) 
and net thrust (see Figs. 4.18, B.14 and B.25) for the three optimal cases mentioned 
above are nearly identical. The pulse energies (Figs. 9.20, B.16, and B.27), and the 
pulse repetition frequencies (Figs. 4.19, B.15, and B.26) fall in similar ranges also. 

In summary, while much is known about laser-induced breakdown and impulse- 
generation within a static atmosphere, minimal experimental data exists on laser 
impulse coupling within a gasdynamic environment. Although this “refresh” aspect 
of the ERH thruster has a certain similarity to all airbreathing engines, a clear funda- 
mental understanding must await critical proof-of-concept experiments. This chapter 
reveals analytical performance projections for the ERH thruster cycle that (taken as 
a whole) look exceedingly promising. The key parameters are clearly layed out such 
that a well conceived experimental program can now be planned. 
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Chapter 10 


LASER-HEATED ROCKET MODE 


As presently conceived, the LTD vehicle could be constructed using state-of-the-art 
components derived from current liquid chemical rocket engine technology, advanced 
composite structures, and high power laser optics developed for SDIO applications. 
In the rocket mode, the regeneratively cooled, 1 m diameter primary optic will also 
serve as a plug nozzle; hence, the parallel with liquid chemical rocket technology 
becomes obvious. 

The LTD combined-cycle engine would start off in an airbreathing mode, then 
“shift” into the rocket mode at Mach 5 and 100 kft. Fourteen kilograms of LN 2 
would be expended in the airbreathing mode (to cool the primary mirror), leaving 
126 kg for the rocket mode. The rocket specific impulse will be > 875 s (coupling 
coefficient, CC ~ 100 N/mW), using LN 2 as the reaction propellant. This single- 
stage-to-orbit vehicle has a mass ration of 2.0, and the entire final mass of 124 kg 
(i.e., 120 kg dry mass, plus 2 kg of LN 2 and 2 kg of He ullage) becomes the payload. 

Soon after arriving in orbit, the remaining 2 kg of LN 2 evaporates, helping to cool 
down the vehicle, while maintaining adequate propellant pressure. This compressed, 
cold N 2 &: He mixture is then utilized by a simple 3-axis rocket attitude control 
system to accomplish fine pointing required for sensor satellite functions throughout 
its life-time. 

A cutaway view in Fig. 10.1 shows the internal arrangement of hardware com- 
ponents for the LTD machine; the projected mass breakdown is given in Chapter 2. 
Note that the two most massive components are the 28 in. diameter LN 2 tank and 
the annular shroud engine. The LN 2 tank is a filament wound pressure vessel similar 
to the ones made by Brunswick Defense in Lincoln, Nebraska. 

Figure 10.2 displays the range of primary optic (PO) contours considered for the 
LTD machine. These contours are generated by rotating a parabola about an off- 
center axis. Note that the annular shroud airfoil has a circular arc cross-section with 
a fiat bottom, and a chord (C) length of 22.5 cm. Five PO contours are displayed with 
focal rings positioned at 16.67 %C , 22.0 %C, 33.33 %C, 41.67%C, and 50.0%C. The 
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Figure 10.1: Cutaway View of LTD Aerospace Vehicle 

longest PO afterbody has its focus at the 16.67%C chord point. As mentioned earlier, 
one wants to place the focus as far forward on the shroud as possible, yet at the same 
time minimize the PO weight penalty. For the LTD machine, this compromise was 
reached by placing the focus at 33.33%C. Note that the PO contours in Fig. 10.2 
bear a striking resemblance to isentropic spike rocket engine nozzles. 


10.1 Plug Nozzle Rocket Engines 

Graham and Bergman were the first to build plug nozzle rocket engines back in the 
late 50’s and early 60’s while working for General Electric. Figure 10.3 shows their 
16,000 lb (71.2 kN) thrust, H 2 0 2 uncooled plug nozzle thrust chamber mounted in 
a test cell. (Incidentally, this engine developed a thrust level roughly equivalent 
to the LTD engine’s maximum thrust in the airbreathing mode.) As indicated in 
Fig. 10.4, six H 2 0 2 decomposers were arranged around the base of the plug to provide 
exhaust gases for hot tests. Two, of the many external expansion nozzles which were 
tested, are displayed in Figs. 10. 4A and B; an isentropic plug, and a 20° foreshortened 
plug. This rocket engine demonstrated good performance at off-design pressure ratios, 
and proved the feasibility of thrust vector control by selective combustor throttling. 
Subsequent analysis of the test data revealed that full isentropic plug nozzles could 
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Figure 10.2: Primary Optic contours vs. Focal Locations. 


be subjected to substantial truncation/foreshortening, without degrading the thrust 
performance to any large degree. Figure 10.5 indicates the method used to foreshorten 
the plug nozzles, by replacing the tip with a cone at the appropriate tangent location. 

For the first flight-weight demonstration of a regeneratively cooled plug nozzle 
rocket engine, Graham and Bergman decided to completely truncate the isentropic 
spike nozzle, replacing it altogether with the 42° half-angle cone shown in Fig. 10.6. 
Calculations were carried out to determine the minimum size of plug nozzle that could 
cool itself with RP-l/LOX propellants. The result was the 50,000 lb thrust stainless 
steel engine displayed in Fig. 10.7. As indicated in Fig. 10.8, the plug nozzle had an 
exit nozzle diameter of 1.07 m (which is identical to that of the primary optic of the 
LTD machine). The overall maximum diameter of the 50,000 lb thrust engine was 
1.34 m. 

Figure 10.9 shows one of the eight cooled segments of which the engine is com- 
posed; each segment could be independently throttled to accomplish thrust vectoring. 
Pictured in Fig. 10.10 is an enlarged view of the lower cone segment, revealing RP-1 
coolant flow passageways. 

The lower cone segments were fabricated from two thin stainless steel faceplates 
separated by numerous vertical ribs that formed the coolant passageways, all fused 
together by an electron-beam welding machine. Each of the eight lower cone segments 
had a mass of 2.077 kg and an area of 1089 cm 2 ; this translates to a mass penalty of 
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Figure 10.3: 16 Klb H 2 0 2 Plug Nozzle Development Thrust Chamber in Test Cell 
(After Graham and Bergman, Ref. 1) 
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Figure 10.4: Detail of 16 Klb H 2 O 2 Plug Nozzle Configuration (from Ref. 1) 



170 


CHAPTER 10. LASER-HEATED ROCKET MODE 



Figure 10.5: Isentropic Plug Nozzle, External Expansion (from Ref. 1) 

1.9 grams for every 1 cm 2 of regeneratively cooled surface area. 

It is interesting to imagine how this method of construction might be applied to 
the design of a cooled primary optic for the LTD spacecraft. Similar techniques are no 
doubt being investigated for large cooled optics'in SDI power-beaming applications. 
The shortest primary optic contour in Fig. 10.2 (i.e., focus at 50%C) has a total 
surface area of 17,370 cm 2 ; at 1.9 g/cm 2 , a stainless steel mirror might come in at 
33 kg. Using a thin lightweight metal faceplate with a carbon fiber composite sand- 
wich substructure, perhaps the 33.3%C optic could be built for 22 kg. A schematic 
diagram showing the internal details of an isentropic plug nozzle rocket engine with 
a geometry that is similar in shape to this primary optic design is given in Fig. 10.11. 


10.2 Injector Head Design 

Pictured in Fig. 10.12 is one of the eight ring injectors that forms a complete annulus 
at the top of the 50 Klb thrust rocket engine designed by Graham and Bergman. 
One might envision a similar injector head, reduced in width to 0.5 cm, for the LTD 
spacecraft at the mirroris annular focus region. Laser supported detonation (LSD) 
waves must be ignited at this location with absolute reliability. Therefore, these 
injector heads must be united with a tuned ignition array [2,3] or other mechanism 
to minimize the laser energy (and time) required to ignite LSD waves. Other “hot 
section” surface areas within the LTD engine which must be actively cooled (see 
Fig. 10.1) include the shroud support struts, and entire lower surface of the shroud 
— which is the principal impulse coupling region. The secondary impulse coupling 
surface is, of course, the primary receptive mirror, which serves to turn the rocket 
exhaust flow back towards the axial direction. 
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Figure 10.7: Cooled Engine Prior to Wire Wrapping and Thrust Mount Installation 
(from Ref. 1) 



Figure 10.8: 50 Klb Plug Nozzle Configuration (Dimensions in Inches)(from Ref. 1) 


i 







174 


CHAPTER 10. LASER-HEATED ROCKET MODE 


ENGINE mounting RING 


OXIDIZER MANIFOLD 


w me WRAPPED 
TUBULAR WALL 


retainer ■ 


foam plastic 




RESTRAINING 
BA NO 


FUEL COOLANT 
MANIFOLD 


FUEL COOLANT 
TUBE MANIFOLD 


‘ TUBULAR WALL 


. RESTRAINING BAND 
PLUG NOZZLE 


Figure 10.11: Quarter Sector of 33.33%C Primary Optic (from Ref. 1) 
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Figure 10.12: Injector Model III (from Ref. 1) 
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Chapter 11 


EARTH-TO-ORBIT TRAJECTORY 
SIMULATION 


The trajectory analysis of any launch vehicle becomes a critical step in the overall 
system integration process. Many important engine/vehicle related characteristics 
of the craft must come together for the final product: performance. Launch vehicle 
performance is typically measured in terms of payload capability which for the LTD 
is the entire dry mass and ullage gas, totaling 124 kilograms. This vehicle must be 
capable of attaining low Earth orbit with the available amount of propellant while 
minimizing the total laser energy consumed along the insertion trajectory. 

The trajectory was evaluated using a computer tool called SORT (Simulation 
and Optimization of Rocket Trajectories) which was written by McDonnell Douglas 
Astronautics Co. for NASA to design space shuttle trajectories [1], The code is 
sufficiently general such that it can analyze any trajectory (even those for other 
planets), and model all important environmental parameters that affect the vehicle 
dynamics. Engine performance, vehicle aerodynamics, guidance algorithms, and mass 
histories interact with atmosphere and gravity models. 


The capability of the SORT program is quite extensive in its use of these sophis- 
ticated vehicle and environmental models. The program can iterate on trajectory 
parameters to optimize performance, achieve a desired criteria, or constrain the so- 
lution to avoid some specified limit. These capabilities make SORT one of the most 
powerful trajectory analysis tools available in the world. 

Even with all of the generality built into SORT certain modifications were required 
for the Lightcraft Technology Demonstrator vehicle. The most significant software 
modification involved the unique energy source, a laser. A new vehicle steering option 
was encoded so that the Lightcraft could always point at the ground-based laser power 
station or laser satellite which initially passes overhead (see Fig. 11.1). Since the use 
of laser propulsion was not anticipated by the authors of SORT, software changes 
were also needed in the engine model in order to evaluate the total amount of laser 
energy expended along the trajectory. After these modifications to the code were 


177 


178 


CHAPTER 11. EARTH-TO-ORBIT TRAJECTORY SIMULATION 


included, SORT was capable of modeling Lightcraft performance to a high degree of 
accuracy. 

11.1 Vehicle and Environmental Models 

11.1.1 Aerodynamics 

For this preliminary performance analysis, vehicle drag was the only aerodynamic 
force (or moment) that was modeled. Since the Lightcraft is axisymmetric, it produces 
no lift at zero angle-of-attack. During flight, the vehicle’s angle-of-attack stays fairly 
small due to the restrictive steering requirements and high thrust to weight ratios. It is 
suspected that aerodynamic lift can possibly be used to benefit Lightcraft performance 
along the trajectory. 

The vehicle drag profile is shown as the dashed line in Fig. 11.2 (from Ref. 2). The 
conservatively high drag coefficients indicated for the LTD are the result of the blunt 
vehicle configuration. In the analysis of hypersonic airbreathing launch vehicles, it 
is not always immediately evident as to what should be classified as “engine” versus 
“airframe” — especially when large portions of the vehicle forebody provide the hy- 
personic inlet function. Hence careful inventory must be kept to avoid penalizing the 
engine for what might be more appropriately accounted against the vehicle airframe. 

In the LTD, the annular shroud is classified as the “engine,” and the double-coned 
center body is treated as the “airframe.” This is a conservative approach that results 
in very high vehicle drag coefficients. As portrayed in Fig. 11.3, the LTD conical 
forebody serves to turn the external inlet airflow so that it “blows” in a direction 
parallel to the annular shroud. Here it is assumed that the unstalled annular shroud 
airfoil produces negligible lift and drag. In sharp contrast, flow over the conical 
centerbody is assumed to separate at the sharp outer edge (see Fig. 11.3), so that 
the entire afterbody primary optic is immersed in a turbulent wake. Hence, the 
theoretical and experimental drag coefficient data for a 30°, conical head (half-angle) 
presented in Fig. 11.2 would be appropriate. Note that the subsonic drag coefficient 
was taken to be 0.45 which closely agrees with the data in Fig. 11.4, also from Ref. 
2, for a 30° half-angle cone. 

11.1.2 Engine 

Detailed discussions for the two engine modes are presented in earlier chapters of this 
report. The flight performance for the airbreathing mode was calculated for a wide 
range of altitudes, Mach numbers, and power settings. These values were then loaded 
into trivariant tables accessible to the program. The SORT program interpolated this 
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Figure 11.1: Ground-Based Laser Launch 
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Figure 11.2: Theoretical and experimental drag coefficients of various conical heads 
at transonic and supersonic Mach numbers (from Ref. 2). 
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Figure 11.3: Classification of “engine” and “airframe” components, for the purpose 
of determining vehicle drag profile. 
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Figure 11.4: Drag coefficients of wedges, cones and similar shapes as a function of 
their half-vertex angle (from Ref. 2). 
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data to find the thrust, specific impulse and power consumption at any point along 
the trajectory. The engine model also calculates the cumulative laser energy used 
throughout the launch trajectory, which also can include atmospheric attenuation of 
the laser beam. 

11.1.3 Steering 

The Lightcraft Technology Demonstrator must point directly at the ground-based 
laser (GBL) or the relay satellite at all times during powered flight. Two laser power 
transmission geometries were assumed for the four trajectory simulation cases which 
have been carried out on the LTD. The geometry used in Case I is portrayed in 
Fig. 11.5. The stationary ground-based laser is assumed to have an elevation of 3 km 
(density altitude). The LTD is accelerated to an initial velocity of 100 m/sec by 
a compressed air cannon, which ejects the vehicle at an initial launch angle of 30°. 
One-half second later the GBL is switched on, and the “shroud-lift” ERH thruster 
begins the laser powered ascent. 

As shown in Fig. 11.6, a boost reference angle of 30° is maintained until 50 seconds 
into the flight, whereupon this restriction is relaxed. When orbit velocity (8 km/sec) is 
attained, the boost reference angle has fallen to 19.5°, and the power beam is switched 
off. This launch trajectory (Case I) will place the vehicle into a highly elliptical orbit 
that might cause the vehicle to re-enter the atmosphere before one revolution of Earth 
is completed. Hence, this launch option may require the addition of a small chemical 
rocket, as shown in Fig. 11.7, for orbit circularization. The Case I launch geometry, 
with the chemical “kick” rocket, should enable circular satellite altitudes of up to 
2000 km. 

For the last three trajectory simulation cases, a different laser power transmission 
geometry and steering algorithm were used (see Fig. 11.8). As before, the LTD is 
ejected by the compressed air cannon at an initial velocity of 100 m/sec, density 
altitude of 3 km and initial launch angle (measured from horizontal) of 30°, 45° or 
60° for Cases II, III, and IV, respectively. Again, the laser powered boost begins 
0.5 seconds later, and the initial LTD pitch angle is maintained constant until the 
airbreathing engine phase terminates at Mach 5-7 and 100-115 kft density altitude. 
Thereafter, the craft is allowed to coast in a ballistic trajectory until apogee is reached 
(« 150 km), whereupon the craft is pitched over to roughly 0°, measured with respect 
to the horizon (see Fig. 11.9). At this point, the GBL is turned back on, and the beam 
is bounced off a laser relay satellite (200 km altitude orbit) and redirected toward the 
LTD, as shown in Fig. 11.8. Upon arrival of the beam, the LTD’s laser heated rocket 
mode is engaged, and the LTD accelerates to orbital velocity (8 km/s) in a matter of 
80-100 seconds of firing time. Currently, the steering model works only for satellites 
passing directly over the launch site, and further improvements should permit the 
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Figure 11.6: Boost Reference Angle vs. Time for Case I 



Figure 11.7: Chemical “Kick” Rocket Option (for orbit circularization) 
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inclusion of relay satellites with a more general ground track. 


11.2 Simulation Parameters 

The preliminary weight breakdown for the LTD can be found in Chap. 2. The SORT 
program can model variations in the center-of- gravity (CG) location, but the present 
analysis placed the vehicle CG on the centerline (along with the engine) to simplify 
the thrust-pointing requirements. 

A complete Earth gravity model exists in the SORT code which includes J2, J3, 
and J4 effects. The program continuously calculates the orbital parameters through- 
out the boost trajectory. 

The atmospheric data loaded into SORT was the ARDC Model Atmosphere [3]. 
Although wind can be a significant aerodynamic perturbation, none was included at 
this stage of the analysis. 

The SORT program treats the trajectory as a specified sequence of maneuvers of 
phases. Each event is defined by the user at appropriate places in the simulation so 
that new inputs can be assembled into the trajectory. The launhh vehicle flies from 
event to event until, finally, it obtains the orbital parameters necessary to achieve .the 
desired orbit altitude. The craft then coasts until it reaches apogee, and performs 
a final orbit-insertion bum. The SORT program has iteration capability that can 
be used to optimize total propellant weight (or laser energy) used during ascent. A 
typical trajectory sequence for Cases II, III, and IV is described below. 


11.3 Trajectory Sequence 

The LTD is given an initial velocity of 100 m/s at a precise launch pitch angle, by a 
compressed air cannon. While accelerating within the cannon, a lightweight sabot is 
used to cradle the LTD and also to prevent the escape of compressed gasses around 
the LTD. After emerging from the barrel, the sabot quickly separates from the LTD, 
leaving it free and clear — ready to accept the laser power beam. 

At 0.5 sec, the shroud lift ERH thruster is engaged and the LTD climbs out at a 
fixed boost reference angle. Since the ERH thruster mode is completely airbreathing, 
only a small amount (% 10%) of the liquid propellant is consumed as expendable 
coolant for the primary mirror and engine hot sections. This mode effectively initiates 
the acceleration portion of the flight. As the vehicle climbs out, the pulse frequency 
of the ERH engine is increased to accelerate the craft though the transonic region, 
toward Mach 5. As the ram drag increases, net effective thrust of the ERH thruster 
decays. At 100-115 kft, the engine is shut off. 
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Figure 11.8: Use of Laser Relay Satellite for Orbit Circularization (Cases II, II and 
IV) 
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Figure 11.9: LTD Pitch Angle vs. Time (for Cases II, III, IV) 


The vehicle continues to coast upward along a ballistic trajectory toward apogee. 
At the specified (or optimized) altitude, the vehicle pitches over just prior to receiving 
power from the low altitude relay satellite. The exact zenith angle (measured from 
vertical at the GBL site) at which to engage the relay satellite is determined itera- 
tively to achieve maximum performance. Severe atmospheric absorption is suffered 
by high power laser beams transmitted beyond zenith angles of 70.5°; i.e., three ver- 
tical atmosphere path lengths (at a 60° zenith angle, two atmospheric path lengths 
are incurred). Hence, for a 200 km relay satellite, the maximum slant range (GBL to 
relay satellite) is ^ 600 km; this equates to an orbital path length of about 565 km 
to the 0° zenith position. At 8 km/s, the relay satellite will take 70 seconds to arrive 
at a position directly over the GBL. 

The rocket mode is started as the vehicle reaches apogee. This bum is terminated 
as the inertial velocity reaches the necessary value for circular orbit at that altitude. 
The available laser power of 250 MW, translates to a rocket thrust of ^ 17.5 kN 
( CC = 70 N/MW), and a firing time of 80-100 seconds to reach 8 km/sec orbital 
velocity. Also, as mentioned earlier, the LTD optical axis must be aligned with the 
relay satellite to an accuracy of ± 0.5° or better. The LTD launch window must 
be exactly synchronized such that the LTD is approaching apogee when the relay 
satellite is at a zenith angle of 70.5° (or 60°, etc.). 

If the 200 km relay satellite is engaged when it is too high overhead (e.g., zenith 
angle of 0°) and the LTD is still at a much lower altitude (e.g., 150 km), the rocket 
mode will drive the LTD back into the atmosphere before it has a chance to accelerate 
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to orbital velocity. The trajectory will not climb high enough to prevent increased 
aerodynamic drag, which is associated with additional expenditures of fuel or laser 
energy, from terminating the mission. 

The final velocity desired at rocket termination is a function of the relay satellite 
and LTD apogee positions, altitude, LTD flight path angle and the drag experienced 
(i.e., AF) as the vehicle coasts up to the desired apogee altitude to make the final 
rocket bum. There is no easy analytical solution to this targeting problem. Many 
trajectories are run in an optimizing attempt. Typically the performance indicator 
“propellant weight” was optimized, although total laser energy could also be run. 
Analysis in the future will focus on the minimization of laser power. The key pa- 
rameters which were varied in the attempt to optimize propellant weight were the 
initial position of the relay satellite, power available to the engines (throttle setting), 
and apogee altitude for the pitchover maneuver. Trajectory optimization is a never- 
ending task (e.g., engineers are still fine-tuning Space Shuttle trajectories to discover 
more payload capability), but the results so far give a baseline indication of the LTD 
microspacecraft performance potential. 


11.4 Results 

The trajectory analyses to date using SORT have just barely scratched the surface of 
the comprehensive performance optimization problem. However these simulations do 
give highly accurate answers for the selected models and input data. Assuming the 
vehicle models (engine, drag, etc.) are fairly accurate, the resultant performance re- 
ported here would certainly be attainable. Although not optimal, the current results 
indicate the types of problems to be encountered, and present conservative perfor- 
mance capabilities (within the accuracies of the models and assumptions). Further, 
more detailed analyses can be expected to improve upon these results. 


11.4.1 Case I 

As mentioned above, Case I doesn’t use a laser relay satellite, and was run at a 30° 
initial launch angle with the craft switching from the “shroud lifter” ERH thruster 
mode directly to laser-heated rocket mode at about 65 seconds into the flight. This 
transition occurred at about Mach 7 and 35 km altitude. The rocket specific impulse 
was set at 875 seconds, or in the middle of the range projected for nitrogen propellant 
(e.g., 725—1025 sec). Figures 11.10 through 11.17 show the projected performance of 
the LTD machine for this case. Note that the rocket “bum” lasts for approximately 
85 seconds, and accelerates the vehicle to 8 km/sec at t = 150 seconds. The trajectory 
in Fig. 11.12 indicates that this point is reached at a down range distance of 500 km 
and an altitude of about 185 km, where the boost reference angle is 19.5°. If the 
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Figure 11.10: Thrust vs. Time (Case I) Figure 11.11: Altitude vs. Time (Case I) 


LTD is equipped with a small chemical rocket for orbital insertion (see Fig. 11.7), 
this mode could possibly enable circular orbits of up to 2000 km in altitude. Note 
in Fig. 11.14 that acceleration varies from 2 to 3 G’s in the airbreathing mode, and 
reaches a peak of 10 G’s at the end of the rocket mode. In Fig. 11.16, the peak power 
requirement is 350 MW which occurs in the ERH thruster mode and significantly 
exceeds the 250 MW baseline power level set for the rocket mode. 

The coupling coefficient reaches a peak of 525 N/MW early in the airbreathing 
mode and falls to a low of 70 N/MW in the rocket mode (see Fig. 11.17). 



Figure 11.12: Altitude vs. Downrange Distance (Case I) 
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Figure 11.13: Flight Mach Number vs. Time 
(Case I) 



Figure 11.15: Velocity vs. Altitude 
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Figure 11.14: Acceleration vs. Time 
(Case I) 



Figure 11.16: Power vs. Time 
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Figure 11.17: Coupling Coefficient vs. Time 
(Case I) 

11.4.2 Cases II, III and IV 

As before, cases II, III, and IV utilize the laser relay satellite, and assume initial 
boost reference angles of 30°, 45° or 60°, respectively. These cases are useful for low 
Earth orbit missions (e.g., around 200 km) and depend upon the laser relay satellite 
for orbited insertion. Figures 11.18 through 11.25 present the LTD performance for 
these cases. The rocket specific impulse was set at a conservative value of 650 seconds 
which allowed a coupling coefficient of 125 N/MW (see Fig. 11.25). 

Note from Fig. 11.18 that the unpowered coasting phase lasted from 100-150 sec- 
onds in duration, which extended the total launch duration to 250—300 seconds total. 




Figure 11.18: Thrust vs. Time 
(Cases II, III, &IV) 
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Figure 11.19: Altitude vs. Time 
(Cases II, IE, &IV) 
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Figure 11.24: Power vs. Time Figure 11.25: Coupling Coefficients vs. Time 

(Cases II, III, &TV) (Cases II, III, &TV) 


As indicated in Fig. 11.20, Case II (30°) resulted in a retrograde orbit and would 
immediately re-enter the atmosphere; in contrast, Case III (45°) produced a circular 
orbit at about 175 km altitude, and Case IV (60°) — a near-circular orbit at 265 km. 
From Fig. 11.23, it is apparent that, although Case III suffered the greatest AV loss 
due to drag in the coasting phase, 60° appears to be an excellent launch angle. Fig- 
ure 11.24 indicates that the peak airbreathing engine laser power required for the 60° 
case is about 310 MW; for the 45°, 330 MW. 


11.5 Summary 

The trajectory analysis will need to continue as the engine and aerodynamic models 
evolve and the final results may end up considerably different than those presented 
here. This does not imply that the current analysis is crude or inaccurate, because it is 
quite precise within the confines of the models and inputs. The predicted performance 
for the Lightcraft Technology Demonstrator is undoubtedly revolutionary, and further 
analysis will continue to define this new frontier. 


i 



11.5. SUMMARY 


195 


REFERENCES 


[1] Beming, M., “Version IV User’s Guide for the Simulation and Optimization of 

Rocket Trajectories Program,” McDonnell Douglas Astronautics Co., Houston, 
TX, Transmittal Memo No. 1.2-TM-FM86028-53, December 1985. 

[2] Hoemer, S.F., Fluid Dynamic Drag , Published by the Author, 1965. 

[3] Hill, P. and Peterson, C., Mechanics and Thermodynamics of Propulsion. Addison- 

Wesley Publishing, 1965. 



196 


CHAPTER 11. EARTH-TO- ORBIT TRAJECTORY SIMULATION 



Chapter 12 

LIGHTSAT SYSTEMS DESIGN 


As mentioned earlier, the Lightcraft Technology Demonstrator (LTD) is a microspace- 
craft that will not only demonstrate laser propulsion technology, but also serve as a 
lightsat. After attaining orbit, the LTD will become a low altitude (e.g., 200 km - 
2000 km) sensor satellite, and have a nominal pointing accuracy of ± 0.5 deg (see 
Chapt. 6). Satellite hardware must be chosen to fulfill the LTD’s role as a lightsat; 
i.e., weight, cost and size of electronics, attitude control systems, power supply, etc' 
must be minimized. 


12.1 Lightsat Background 

Small, lightweight (several kg), unmanned missions for specific applications have re- 
ceived increased attention as both equipment and launch costs continue to rise [1], 
This was one of the underlying motivations behind the development of “get-away 
special” canisters for the Space Shuttle. A recent proposal by NASA and the Univer- 
sity Corporation for Atmospheric Research involves housing scientific apparatus in 
the unpressurized sections of discarded Space Shuttle external fuel tanks [2], Experi- 
ments could be housed in a cylindrical section 22.5 ft. long and 27.5 ft. in diameter; 
the apparatus would function like a high-flying sounding rocket to collect data for 
about an hour in suborbital space before the tanks re-enter the atmosphere and bum 
up. The program, known as Space Phoenix, is intended to recover some portion of 
the $40 million cost for each shuttle fuel tank. Efforts such as these illustrate the im- 
portance of using small components for multiple mission purposes in order to provide 
economical access to space. 

Many significant experiments may be performed in space without resorting to 
large, relatively expensive payloads. These so-called “lightsats” or “microspacecraft” 
have been considered for many diverse applications, including communications satel- 
lites, lidar systems, and satellite laser range finding. Several less obvious applications 
have been presented in Ref. [1, 2, 4], Specifically, NASA’s proposed Mission to Planet 
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Earth calls for a global network of small surface or atmospheric sensor satellites to 
monitor globed climate trends such as ozone depletion, pollution, or continental drift. 
Another NASA program, Project Pathfinder, has suggested the use of many mi- 
crospacecraft at different locations to act as distributed (phased) arrays for either 
radio or optical signals. Such an array directed towards space could provide the most 
accurate triangulation measurements of our galaxy ever made by astronomers. As 
microelectronics technology continues to reduce payload size and weight, lightsat sys- 
tems have been recognized as the means to fill a vital market niche in space science 

[4]- 

Although the missions described above do not specifically require laser boosted 
spacecraft technology, it has been noted [4] that the development of appropriate 
launch systems is important to the success of lightsats. With these applications as 
motivation, both lightsats and laser propulsion are mentioned as future alternatives 
in the report of the National Commission on Space [5]. Small, multi-purpose payloads 
may also encourage the involvement of universities and the private sector in space 
technology (similar to the get-away specials) as the launch costs decrease. Although 
cost projections vary, Canavan [6, 7] has provided some interesting speculation in 
this regard. In a recent report on the laser deployment of sensor satellites [6], it 
is assumed that a ground based laser system could launch 1 kg per 1 MW of laser 
power, and that sensor satellites could be produced for as little as $1, 000/kg (ad- 
vanced sensor hardware, however, can cost up to $200, 000/kg [4]). Taking the lower 
estimate, this implies performance of $ 1,000/MW; a 10 MW laser facility is projected 
to cost $100 million [6, 7], and would be able to launch a $10,000 satellite under these 
assumptions. In order to recover its initial cost at this rate, the facility would have to 
launch 100,000 kg total payload. By comparison, a Space Shuttle flight to low earth 
orbit (250 km) costs about $3000/kg [7], so that using the shuttle to launch the same 
total weight would cost about $300 million - roughly three times the amortization cost 
of a ground based laser facility! Also, we have not made any assumptions concerning 
the launch rates required; a four shuttle fleet may launch 14 flights in a single year 
at best [6, 7], while the lightsats could launch at a much higher rate. It has been 
noted [7] that over three years, the laser facility would amortize to about $l/second 
of launch time, making the technology highly competitve. 


12.2 Orbital Positioning 

For an Earth orbiting satellite to accomplish its objectives, it must maintain a specific 
attitude in relation to the planet. The spacecraft must also be able to determine its 
own position in orbit. Finally, it must then have the ability to correct its attitude, 
and possibly position, from sensed orbital deviations. 
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12.2.1 Attitude Determination 

Since the LTD is to be a low Earth orbit (e.g. 200 km - 2000 km) satellite, a pointing 
accuracy of ± 0.5 degree will be sufficient for the mission. Attitude sensors which 
have been considered are star trackers, sun sensors, and Earth horizon sensors. 

Star trackers employ a detector which searches for visual signals above a given 
threshold value and then tracks the signal until the star is lost, or until it resumes 
searching after a predetermined time interval. Although these sensors were origi- 
nally considered for use in the LTD, after further investigation this type of sensor was 
deemed too accurate for mission requirements in relation to the cost of the component. 
The crudest model available is accurate to one arc minute, and its use would necessi- 
tate periodic compilation and uplinking of an ephemeris 1 to the satellite. Therefore, 
a star tracker would not be cost efficient. 

Another useful method of attitude determination is the sun sensor. In this cate- 
gory there are three basic types: 1) an analog sensor, which continuously calculates 
the sun angle; 2) a sun presence sensor, which calculates the sun angle when the 
sun is in the field of view; 3) and a fine sun sensor, which is the most precise of the 
three. For LTD purposes, the analog sensor proved to be the most beneficial because 
continuous sun tracking is often useful, if not necessary, when utilizing a solar array. 
An example of an acceptable sensor is made by Lockheed Missiles and Space Com- 
pany, located in Sunnyvale, California. Their wide angle sun sensor (WASS) weighs 
0.155 kg and measures 6x5x5 cm. It is accurate to 0.5 degrees and has a 27 t steradian 
field of view. Using the most popular output of 620 fiA, almost no power would be 
consumed. The voltage type of output would also use minimal power. 

The third type of attitude sensor is the Earth horizon sensor. It uses i nfra red 
sensors to scan the Earth and measure rotation angles, thus defining the spacecraft 
attitude relative to the Earth from a constant altitude. The detector senses radiation 
between 13.5 to 25 /im in the spectral band. A small sun detector is also utilized 
to identify sunlight intrusion. A circular orbit requires one horizon sensor while an 
elliptical orbit needs two sensors. These attitude sensors are cheaper than star trackers 
and do not need the added expense of an ephemeris. A representative example is made 
by Ithaco, located in Ithaca, New York. Their horizon crossing indicator (HCI) has a 
head that is 10.0 cm long and 5.9 cm in diameter with an electronics box measuring 
12.8x10.3x4.3 cm. The combined weight of the system is 0.65 kg. The HCI has an 
accuracy of 0.1 degree and a power requirement of less than 0.7 watts. 

Therefore a combination of an analog sun sensor and one or two Earth horizon 
sensors will be used in the LTD for attitude determination. The analog sun sensor 
will be the primary device while the Earth horizon sensor will be used for back up. 


1 An ephemris is an astronautical almanac or table of predicted positions of celestial bodies. 
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12.2.2 Orbital Position Determination 

In addition to the attitude sensors for solar array positioning and antenna point- 
ing, the satellite will need a method for determining its orbital position above the 
Earth’s skies. Three co mm on methods will be discussed: 1) the Global Positioning 
System (GPS); 2) the Tracking and Data Relay Satellite System (TDRSS); 3) and 
the Satellite Tracking and Data Network (STDN). 

GPS employs several satellites, each of which continuously broadcasts its own 
ephemeris and time. A GPS receiver can then determine its own position from four 
separate simultaneous GPS observations. The system can provide receivers with a 
positional accuracy of 20 m, velocity accuracy of 0.06 m/s, and time accurate to 
10 ms. It was determined, however, that the hardware required weighed far too much 
for the LTD. 

.Although TDRSS and STDN are primarily communications options 2 , they each 
can be used for determining a satellite’s orbital position. TDRSS has a tracking 
system much like that of the GPS, except that it involves a real time cost since 
calculations are done by TDRSS and not by one’s own satellite. The TDRSS consists 
of only two satellites and one receiving station on the ground (see Figure 12.1). The 
STDN is a network of ground stations (see Fig. 12.2) used mainly for communication 
purposes, but also from which ranging can be performed. Ranging is the process 
of sending radar waves to a satellite and having it reflect the waves back to the 
station. Knowing the wavelength and the velocity of the wave and the time it takes 
to travel to and from the satellite, the range (position) and range rate (velocity) can 
be determined. Both the TDRSS and the STDN require only a lightweight antenna 
as hardware. 

Therefore the LTD will use the TDRSS as its primary means of determining its 
orbital position and velocity, because of its simpler system. When the orbit is less than 
1200 km and the satellite falls within TDRSS’s zone of exclusion (see Figure 12.3), 
STDN will be used. In any case, STDN will be used as a backup. 

12.2.3 Attitude and Orbital Control 

Once corrections are determined to be needed, whether because of deviations or com- 
mand decisions, there are three methods to change the attitude of a satellite: mo- 
mentum reaction wheels, magnetic torque bars, and thrusters. For communication 
applications, only the yaw and pitch axis need be controlled since the LTD is axi- 
symmetric about its roll (longitudinal) axis. For satellite photography and remote 
sensing, however, “roll” stabilization about the longitudinal axis must be provided. 

2 Using TDRSS k STDN would give the added bonus of easy satellite housekeeping downlink and 
command uplink. 
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TDRS-EAST 

LONG. 41* WEST ANTENNA STEERING TDRS-UEST 



Figure 12.1: Tracking and Data Relay Satellite System 

Finally, for orbital translations, only thrusters will work. 

Momentum reaction wheels would not be a feasible option for the LTD because, 
although some of the more modem ones are small, they are still too heavy for a 
lightsat. 

Ithaco of Ithaca, New York, manufactures a magnetic torque bar system called 
Torqrod. The Torqrod uses the readings from a magnetometer (which gives a three- 
axis measurement of the Earth’s magnetic field) in determining which rod to energize 
to produce the desired attitude correction. Torqer bars, used in the Torqrod system, 
are inexpensive, weigh as little as 0.375 kg each, and measure less than 40 cm long. 
Two Torqer bars are needed; one in the yaw axis and one in the pitch axis. The 
major advantage of employing the Torqrod system is that the mission life becomes 
component dependent, because the system uses very little power and no consumables. 
The Torqrod system, however, compensates slowly. 

Cold gas thrusters should be used for attitude control on the LTD for several 
reasons. First, although hot gas thrusters generate greater specific impulse, cold 
gas thrusters operate more consistently, with greater precision and are substantially 
cheaper. Also, no extra fuel need be lifted to orbit by the LTD because 2 kg of cold gas 
helium (and some LN 2 ) will already be on board, having been used to pressurize the 
liquid nitrogen propellant tank prior to orbit. Only four thrusters will be necessary, 
two in each direction along both the yaw and pitch axis. 



ORKINAt PAGE IS 

OF POOft QUALITY 


CHAPTER 12. LIGHTSAT SYSTEMS DESIGN 



BUR Ajjf TAN 


! 1 ' 41 1 T^ 1 

120° 


40° 

1 

40° 


120° 1 


lonC'tuoi to 


ASCtNSlON 'SL ANO 1 »CNI 

V~ USi 

- l*sr 1 f 33 

3«S«40 r 79 »T‘ 

Santiago Ch'lI f*GO' 

V- US* 

VM< G A A A M 

IN T ( M* { «0M( T t a 

- 33*09 03 SJ 

- ]]«M 04 OC 

- ) r°o* s* »a 

799° 70 01 0* 
78*°7COl 07" 
Kr 

• (•muO* '>0*1 

V- US* 

1 A 0-4 H aO A« 

1 

3 7*?i OS 00 
37*70 S3. ON 

?9S®?0 3* *4 
799*704) W 

C"*NC Can A« Y 
• SLANO 'C**' 

V- us* 

1 1 

;t® 4* s' S' 

3*4*71 91 M 

t NG'N { f MtNG TMAiN'NG 
C(N11M MAM*L AND ■( TO 

*- uss 

V- US* IIMTSl 
>NT f m * f MOMf T [ » 

jr»s> S4 m 

JfM V4 0 * 

1 )>°s« s» n 

793*09 74 73 
793*09 79 7' 
793*09 31 i* 

GOiOSTONt 
C*Li» Q«N< a IGOV 

}br~ US* 

V- US* i MTS- 

3 S f 70 n u 

3S C ?0 79 M 

743*01 3S OS 
7* 3*0 T J) 4S 

Guam iGvnm 

v- us* 

1 JO'S 3* 7S 

1444*4 I 7 S3 

HAMAH |Htw 

v« ust 

* AS <1 MADAM 

1 r.-^CT 3* «4 

j 77*0* 74 3* 

700® 70 ON 43 
700® 70 04 c: 

maom'O $/ain <maoi 

TV- USS 

! *C® 7 7 H ST 

ns*49 SJ S* 

ut mmttt island 
KOA'O* 'M'l' 

V-* US» ho 1 
V- US* NO 7 

1 

1 79° JO 79 T9 

tiojojt*’ 

1 

779®'* 73 rs 
1 J T9* i ( 73 C. 

AuSTMA^tA >0° *> ' 

INM M» | M OM ( It A 

i 3S® 3 T 37 '9 

1 

! '*9*ST 'S IS 

Quito fcu*oo» 'Ou»» 

INTfM# |PON*|T|M 

t 

- OCJOJT 77 04 

7*'®29 IS '0 

«OSman 

4 J*** uS* 

[ ]!<"'«« 

I T 7*07 74 94 

NQM Th C AROl 'N* IMOS- 

VMl GMARM 

j IS®' * 47 07 

337*07 ?* *) 


Ml *U*L 'C IAN 

VMl G* AMR 1 

INTfM! {MOM( T E a 

7 MS is MAO* M 

'9*0' 'S >* 

- t9°00 3' W> 

- >9°00 OS ST 

1 47® 1* 11 93 

*)®i 7 S9 n 
47®»| S3 44 

#A'M*AN*S ALASKA 

v us* 

M°S» 19 70 

712*79 13 39 

<Ul *> 

VMl gmamm 
INTfM* (ROME T I ■ 

M°S* I T SO 
<m*°SI 3S 9* 

717*79 '9 '? 
7i 7*7* 3' *9 

MMNKflfLO IHGSAHO 

iNTEMf (ROME Tf M 

S'°7S 44 1? 

3S9*1* 09 13 


»*OV| 

( LU^SOlQ 

iwitCM) 


* KMNlrt SIGN (-1 iNO'CATfS SOUTH LATlTuOl 


Figure 12.2: Locations and Geodetic Coordinates of STDN Stations 
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Although orbital translations are not considered in the mission parameters, it has 
been suggested that the injector heads could be used in emergency orbital maneuvers. 

irst the satellite would be pointed in the right direction, then the pressurized helium 
would enter the injector (see Chapt. 4) heads by way of the liquid nitrogen tank. 

summary, for primary attitude control, a magnetic torque bar system is chosen 
with cold gas thrusters being used when immediate corrections are needed. For any 
unforeseen translational motion required, the thrust could be provided by the injector 
heads used as thrusters. 


12.3 Communications and Sensing 

An effective sensor lightsat must be able to communicate with Launch Command 
to transmit its status, to send data it has gathered (from another source), and to 
receive command instructions. To accomplish this goal, specific orbital parameters, 
communication frequencies, and antenna options must be explored. 


12.3.1 Orbital Parameters 

There are three major orbital parameters associated with sensor satellite missions: 
) coverage area, or the portion of the Earth’s surface that can be viewed by the 
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Satellite 



Figure 12.4: Satellite- Earth Angles 

satellite, 2) slant range, or the distance from a fixed point on the Earth to the satellite, 
and 3) length of time a satellite can observe with a prescribed elevation angle (i.e., 
“window of opportunity”). The elevation angle is essential since a sensor satellite’s 
“vision” can be impaired if it is at too low an elevation angle, close to the horizon. 
An example of satellite-Earth angles is shown in Figure 12.4 (<f>i = elevation angle 
from ground station, A cov = shaded area)(see Ref. 8). 

The coverage angle, A ccw , with an elevation angle, <£i, can be obtained from: 

A cov = 2 itR 2 e (1 - cos<f> E ) (12.1) 

where, 

<j> E = arccos f + h^j — <f> l (12-2) 

R e is the radius of the Earth; and h is the altitude of the satellite above the Earth. 
Since the Earth’s total surface area is 4x R E , A cov can be rewritten as a fraction of 
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the Earth’s total surface: 

= 0.5(1 — cos<j> E ) (12.3) 

For example, with a <f> E of 0° and altitude of 1000 km, the maximum percentage of 
the Earth that the LTD can view is 6.78%. 

The slant range between a point on Earth and a satellite at altitude h and angle 
<f > i is calculated from: 


z = \J (Re sin 0x) 2 + TRsrh + h 2 — R E sin <j> i 


(12.4) 


Before the LTD’s sensor “window of opportunity” can be determined, one must first 
find the lightsat’s orbital period t,. Based on Kepler’s 3rd law 


2iry/aA 

t,= vt 


(12.5) 


where a is the semi-major axis of an elliptical trajectory, and fj, is the product of the 
universal gravitational constant G and the sum of the mass of the Earth m. E and 
the satellite m,. Since the mass of the satellite is negligible compared to that of the 
Earth, 

fj. = G(m E + m.) = Gm E (12.6) 

For an orbit of 1000 km, Earth radius of 6378.155 km, and /i of 3.986013 x 10 s km 3 /s 2 , 
the LTD’s orbital period (t,) is approximately 105 minutes. 

For a satellite at altitude h, the “window of opportunity” (t p ) is 


t 


v 




(12.7) 


where t E is the rotation period of the Earth (one sidereal day, which equals 23 h, 
56 min, and 4.09 s, or 86164.09 s). The ± sign depends on whether the satellite is 
in a prograde (same direction with respect to Earth) orbit, or a retrograde (opposite 
direction) orbit. For a <f> E of 30.18° and an altitude of 1000 km, t p is approximately 
72 s. 


12.3.2 Communication Frequencies 

The typical frequencies allocated for satellite communications fall in the super high 
frequency (SHF) and extremely high frequency (EHF) bands (see Table 12.1). The 
frequency spectrum is divided into subbands which are listed in Table 12.2. The 
LTD’s communications system will be chosen for low power requirements, mi nimal 
propagation distortions, and reduced noise and interference effects. The transponder 
must also be designed to function both as a receiver and transmitter for TDRSS and 
STDN purposes. 
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Table 12.1: Electromagnetic Frequency Spectrum 


Frequency 

Wavelength (m) 

Designation 

3 Hz- 30 kHz 

10 8 - 10 4 

Very Low Frequency (VLF) 

30 kHz- 300 kHz 

10 4 - 10 3 

Low Frequency (LF) 

300 kHz- 3 MHz 

10 3 - 10 2 

Medium Frequency (MF) 

3 MHz- 30 MHz 

10 2 - 10 

High Frequency (HF) 

30 MHz- 300 MHz 

10- 1 

Very High Frequency (VHF) 

300 MHz- 3 GHz 

1- 10' 1 

Ultra High Frequency (UHF) 

3 GHz- 30 GHz 

10" 1 - 10~ 2 

Super High Frequency (SHF) 

30 GHz- 300 GHz 

10~ 2 - 10~ 3 

Extremely High Frequency (EHF) 

10 3 GHz- 10 7 GHz 

3 X 10" 5 - 3 x 10~ 9 

Infrared, Visible Light, 
and Ultraviolet 


Table 12.2: Satellite Frequency Spectrum 


Frequency 

Rang< 

(GHz) 

L 

1 - 

2 

S 

2 - 

4 

c 

4 - 

8 

X 

8 - 

12 

Ku 

12 - 

18 

K 

18 - 

27 

Ka 

27 - 

40 

Millimeter 

40 - 

300 


I 
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Sidelobes 



Figure 12.5: Antenna Gain Pattern (see Ref. 8) 


12.3.3 Antenna Options 

The most important antenna characteristics are its gain (maximum value of the gain 
pattern), beamwidth (a measure of the angle over which most of the gain occurs), and 
sidelobes (amount of gain in off-axis directions). Figure 12.5 shows a typical antenna 
gain pattern, with its essential features. The most desirable antenna patterns are 
those with high gain concentrations over a narrow beamwidth and small sidelobes. A 
list of typical antenna patterns are shown in Table 12.3. The LTD will be equipped 
with an erectable light-weight antenna for communicating with launch command, 
accessing the TRDSS (and STDN) network, and relaying data gathered from LTD 
sensors (see Fig. 12.6). 

The maximum gain, g, and half power beamwidth, <f> b , in radians are given as: 


9 = 



( 12 . 8 ) 


4>b ~ 


A 


(12.9) 


where A is the antenna area, d. is the aperture diameter, and p a is the antenna 
efficiency factor. 


12.4 On Board Computer 

The on-board-computer (OBC) will be located in the electronics payload compart- 
ment of the LTD. It will be encased in a radiation-shielded and thermal protective 
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Antenna tvpe 


Half power 
beamwtdth 


Long dipole 


dimension d 


Horn 

direction 






r 

17 6 dB i W 

Circular 

reflector ! 

( ° 

o , =22->m, 

i 

\x 

Y 



Reflector 

direction 


Table 12.3: Typical Antenna Gain Patterns (see Ref. 8) 
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LIGHTWEIGHT ANTENNA EXTENDED 
BY TWO ASTROMASTS 


Figure 12.6: LTD will deploy an antenna for communications during satellite mode. 

casing. The OBC will control the specific actions of all electrical and mechanical 
components. It will manage the antenna, transponder, temperature sensors, solar 
array, batteries, attitude system, and LN 2 propellant delivery. 

The OBC will begin operation at liftoff and remain in operation throughout the life 
of the LTD. For example, in the boost mode temperature sensors that are distributed 
at critical areas on the mirror and thrust surfaces will advise the OBC so that adequate 
coolant flow rates are administered. Later, once orbit is achieved, the OBC will direct 
the LTD forebody to open up into four “petals,” expose the solar array’ and extend 
the antenna. These “petals” will be actively positioned by mechanical actuators that 
are programmed (and controlled by the OBC) to maxi miz e solar power collection (see 
Fig. 12.7). 


12.5 Power Systems 

The LTD will primarily be powered and recharged by solar energy gathered by an 
array of solar cells. During launch and when the satellite is in the Earth’s shadow, 
however, the LTD will have to rely on its batteries to power its electrical systems. 
Therefore, to design an adequate power system for the LTD, several areas must be 
considered; the shadow factor, solar array, and rechargeable batteries. 



210 


ORIGINAL PAGt 

black and white photograp 


CHAPTER 12. LIGHTSAT SYSTEMS DESIGN 



Figure 12.7: Solar cell array exposed during LTD satellite mode. 

12.5.1 Shadow Factor 

In Low Earth Orbit (LEO), a certain portion of the LTD’s trajectory will remain 
unilluminated by the sun. Clearly, the craft will receive no solar energy while in the 
Earth’s shadow. The eclipse period can be calculated from the maximum shadow 
angle: 

4> max = 180° - 2arccos (12.10) 

For example, at an altitude of 1000 km, (f>max is 119.6°. Hence, 1/3 of the satellite s 
time is spent in the Earth’s shadow. With a period of 105 minutes, the satellite would 
be in the sun for 70 minutes. 


12.5.2 Solar Array 

Since the LTD will be illuminated by the sun for only 2/3 of its orbital period, the 
actual power requirement must be increased by 50% to a total of 60.8 W from an 
initial requirement of 40.5 W (see Table 12.4). By allowing for an additional 25% 
safety margin, the solar array must be able to generate approximately 76 W of energy. 


i 
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Table 12.4: Power Requirements for Solar Array 


Solar Array Actuators: 

20.0 W 

Horizon Sensor: 

0.5 W 

On Board Computer (OBC): 

5.0 W 

Regulators- (battery, gas /pressure), 
Sensors (temperature, optical), and 
Heaters: 

15.0 W 

Total: 

40.5 W 


The solar array surface area can then be calculated from the following equation: 


A = 


P 

S *tj * F * cos r 


( 12 . 11 ) 


where P is the power load requirement (76 W), S is the average solar intensity' (1353 
w/m ), r is the compound angle between the sun line and the average array normal 
(45°), T) is the solar cell efficiency (0.1771 at 64° C), and F is the sum-total of all array 
design and degradation factors (0.8). The compound angle, T, is an approximation of 
the average power angle that occurs over the array surface. The solar cell efficiency 
is obtained from: 

-8{T -T nf )\ (12.12) 

where is the initial solar cell efficiency (0.22 for a silicon solar cell), (3 is the 
thermal calibration coefficient (0.005° C), T is the average temperature on a flat 

solar panel at 1000 km (64° C), and T n j is 25° C. The needed solar array surface 
area is found to be 0.56 m 2 . 


12.5.3 Battery 

The LTD battery must be rechargeable and be able to store enough energy to supply 
power needs through out the 35 minutes of darkness. From an initial power require- 
ment of 40.5 W, less power needed for the solar cell actuators, approximately 20 W 
will be required during this period of darkness, which equates to a 11.7 W-h battery. 
Allowing for a 80% depth of discharge, a 20 W-h battery should suffice. 

Microprocessors commonly operate at 5 ± 0.5 V, or about 6 V. The use of Nickel- 
Hydrogen batteries would be a logical choice, since they provide much greater energy 
densities, and volume reductions than conventional battery types. A weight of 0.5 kg, 
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length of 10 cm, diameter of 9 cm, and a capacity of 1 A-h are likely battery charac- 
teristics. 


12.6 Conclusion 

As mentioned above, once in orbit, the LTD will function as a light-weight sensor 
satellite, complete with all the requisite on-board electronics, communications, and 
attitude control systems. The LTD senses its attitude with a sun sensor and a Earth 
horizon sensor. To change attitude, it uses a combination of cold gas jets and/or 
magnetic torque bars. It will be equipped with an erectable, light-weight antenna for 
communicating with Launch Command, accessing the TRDSS (and STDN) networks, 
and relaying data gathered from LTD sensors. An on-board-computer will manage 
satellite electrical, mechanical and attitude control systems, coolant flow rates, and 
solar array actuators. The solar array is mounted onto the inner forebody surface, 
which opens into four “petals” to collect solar energy for recharging the on-board 
nickel- hydrogen battery. The battery is used primarily when the satellite traverses 
the Earth’s shadow. This combination of satellite components is chosen for minimal 
cost and weight, so that the LTD can fulfill its primary objective of a self-launched 
lightsat. 


REFERENCES 

[ 1] French., J. R., “Very Small Spacecraft for Planetary Exploration missions,” USRA 
Lightsat Conference (1988). 

[ 2] Bell, T. E., “Jettisoned Shuttle Fuel Tanks to Harbor Atmospheric and Space 
Science Experiments,” The Institute (supp. to IEEE Spectrum), p 4, May 1989. 

[ 3] Browde, A., “Shrinking Electronics Boost Satellite Utility,” Aerospace American, 

pp 22-24, 1988. 

[ 4] Fuhs, A. E., Mosier, M. R., “A Niche for Lightweight Satellites,” Aerospace 
American, pp 14-16, 1988. 

[ 5] Cohen, B. N., editor, Pioneering the Space Frontier: The Report of the National 
Commission on Space , Bantam Books, Inc., N.Y., 1986. 

[ 6] Canavan, G. H., “Laser Propulsion of Sensors and Interceptors,” LA-UR-87- 
3208, Los Alimos National Laboratory, Los Alimos, NM, 1987. 

[ 7] Canavan, G. H., “Decoy Deployment with Laser Propulsion,” LA-UR-88-1246, 
Los Alimos National Laboratory, Los Alimos, NM, 1988. 


I 



12.6. CONCLUSION 


213 


[ 8] Gagiiardi, R. M., Satellite Communications, Lifetime Learning Publications, Bel- 
mont, California, 1984. 

[ 9] Ha, T. T., Digital Satellite Communications, Macmillan Publishing Company, 
New York, 1986. 

[10] Husson, J. C., “Introduction to Space Systems,” Electrical Communication, 
Volume 62, Number 1, 1988, pp 4-11. 

[11] Product Literature, Micro Mo Electronics, Inc., 742 Second Ave S., St. Peters- 
burg, Florida, 33701. 

[12] O’Lone, Richard G., “U.S. Planning New Emphasis on Lightweight Satellite 
Systems.” Aviation Week &z Space Technology, August 10^ 1987. 

[13] Patterson, R. E., Luft, W., and Dunlop, J. D., “Development of Spacecraft 
Power Systems Using Nickel- Hydrogen Battery Cells,” Paper 76-288, AIAA/CASI 
6th Communications Satellite Systems Conference, Montreal, Canada, April 5- 
8, 1976. 

[14] Rauschenbach, H. S., Solar Cell Array Design Handbook, Van Nostrand Reinhold 
Company, New York, 1980. 

[15] Product Literature, Solarex Corporation, 1335 Piccard Drive, Rockville, Mary- 
land, 20850. 

[16] Wertz, James R., Spacecraft Attitude Determination and Control,. 

[17] Product Literature, Whittaker- Yardney Power Systems, a division of Whittaker 
Technical Products, Inc., 82 Mechanic Street, Pawcatuck, Conneticut, 02891. 

[18] Williams, A. F., The Handbook of Photovoltaic Applications, The Fairmont 
Press, Inc., Atlanta, Georgia, 1986. 



214 


CHAPTER 12. LIGHTSAT SYSTEMS DESIGN 


I 



Chapter 13 


FUTURE DIRECTIONS 


For the following year, the Lightcraft Technology Demonstrator’s combined- cycle en- 
gine will remain as the central focus of further research efforts. The preliminary 
conceptual design of this laser-boosted microspacecraft was performed last year. The 
study involved an analytical performance analysis of the LTD engine in its two modes 
of operation, including an assessment of propellant and tankage requirements, and 
a detailed design of the inner structure and external body shape. Also, a structural 
vibrational analysis was performed on the annular shroud pulsejet engine. Finally, 
the sensor satellite mission was examined to identify the requisite subsystem hard- 
ware (e.g., electrical power supply, optics and sensors, communications and attitude 
control systems). 


13.1 Computer Modeling Efforts 

In the upcoming year, efforts will be focused on the LTD’s propulsion system and 
vehicle aerodynamics rather than on its mission as a sensor satellite. A 3-D compu- 
tational fluid dynamics (CFD) model of the LTD’s external inlet will be assembled, 
and run for a variety of different shapes and flight speeds. From these computer 
simulations, the boundary layer thickness, external surface pressure distribution, and 
transition point (from laminar to turbulent flow) can be determined. A 2-D blast 
wave code provided by the Naval Research Laboratory (NRL) may also be employed 
to model the impulsive thrust generated by laser-induced blast waves. In addition, 
a computer analysis of the radiation/ convection heat transfer will be conducted to 
find the engine/ vehicle thermal profile. All of this will enable the intelligent choice 
of appropriate spacecraft materials (specifically tailored to withstand the hypersonic 
transatmospheric environment) and also the design of active thermal cooling systems 
for the primary optics and engine hot-sections. Finally, a more detailed vibrational 
analysis will be performed using finite element models that will include not only the 
shroud/strut assembly, but also the spacecraft internal support structure. 
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13.2 Wind Tunnel Experiments 

Several experiments will be performed to validate theoretical vehicle and engine inlet 
performance simulations. Hot and cold flow wind and shock tunnel tests are planned , 
with two different scale models ( 5 in. dia. and 1.25 in. dia.). The test flow con- 
ditions will be varied from subsonic to hypersonic. Flow visualization techniques, 
including smoke injection and schlieren photography, will provide details of flow dy- 
namics around the vehicle, along with bow shock wave structure. Static and total 
pressures will be measured at critical positions on and around the test models. Heat 
transfer data will be obtained in high temperature flow experiments, using the RPI 
24 in. diameter hypersonic shock tunnel. 


13.3 Laser Impulse Experiment 

Laser impulse experiments will also be conducted on a full-scale, 1 / 15th segment of 
the annular engine, in both static and dynamic wind conditions. These high power 
laser experiments will yield engine impulse and heat transfer data to prove technical 
feasibility of the propulsion concept. The proposed experiment is designed to verify 
• the analytical performance for the airbreathing combined-cycle engine mode, referred 
to as the External Radiation- Heated (ERH) thruster (see chap. 9). 

The experiment is specifically designed for the one kilojoule pulsed CO 2 laser 
shown in Figure 13.1, which was loaned to Rensselaer by the Naval Research Labora- 
tory in Washington, DC. This laser is similar to the Lumonics 624-TEA commercial 
pulsed CO 2 laser. The high power NRL laser (now at RPI) will be set up with unsta- 
ble resonator optics to produce a laser pulse duration of 50 ns. The 20 cm x 20 cm 
output laser beam will be focused by a special parabolic mirror to intensities exceed- 
ing 5 x 10 8 W/cm 2 at the impulse pressure plate, as shown in Figure 13.2 (note that 
the laser-induced clean air breakdown threshold is greater than 10 fl W/cm ). The 
resultant LSD wave will be roughly 20 cm long and 0.5 cm in diameter. The impulse 
pressure plate pictured in Figure 13.2 will measure 25 cm x 25 cm and represent a 
full-scale, 1 / 15th. segment of the annular LTD engine. Finally, it should be noted that 
the four module NRL CO 2 laser will be set up to run at an effective PRJ 1 of 1200 
to 12,000 Hertz — for two pulses only. With the latter interpulse time, the effective 
time- average laser power will be 5 megawatts. Only two pulses are required to exactly 
simulate the propulsive physics of a repetitively-pulsed engine. 

Several statements regarding the diagnostics and instrumentation for the experi- 
ment are in order. First, schlieren photography and a TRW image converter camera 
with nano-second resolution will be employed to record the process of LSD wave igni- 
tion and the evolving cylindrical shock wave patterns as a function of time: framing 
rates of 50-100 nano-seconds are possible. Efforts will be made to ascertain the min- 
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imum flux and fluence required to start the LSD wave at a metallic “ignitor” plate. 
Also, the laser pulse shape will be correlated with temporal variations in the laser flux 
upon the LSD wave, and their effect upon the LSD wave velocity will be assessed. 

Second, the 25 cm x 25 cm ERH thruster plate will be instrumented with at least 
six pressure transducers, and a number of thin (and thick) film platinum heat gauges 
for measuring convective and radiative heat transfer. This data will be recorded on a 
4-channel Nicolet digital osdlliscope with nano-second resolution. Because of the high 
temperature engine environment, radiation losses from the cylindrical, high pressure 
air plasma must be studied. The chemistry and hydrodynamics of the pulsed thruster 
exhaust will be examined to assess the effects of “frozen flow” losses. 

Also, the gas dynamics of refreshing inlet air will be documented. These ex- 
periments will provide clear insight to the pulse-by-pulse interaction process during 
impulse generation. As mentioned above, the experimental apparatus is designed 
to convect air across the thruster plate to enable the determination of the optimum 
“wind” velocity for adequate “clearing-time” between laser pulses — such that max- 
imum impulse is delivered to the ERH thruster plate. (Inddently, to calibrate the 
pressure-integrated impulse data system, the “pendulum” method of direct impulse 
measurement will be used for the initial static tests only.) 

Portrayed in Figure 13.3 are the details of the air supply system (i.e., the wind 
tunnel) for the test section. Fresh air from the outside enters the large flow area 
“dryer,” which is plugged into a 4 ft x 8 ft window in the laboratory. The air then 
accelerates through a converging/diverging nozzle and enters the test section which 
is supported by two adjustable jack stands. After flowing through the test section, 
the spent heated air is vented to a large dump tank, which turns the air flow 90° and 
then feeds it into a converging duct. The duct is connected directly to the vacuum 
line as shown in Figure 13.3. 

Figure 13.4 gives an interior view of the dump tank and test section. Note that 
the laser beam enters from the left through a high power infrared “salt” window, and 
then fails upon the parabolic metal (uncooled) front surface mirror. This mirror is 
suspended by a three point adjustable mount, which can be aligned from the outside. 
The impulse plate is designed for easy removal and servicing. It is anticipated that 
numerous configurations will have to be tried before the best one can be identified. 
Note also that the test section is supplied with removable nozzle blocks so that the 
air inlet Mach number can be changed at will. 


13.4 Systems-Integration Issues 

Finally, additional systems-integration questions that were not addressed in this year’s 
final report will be explored. Several critical areas include LTD range safety during 
launch, manufacturing cost analyses, and possible high power laser attenuation by 
the engine exhaust plume. 
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Figure 13.1 Layout of RPI Laser Propulsion Laboratory in Ricketts Building. 

(showing proposed location of NRL CO 2 Laser and Experiment Test 
Sect ion) 



Figure 13.2 The impulse pressure plate will represent a 

full-scale, l/15th segment of the annular LTD engine. 
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Figure 13.3 Assembly drawing showing details of air supply system. 

(Side View) 



Figure 13. 4 Test Section, showing impulse plate and parabolic mirror. 
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Appendix A 


RESULTS OF COMPUTER GENERATED RAY 
TRACES 


This appendix summarizes the results of a computer ray trace study performed on 
the LTD primary mirror. Calculations were performed using the program “Beam 2” 
from Stellar Software; the graphical ray traces were rendered in the CATIA interactive 
computer graphics program. In the following tables, individual rays were numbered 
for easy reference to the traces; ray 1 is defined as incident along the centerline of the 
LTD craft, and rays are consecutively numbered proceeding outward from the vehicle 
centerline on either side. A total of 11 rays were traced for each side of the LTD in 
the extreme cases of 0.1 and 1.0 degrees vehicle incidence. Only the odd numbered 
rays were traced for the remaining cases. 

Tables 1 and 2 represent the perpendicular distances between incident rays and the 
ideal focus for various oblique angles of incidence. The corresponding rms blur is then 
calculated. Figure 1 illustrates graphically that the rms blur is linearly proportional 
to incident angle; the equation of the fitted line is also given. The remaining figures 
define the ray trace geometry, and present magnified views of the focal planes on 
either side of the vehicle for incident angles of 0.1, 0.2, 0.5, 0.7, and 1.0 degrees. All 
angles are measured counter-clockwise from the vehicle centerline; because the LTD 
is rotationally symmetric about this centerline, the same meridonal ray traces are 
valid for any cross section view of the craft. 

REFERENCES 

[1] CATIA Basic Three-Dimensional Design, User Manual, Vol. 1, prog. no. 5668- 
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Figure A.l: RMS blur vs. angle of incidence. 
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Figure A.2: Geometry for ray trace analysis, indicating the focus of parallel incident 
light rays. 
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0.1 degrees 


1 .0 degrees 


ray 

Left 

Right 

Left 

Right 

1 

0.079133 

0.085295 

0.826300 

0.840476 

2 

0.066513 

0.073328 

0.709123 

0.71 5295 

3 

0.055705 

0.061820 

0.595425 

0.601060 

4 

0.045750 

0.051470 

0.492471 

0.497744 

5 

0.036709 

0.042435 

0.402163 

0.405508 

6 

0.029108 

0.034038 

0.319032 

0.323672 

7 

0.022621 

0.027244 

0.248805 

0.253201 

8 

0.016225 

0.020690 

0.217250 

0.189286 

9 

0.011627 

0.016520 

0.139714 

0.144422 

If) 

0.007704 

0.012534 

0.102291 

0.105632 

11 

0.007278 

0.008480 

0.076231 

0.077365 

rms 





blur 

0.034397 

0.039441 

0.375346 

0.377605 


Table A.l: Perpendicular distance from ray intercepts to focal point at oblique inci- 
dence, and RMS blur (all values in cm) 


0.2 


0.5 


0.7 


1 0.162961 0.169128 

2 0.114840 0.121739 

3 0.077013 0.082779 

4 0.046950 0.052352 

5 0.025825 0.030647 

6 0.014932 0.016135 


0.414376 0.420642 

0.295197 0.301620 

0.196428 0.203807 

0.120776 0.127849 

0.068461 0.073368 

0.037894 0.039097 


0.582036 0.588282 

0.415029 0.421443 

0.278725 0.284490 

0.172485 0.177886 

0.046833 0.101790 

0.053057 0.054404 


rms 0.073753 0.078796 

blur 


0.188853 0.194394 0.258027 0.271466 


Table A. 2: Perpendicular distances and RMS blur for intermediate angles. 
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Figure B.l: Net coupling coefficient vs. Mach number and altitude using = 

0.75 for Mqo > 0.6 (the only change from baseline case). 
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Figure B.2: Net coupling coefficient vs. Mach number and altitude, using 
0.5 for Mqo > 0.6 (the only change from baseline case). 



Figure B.3: Net coupling coefficient vs. Mach number and altitude, using = 

0.25 for Mqo •> 0.6 (the only change from baseline case). 
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Figure B.4: Net coupling coefficient vs. Mach number and altitude, with the inlet 
gap reduced by 1 cm (the only change from baseline case). 
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Figure B.5: Net coupling coefficient vs. Mach number and altitude, with the inlet 
gap increased by 1 cm (the only change from baseline case). 
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Figure B.6: Net coupling coefficient vs. Mach number and altitude, with t^sd 
changed to 0.1 cm at t p (the only change from the baseline case). 
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Figure B.7: Net coupling coefficient vs. Mach number and altitude, with r LSD 
changed to 0.25 cm at t p (the only change from the baseline case). 
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Figure B.8: Net coupling coefficient vs. Mach number and altitude, with ti$d 

changed to 0.75 cm at t p (the only change from the baseline case). 
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Figure B.9: Net coupling coefficient vs. Mach number and altitude, with t^sd 

changed to 1 cm at t p (the only change from the baseline case). 
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Figure B.10: Net coupling coefficient vs. Mach number and altitude, with r LSD =0.1 
cm and $ LSD = 3 x 10 9 W/cm 2 at t p . 
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Figure B.ll: Net coupling coefficient -vs. Mach number and altitude, with t lsd = 0.25 
cm and $lsd = 3 x 10 9 W/cm 2 at t v . 
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Figure B.12: Net coupling coefficient vs. Mach number and altitude, with r LSD = 
0.5cm and $lsd = 3 x 10 9 W/cm 2 at t p . 
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Figure B.13: Net coupling coefficient vs. Mach number and altitude, withr^p = 1cm 
and $lsd = 3 x 10 9 W/cm 2 at t p . 
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Figure B.14: Net time-averaged thrust vs. Mach number and altitude, with t^sd — 
0.25 cm and $ = 3 x 10 9 W/cm 2 at t v . 
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Figure B.15: Pulse Repetition Frequency vs. Mach number and altitude, with tl S d 
= 0.25 cm and $ LSD = 3 x 10® W/cm 2 at t p . 
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Figure B.16: Pulse energy vs. Mach number and altitude, with t LS d = 0.25 cm and 
= 3 x 10 9 W/cm 2 at t p . 
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Figure B.17: Laser pulse time vs. Mach number and altitude, with tlsd = 0.25 
and $lsd = 3 x 10 9 W/cm 2 at t p . 
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Figure B.18: Time-averaged power vs. Mach number and altitude, with r LSD = 0. 
cm and $lsd = 3 x 10 9 W/cm 2 at t p . • 
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Figure B.19: Net coupling coefficient vs. Mach number and altitude, with t^sd 
0.1 cm and $ LSD = lO 10 W/cm 2 at t p . 
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Figure B.20: Net coupling coefficient vs. Mach number and altitude, with r LSD = 
0.2 cm and $lsd = 10 l ° W/cm 2 at t p . 
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Figure B.21: Net coupling coefficient vs. Mach number and altitude, with tlsd — 
0.5 cm and $lsd — 10 10 W/cm 2 at t v . 
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Figure B.22: Net coupling coefficient vs. Mach number and altitude, with r LSD = 1 
cm and § LSD = lO 10 W/cm 2 at t p . 
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Figure B.23: Laser pulse time vs. Mach number and altitude, with r LSD = 0.2 cm 
and $£,s£> = lO 10 W/cm 2 at t p . 
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Flight Mach number 

Figure B.24: Time-averaged power vs. Mach number and altitude, with t^sd = 0.2 
cm and $lsd = 10 10 W/cm 2 at t p . 



Figure B.25: Net time-averaged thrust vs. Mach number and altitude, with t^sd — 
0.2 cm and $lsd — 10 10 W/cm 2 at t p . 
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Figure B.26: Pulse Repitition Frequency vs. Mach number and altitude, with r isd 
= 0.2 cm and $lsd = 10 l ° W/cm 2 at t v . . 
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Figure B.27: Pulse energy vs. Mach number and altitude, with tlsd = 0-2 cm and 
$lsd = lO 10 W/cm 2 at t p . 
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